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Glossary 


a Lift curve slope (1/rad) 

a Speed of sound (ft/sec) 

b Wing span (ft) 

c Wing chord (ft) 

c Mean aerodynamic chord (ft) 

D Drag force (lbs) 

F Total force (lbs) 

F x Force component along the x-axis (lbs) 

F y Force component along the y-axis (lbs) 

F z Force component along the z-axis (lbs) 

Gravitational acceleration (ft/sec 2 ) 

Hinge moment (ft-lbs) 

Identity matrix 
Inertia matrix (slugs-ft 2 ) 

Moment of inertia about the x-axis (slugs-ft 2 ) 
Moment of inertia about the y-axis (slugs-ft 2 ) 
Moment of inertia about the z-axis (slugs-ft 2 ) 
Lift force (lbs) 

Rolling moment (ft-lbs) 

Vehicle mass (slugs) 

Mach number (dimensionless) 

Total moment (ft-lbs) 

Moment about the vehicle x-axis (ft-lbs) 
Moment about the vehicle y-axis (ft-lbs) 
Moment about the vehicle z-axis (ft-lbs) 
Pitching moment (ft-lbs) 

Yawing moment (ft-lbs) 
p Roll rate (rad/sec) 
q Pitch rate (rad/sec) 
q Dynamic pressure (psi) 

r Yaw rate (rad/sec) 

S Surface area (ft 2 ) 
t Time (sec) 

T Thrust force (lbs) 

u Velocity component along the x-axis (ft/sec) 
v Velocity component along the y-axis (ft/sec) 

V Aircraft velocity vector (ft/sec) 

V Velocity (ft/sec) 

Vh Horizontal tail volume (dimensionless) 
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Velocity component along the z-axis (ft/sec) 

Vehicle weight (lbs) 

Force along the x-axis (lbs) 
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Chapter 1 
Introduction 


The objective of this course is to develop fundamental understanding on the subject of stability, control and 
flight mechanics. Familiarities with the basic components in aerodynamics of wing and airfoil section arc 
expected; namely definition of lift, drag and moment of wing section, and physical parameters that govern 
these aerodynamic forces and moments suchas freestream velocity, density, Mach number, Reynold number, 
shape of the airfoil (camber, thickness, aerodynamic center), wing configuration (wing span, reference area, 
mean aerodynamic chord, taper ratio, sweep angle), angle of attack, dynamic pressure, etc- • -. It is not the 
intent of this course to provide all these relevant background materials, although we will define the relevant 
ones as we encounter them in our problem formulation. 

Starting from known forces and moments generated on a given wing, fuselage and tail configuration, we 
will develop airplane static and dynamic model to study its behavior under different flight regimes. Mass 
properties, wing, fuselage and tail configurations of the airplane are therefore assumed known and given 
a-priori. Concepts of static stability and dynamic stability arc introduced in this course. General equations 
of motion for a rigid-body airplane arc derived. Basic motions of the aircraft separated into longitudinal and 
lateral modes arc discussed in details. Effects of aerodynamic stability derivatives upon the behaviour of the 
perturbed equations of motions are studied. Flying qualities of the uncontrolled airplane can subsequently 
be assessed. Analysis of the airplane dynamic responses to initial changes in its basic motion variables (e.g. 
angle of attack, pitch attitude, roll angle, sideslip, etc- • •), tocontrol inputs and external gust inputs iscovered 
using Laplace transform techniques and time simulation. It is expected that students arc familiar somewhat 
with the use of the MATLAB software for analysis. 

There is an extensive number of references that cover the subject of aircraft stability, control and flight 
mechanics. Listed below arc some that provide good reference materials: 

1. John D. Anderson, Jr., Introduction to Flight, Third Edition, McGraw-Hill Book Company, 1989. 

2. Arthur W. Babister, Aircraft Dynamic Stability and Response, First Edition, Pergamon Press, 1980. 

3. John H. Blakelock, Automatic Control of Aircraft and Missiles, Second Edition, John Wiley & Sons, 
Inc., 1991. 

4. Bernard Etkin, Dynamics of Flight: Stability and Control,, Second Edition, John Wiley & Sons, Inc., 
1982. 

5. Barnes W. McCormick, Aerodynamics, Aeronautics, and Flight Mechanics, John Wiley & Sons, Inc., 
1979. 
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x-body axis 



6. Courtland, D. Perkins and Robert E. Hage, Airplane Performance, Stability and Control, John Wiley 
& Sons, Inc., 1949. 

7. Jan Roskam, Airplane Flight Dynamics and Automatic Controls, Roskam Aviation and Engineering 
Corporation, 1979. 

8. Edward Seckel, Stability and Control of Airplanes and Helicopters, Academic Press, 1964. 

9. David R. Hill and Clever B. Moler, Experiments in Computational Matrix Algebra, Random House, 
First Edition, 1988. 

As an introduction, let’s first examine the development of the following three basic equations corre¬ 
sponding to motion in the vertical plane (i.e., longitudinal set). They correspond to the lift, drag and moment 
equations respectively. 

1.1 Lift Equation 

Let’s consider the point mass system shown in Figure 1.1. Here we idealize the airplane as a lumped system 
with mass m andmomentof inertia aboutthe y-axis as I yy . Notethat theflight pathisalwaystangentialto the 
velocity vector V. The aerodynamic forces applied to the center of mass of the vehicle can be decomposed 
into the lift and drag components. The lift force is by definition perpendicular to the velocity vector V 
while the drag force is parallel to the velocity vector V and is pointed in the opposite direction. Gravity force 
mg and the engine thrust T constitute the remaining forces exerted on the vehicle. The pitching moment 
M v about the airplane center of gravity (eg) is mainly due to aerodynamic and propulsion forces and has no 
contribution from gravity. 

The equation of motion along the z-body axis is given by 

m (w — qu ) = XF, 


(1.1) 





1.2. DRAG EQUATION 
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where u is the component of velocity along the x-bcdy axis of the vehicle, w is the component along the 
z-body axis and q is the pitch angular velocity about the y-body axis. From Figure 1.1, we have 

FT- — F 7 (aerodynamics 

where 

F'z(aerodynamics ) T ^7.(propulsion ) 

Fz(gra vity ) = 

where W is the weight of the vehicle. Fet’s rewrite the velocity component w as w = Vsin a and thus, 

w = Vsin a + Vacosa 

= Va + Va (for small a) 

Generally, we notice that the product Va is much smaller than Va. Hence, equation (1.5) is simplified to 

w = Va (1.6) 

Combining equations (1.1), (1.2), (1.3), (1.4) and (1.6), we obtain the following equation in the z-body 
direction, 

m(Va -9V) = -L + (T - D)a + W (1.7) 

since q — 0 and u = Vcosa = V. Usually, the terms T — D and a are small and hence we can drop the 
product (T — D)a in the above equation. Thus, 


I T F z ( propulsion ) T F z(gra nily) 


( 1 . 2 ) 


= —Lcosa + (T — D)sin a 
= —L + {T — D)a (for small a) 

mgcos 6 

W (for small 9) 


(1.3) 

(1.4) 


mV (a — 9) = W — L 


( 1 . 8 ) 


Note that the flight path angle is defined as y = 9 — a, equation (1.8) can be rewritten as 

mV y — L — W (1.9) 

Thus, change in flight path occurs when L — W ^ 0 and the correspondingflight trajectory would be curved. 
For a constant flight path angle (i.e. y = y„ =constant), we must have y = 0 and L — W = 0. 

1.2 Drag Equation 

Again we refer to Figure 1.1, the equation of motion in the v-body direction is as follows, 

m(u + qw) = YjF x (1-10) 

since we are limited to motion in the vertical plane only. The force components in the v-body direction arc 
only consisted of Y.F X = F x(aerodynamics , + F x(propu i sion } + F x(gravity) . Each of these components can again 
be written in terms of the lift L, drag D, thrust T and gravity W forces according to Figure 1.1. Namely, 


F x (aerodynamics ) T ^ x ( propulsion ) 


Lsin a + (T — D)cosa 

La + (T — D) (for small a) 


( 1 . 11 ) 
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and 

Fx(gravity ) = —mgsin 6 ~ — W8 (for small 9) (1.12) 

Furthermore, the velocity component u can be rewritten as u = Vcosa. Hence, its time derivative becomes 


u = Vcosa — Vasina 

= V — Vaa (for small a) 


(1.13) 


Substituting equations (1.11), (1.12) and (1.13) into equation (1.10), we obtain the following 

mV + ma(V9 -Va) = T-D + (L- W)a - W(9 - a) (1.14) 


or, 

mV + Wy — T — D + a(L — W — mV y) (1.15) 

Using equation (1.9), equation (1.15) is simplified to 

mV + Wy — T — D (1.16) 

Thus from the above equation with excess thrust, i.e. (T — D) > 0, one can have different flight trajectories: 

• Positive flight path angle y > 0 with V — 0. This results in a steady (nonaccelerated) climb. 

• Positive acceleration V > 0 with y — 0. This corresponds to an accelerated straight and level flight. 

• Positive acceleration V > 0 and y > 0. The vehicle speed increases while climbing. 


1.3 Pitching Moment Equation 

Finally, we derive the pitching moment equation for the vehicle shown in Figure 1.1, 


lyy (i — My (aerodynamics ) T My ( propulsion ) (1-17) 

Notice that by definition, gravity would have no moment contribution to the pitching moment equation when 
it is taken about the vehicle center of gravity. Detailed description of the moments produced by aerodynamic 
and propulsive forces will be given later when we examine issues related to longitudinal static stability. It 
suffices to say that static longitudinal stability is predominantly governed by the behaviour of the pitching 
moment as the vehicle is perturbed from its equilibrium state. 



Chapter 2 

Linear Algebra and Matrices 


We deal with 3 classes of numbers: scalars, single numbers without association; vectors, one dimensional 
groupings of scalars (one column, several rows, or one row, several columns); and finally, matrices, which 
for us will be 2-dimensional (rows and columns). 

A vector can be either a row vector such as 


s — |.S| .v2 ■ ■ ■ s n ], 


or a column vector, such as 

■s'i 
*2 

S — 

Sn 


Column vectors are vastly more common. Implied with every vector is a basis (often a physical basis) to 
which each component refers. For instance, the position vector r given by 


r = xx + yy + zz 


is represented with respect to a cartesian basis [x, y, £]. Usually the shorthand [x, y, z ] is used. 

MATLAB will use both row and column vectors. However column vectors arc more often used in its 
functions. Note that in the example below, the first entry in boldface is what the user types, and the second 
corresponds to MATLAB’sresponse. 

• A = [1. 2. 3. 4.] 

A = 

1234 

• A = [l. 2. 3. 4.]’ 

A = 

1 

2 

3 

4 
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. A = [1;2;3;4] 

A = 

1 

2 

3 

4 

• A = [l 
2 
3 
4] 

A= 

1 

2 

3 

4 


A matrix can be thought of as a row of column vectors, or a column of row vectors. 


fll 


All 

012 • ' 

d\ m 

«2 

= 

«21 

an ■ ■ 

dim 

a n 


a n 1 

G)2 

dyim 


One can think of matrices as having a basis in the form of dyadic products of basis vectors, though that will 
be beyond the scope of this course. 

Entry of a matrix in MATLAB is fairly straightforward. It follows along the lines of a vector, but 
remember that the entries arc processed in row fashion: 

• A = [1 2 3 4;5 6 7 8;9 10 11 12] 

A= 

1234 
5 67 8 
9 10 11 12 

• A = [123 4 
5678 

9 10 11 12] 

A= 

1234 
5 67 8 
9 10 11 12 

• B = [A [0 1;2 3;4 5]; 6 5 4 3 2 1; 3 4 5 6 7 8; 1 3 5 7 9 11] 

B = 

123401 

567823 
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9 10 11 1245 
654321 
345678 
1 3579 11 

MATLAB also has facilities for creating simple matrices such as a matrix of zeros or the identity matrix. 
For example, the matrix function zeros(n.m) will create a zero matrix of dimension n by m, and eye(n) will 
create an identity matrix of dimension n. 


2.1 Operations 


Addition and multiplication not only need to be defined for within a certain class, but between classes. For 
example, multiplication of vectors and matrices by a scalar. 



a.si 


anil 

a«i2 • ■ 

0Ld\m 

a * v = 

as 2 

; a * M = 

cra2i 

aa 22 ■ ■ 

OCCl2m 


as n 


aa n i 

rr r/ 2 • ■ 

• • aa nm 


In MATFAB, however, you can add a scalar to every element in a vector or matrix without any special 
notation. 

Adding two vectors occurs on an element by element level. Implied in all of this is that the basis of the 
two vectors is the same: 


V + U = [l>i l>2 ••• V n \ + [U i U 2 ••• U n | 

= t Vi + Ml l>2 + U 2 • • • V n + u n 1 


Multiplication of two vectors is mostly defined in terms of the dot product. While much mathematical 
theory has been expounded on inner product spaces and such, the only item we need know here is the inner 
product of two vectors expressed in a Cartesian coordinate frame, 


v ■ u 


f i’l v 2 


v„] * 


u 1 
u 2 


L u n J 


V\U 1 + v 2 u 2 + • • • + v n u n 


Multiplying a vector by a matrix is equivalent to transforming the vector. In components, the product of 
a matrix and a vector is given by 



an 

a 12 • 

’ * &lm 


Vi 

A *v = 

«21 

a 22 ■ 

®2m 

* 

v 2 


Clnl 

Clnl 

’ * &nm 


Vm 


a\\V\ + a\ 2 v 2 4-h a\ m v m 

a 2 iV\ + a 22 v 2 + • • • + a 2m v m 

O n 1 r 1 1 + Cl n 2 V 2 + • • • + Cl/unVin 
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Of course, the number of columns of A must match the number of rows of v. 

Adding two matrices of the same dimensions would occur on an element by element basis, 



Oil 

012 

CL\m 


*11 

b\2 * * 

b\m 

A + B = 

021 

022 

’'' &2m 

+ 

*21 

b22 ‘ ‘ 

b2m 


O/il 

On2 

ttnm 


*nl 

bn2 

bnm 


Oil 

+ b\ i 

Ol2 + *12 ••• 

Cl\m 

b\ m 


— 

021 

+ *21 

022 + *22 • • • 

Cl2m 

+ b^m 



On 1 

+ b n i 

On 2 T" b n 2 • • • 

ttnm 

“I - b nm 



Multiplying two matrices can be thought of as a series of transformations on the column vectors of the 
multiplicand. The number of columns of the left matrix must be equal to the number of rows of the right 
matrix (left and right arc significant since multiplication is not commutative for matrices in general). 



Oil 

an ••• 

Cl\m 






A * B = 

021 

0 22 

& 2m 

* i b 2 ■ ■ 

■b P \ 





Onl 

On2 ’ ’ ’ 

Q-nm 







011*1 

1 + «12*21 + ■ ■ ■ 

+ a\ m b m \ 

011*12 

+ 012*22 + • ' 


Cl i m b m 2 

= 

021*1 

1 + 022 *21 + • • • 

+ Cl2mbm\ 

021*12 

+ 022*22 + ' 1 


&2mbm2 


Onl*ll + 0/72*21 

+... 

+ Cl n m bm 1 

On 1*12 

+ 0„2*22 + ' ' 


&nm bm 2 


Thestandardsymbols: +, *, and/willhandle all legal operationsbetween scalars, vectors, andmatrices 

in MATLAB without any further special notation. 


2.2 Matrix Functions 

The convenience of matrix notation is in the representation of a group of linear equations. For example, the 
following set of equations, 

anxi +a 12 x 2 H- +a Xn x n 

021X1 +U22X2 H- +02„X„ 

OnlXl To M 2 X 2 T +a nn X n 

can be represented compactly by the relation 

Ax = b. 

Note that we have the number of knowns b equal to the number of unknowns x here. How does one solve 
this? The most simplistic (and computationally efficient) method is to apply a succession of transformations 


= *i 

= b 2 

— bn 
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to the above system to eliminate values of A below the diagonal. Suppose that nn is nonzero then one can 
multiply the first equation by 0 , 21 /a\\ and subtracts it from the second equation. The first term in the second 
equation would be eliminated. 


anxi 

+ 

«12X2 

+ ■■ 

• • + 

n^n 

= b\ 

0 

+ 

(«22 - ai2«2lMl)x 2 

+ ■■ 

• • + 

(^2 n tt\ n Cl2\ / 

= bi b\a 2 \/a\ 


+ 


+ 

: + 


= : 

IX ] 

+ 

^n2^2 

+ 

• • + 

&nnX'n 

= b n 


Suppose that the procedure were repeated for all the other rows, resulting in the removal of the coefficient 
of x\ in these equations. The same procedure is now applied to all coefficients of X 2 for all rows below the 
second row. What one would eventually have is the upper triangular system. (In general, the a,j ’s and bf s 
are NOT the same as the original matrix entries a,;/ and bj, respectively). 


flllXl 

+ 

a 12 x 2 

+ 

•• + 

CL\ n X n 

= b 1 

0 

+ 

CI 22 X 2 

+ 

•• + 

&2 nXn 

= ^2-*l«2l/«ll 

0 

+ 

0 

+ 

•• + 


— 


+ 


+ 

: + 


= : 

0 

+ 

0 

+ 

•• + 

ClnnXn 

— b n 


All values of a,j below the diagonal arc zero. This matrixalso has the interesting propertythat the product of 
the diagonal terms is equal to the determinant. This substantiates the argument that it costs about as much to 
solve a lineal - system as it does to solve for a determinant. Note that one can now solve x n = b n /d nn . Once 
you have x n , you can substitute it into the next equation up and solve for x„_i. This continues until one gets 
to X] (or until some diagonal term dkk is zero). Note that if a diagonal term is zero, then the determinant is 
also zero and the system matrix A is called singular . 

The above method is often called the method of Gaussian elimination with back substitution. MATLAB 
implements a method very much similar to the above for solving a system of linear equations. You can, 
however, obtain an answer from MATLAB with very little effort by just typing the command 
x = A \ b 

In a similarvein to whatis notedabove, the determinant is computedwith the followingcommand syntax 
det(A). 

The inverse can also be computed in a method similar to that above. If one solves for a succession of 
vectors x; (i = 1, n), each one with 



" 1 ' 


" 0 " 

b\ = 

0 

II 

<N 

1 


0 


0 


thenthematrixcomposingofthecolumnvectors[ x\, ■■■ x n ] wo u 1 dbc t he m at r i x in vc rscof A. In MATLAB, 
computation of the matrix inverse is invoked by the command inv(A). 

Itisoftennecessarytocomputeadeterminantoraninverseinsomethingresemblingaclosedform(which 
will be seen in the calculation of eigenvalues). Thus, one introduces the expansion by minors method. A 
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minor M\j (A) of a matrix A is the determinant of the matrix A without its i th row and its j th column. 




Oil 

012 

013 

d\n 



021 

022 

023 

^2 n 

A 

= 

031 

032 

033 

*'* ^3 n 



On ] 

On 2 

On3 

&nn 




( 

021 

023 

n 

M 12 (A) 

= det 


031 

033 

''' ^3 n 




\ 

On 1 

On3 

Mnn 


A determinant is formed from expanding by minors along an arbitrary row i or column j of a matrix: 

n 

By Row/: det (A) = ^ a,-,- My (A) (—1) <+; 

./=i 

n 

By Column j: det (A) = ^ a,-,- M, ? - (A)(—1) I+J 

i=i 

Thus, what one has for an arbitrary matrix is a successionof expansions. First the matrix is broken down into 
a series of minors to determine the determinant, then these minors may need be broken down further into 
their minors to find their determinants, and so on until one gets to a 1 x 1 matrix. For a scalar (i.e., lxl) 
matrix, 

A — [on] =>■ det (A) = ci\\. 

It is also easy to evaluate the determinant of a 2 x 2 matrix by expanding along the first row. 


A - 


«ii fli2 

fl21 022 


det (A) = flnfl 22 — 012021- 


With a little more effort, one can get the determinant for the case of a 3 x 3 matrix. Here we expand along 
the first row. 


Oil 

012 

013 

021 

022 

023 

031 

032 

033 


det (A) = 


Ollfl22fl33 +Ol2fl23 a 31 +Ol3 a 2lfl32 
—011023032 — O 12 O 21 O 33 — O 13 O 22 O 31 


The inverse can also be found through expansion by minors. The recipe for this is a little more complicated 
than for the determinant. First one transposes the matrix, then each element Oy gets replaced by the term 
Mjj (A)(— I)' +/ , where M,y (A) is the minor of A at / and j. Finally, each resulting new element is divided 
by the determinant of the original matrix. For a 2 x 2 matrix, the inverse is 


A - 


011 012 
«21 «22 


inv(A) 


022 —012 


011022 — 012021 


—021 011 
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The inverse of a 3 x 3 matrix is somewhat more complicated. 


an 

a 12 

«13 


a 11 

a 21 

«31 

«21 

«22 

«23 

, A = 

fll2 

a 22 

a 32 

«31 

<?32 

«33 


an 

a 23 

a33 


A = in v(A) = 


1 


det (A) 


«22<333 — a 32^23 
«31«23 — «21«33 
«21«32 — «31 a 22 


«32«13 — «12«33 
«11<333 — fl , 31«13 
«31°12 _ a U a 12 


C112C123 ~ 013^22 
021013 — flllfl23 
O11O22 — 021^12 


Suppose that one applies this to the system 


an 

a \2 

an 


X\ 


£l 

«21 

a 22 

a 23 

* 

X 2 

= 

£2 

«31 

a 32 

«33 


. X 3 . 


. b 3 _ 


where the solution x is given by 


x = A 1 * b 


If one goes through the algebra, the result obtained from the Cramer’s Rule will be 


Xi = -—— 0ifl 2 2O33 + ^2032013 + £3012023 - £1032023 - £2012033 - £3022013) 

det (A) 


: = (etc...) 

This can be expressed more compactly as (with | A| being the shorthand notation for the determinant of A), 


£l 

an 

an 


«11 

£l 

«13 


an 

«12 

£l 

£2 

«22 

«23 


«21 

£2 

«23 


«21 

a 22 

£2 

£3 

«32 

«33 


<Ul 

£3 

a 33 


«31 

«32 

£3 


det {A) ’ 2 det (A) ' 3 det {A) 


2.3 Linear Ordinary Differential Equations 

Note that in this course we encounter mostly linear ordinary differential equations with constant coefficients. 
One can always find an integrating factor for the first-order equation 

x + ax = /, v(0~) = x a 

The term e at is an integrating factor for the above equation. Let’s multiply the above differential equation 
with the integrating factor e at and collect terms, we have the following 


Integrate both sides from 0 to time t. 


t 0 + f e aT f(x) 
Jo- 


dr 
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or, 

* — x 0 e~ at + [ e~ a(t ~ T) f(r) dr 
J o- 

Theabovesolutioncontainsusuallyahomogeneoussolution(frominitialconditions)andaparticularsolution 
(from the forcing function/). Note that solution for the particular paid involves a convolution integral. 

A second-order equation given by 

x + ax + bx = f, x(0~) = x 0 , x(0“) — v 0 

can be written as a system of two first-order equations as follows. Let x\ — x and X 2 = x, the above equation 
can be re-written as 



In fact, the procedure leading to the solution of the above scalar first-order differential equation can be used 
to derive the solution for a system of first-order differential equations. In the latter case, it would involve the 
matrix exponential (i.e., e At ). 

One could solve the above second-order system by first solving the homogeneous system whose solution 
is usually of the form x/,0) = e AI where A is the solution of the resulting characteristic equation. The 
particular solution from the non-homogenous paid can be done either by means of a trial substitution or 
variation of parameters. It turns out that the system of first-order differential equations is more amenable to 
use on the computer when expressed in a matrix state-space form. A formal way of modeling a dynamic 
system is by a set of state equations, 

x = Ax + Bit (State equations) 
y = Cx + Du (Output equations) 

Theinput u isthc particular forcing function driving the system (i.e., u = /). Refer back toour second-order 
system with both (x and x) as outputs. Then, in the above notation we have the following set of system 
matrices 



To generate a time series responseof a linear time-invaidant system inMATLAB, oneneeds to generatethese 
A, B, C, and D matrices. Suppose that one creates a vector of time points where the system outputs arc to 
be computed with the MATLAB command 

T =[0:0.1:10]; 

This vector contains time points from 0 to 10 in steps of 0.1. The forcing function f (often referred to as 
the control input ) is a vector whose entries correspond to the value of this function at those time points. To 
generate the system time responses for the system matrices A, B, C, and D as defined above, at the time 
points defined in the vector T to the forcing function /, we can issue the following MATLAB command 
Y = lsim(A,B,C,D,f,T); 

The vector Y has 2 columns (each corresponding to an output) and 101 rows (each row corresponding to a 
time point). 
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2.4 Laplace Transform 


Solving differential equations can also be done using the Laplace transform. We define 

POO 

L(f(t))= e~ st f (t)dt — f(s) 

J o- 


The function e at transforms as follows, 


L(e~ at 



dt 



~(s+a)t 

e -(s+a)t dt _ - 

s + a 


00 


o- 


1 

s + a 


Note that it is much easier to transform a function than to do its inverse transform which would involve 
from first principles the intricate details of complex variables and contour integration. However, it is much 
faster and easier for an engineer to use Table2.1. Every function used in Laplace transform work is assumed 
multiplied by a unit-step (Heaviside) function at t = 0. That is to say, the function is identically zero for 
t < 0 and equals to one for t > 0. 


f(t) 

f(s) 

m 

1 

1 (unit step at t=0) 

1 

t 

l 

s 2 

t 2 

1 



2 

c 

t „-l 

J_ 

(«—!)! 

s" 

e at 

1 

s—a 

te at 

1 

(s-cr) 2 

t n ~ l e al 

1 

(«-!)! 

(s-a) n 

1 - e~ at 

a 

s(s+cr ) 

sin {cot) 

CO 

s 2 ~hco 2 

cos (cot) 

s 

s 2 +co 2 

e at sin (cot) 

CO 

(s—cr) 2 +co 2 

e a, cos (cot) 

s—cr 

(s—cr) 2 +co 2 


Table 2.1: Laplace Transforms of Some Common Functions 


The Laplace transform is so attractive since 

d r°° d f 00 

L{ dh f{t)) = Jo- e ~ S, ^ f(t)dt = e ~ St + s J 0 _ e ~ S> f dt - -f(0~) + sf(s) 

Thus, one can transform a differential equation into a set of algebraic equations involving the transform 
variable .v. For example, Laplace transform of the first-order differential equation 


x + ax = f{t) 
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is 


sx (s) — x (0 )+ax (s) = f(s) 


Solving for the variable * (,v), we obtain the solution of the above differential equation in the Laplace domain 
as 


* 0 ) = 


x(o~) | m 

s + a s + a 


The first term of the sum corresponds to x(0~)e~ at , the second is not solvable until one specifies fit) (or 
f(s)). However, in general, the product of two Laplace transforms (in this case, /Tv ) and -^) is equivalent 
to the convolution of their time-based functions. This particular case is equal to what has been previously 
demonstrated, 


L- 1 


m i = r 

s + a J J o- 


e - a(t ~ T) f(x) dr 


Consider now the second-order differential equation 


x + ax + bx = / 


Applying Laplace transform to the above equation, we obtain 

.v 2 x(.v) — sx (0“) — x(0 _ ) + asx {s) — ax (0~) + bx{s) = f(s) 
and, solving for the solution x(s) 


x(0 ) + (s + a)x( 0 ) + f(s ) 

X(s) = - J- -—- 

s- + as + b 

The homogeneous parts of this have equivalent time-domain functions that depend on the relation between 
a and b. We distinguish three cases: 

• b — a 1 /4 > 0 

• b — a 1 /A < 0 

• b- a 1 / A = 0 

Case b — a 1 /A > 0: The denominator can be written into the form s 2 + 2o\v + a 2 + u> 2 which corresponds 

to the time-domain functions e~ at sin cot and e~ at coscot where o = —a/ 2 and co = y/b — a 2 /A. 

Solutions to the homogeneous problem (i.e., to initial conditions x(0)) can be obtained directly as 


Xhit) 


e -at/2 sin _ a 2/ 4t) 

*(° ) - /, - + 

y/b — a 2 /A 

x{0~)e~ at I 2 |co5'(yZr — a 2 /At) + — j====sin (yjb — a 2 /A t ) 


Use of Table 2.1 is not always possible with some more complicated forms. Usually one needs to break 
down a complicated polynomial fraction into simpler summands that arc of the forms given in Table 2.1. 
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Suppose that the forcing function fit) is a step input (applied at t=0) whose Laplace transform is simply 
1/s. We will derive the particular' solution as an illustration to the partial fraction expansion methodology. 
The particular solution to the non-homogenous problem is 

1 

; s(s 2 + as + b) 

The right-hand term can be decomposed into 

1 u vs + w 

s(s 2 + as+b ) s s 2 + as + b 

with unknowns u, v and w. Expanding and matching the numerator term, we have 

u (s 2 + as + b) + vs 2 + ws = 1 

Since the coefficients of s 2 and s must be zero, and that ub = 1, we have 

u = l/b , v = —l/b , w = —a/b 


The particular solution is given by 


1 

Xp(t ) = 7 


1— e at ! 2 I cos ijb — a 2 /A t) H- sinjjb — a 2 /At) 

1 V 2jb - a 2 /4 V 


Case b — a 2 /A < 0: In this case, we have two distinct real roots to the equation s 2 + as + b = 0. They are 
given by 

o i = —a/2 + J a 2 /A — b 


oi = —a/2 — -J a 2 /A — b 
Solution to the homogenous problem is simply 

x(0“) + (s +a)x(0~) 


Xh(s) = 


(s - (Ji)(s - ct 2 ) 


or. 


X h it) = 


x(0 ) + (cri + a)x(0 ) . x(0 ) + (cro + a)x(0 ) 


e ait + 


o\ — a 2 

Similarly, for the particular solution we have 


02 - cr i 


,<72t 


x P (s) = 


1 


s(s 2 + as + b) 


In partial fraction expansion 


U V w 

Xp(s) = -1-h 


s s — <7 1 s — a 2 

The unknowns u, v and w are determined from the following equation 


u(s 2 + as + b) + vs 2 — vsa 2 + ws 2 — wscr i = 1 
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For the coefficients of s 2 and s to vanish we must have u + v + w — 0 and ua — voj — wo\ = 0 with ub = 1. 
Thus, we have 


_ 1 _ 1 , _ 1 _ 1 

b o i<r 2 ’ cri(rri — «r 2 ) ’ ct 2 (ct 2 - CTl) 

Hence, the time-domain solution is 




/ 1 1 

V CT]CT 2 CTl (o' 1 — CT 2 ) 


+ 


1 


02(02 — o\) 


H{t) 


where Hit) coiTesponds to the Heaviside (step) function at t = 0. 
Case a 1 / 4 — b: This case is similar to the previous case where 


a 

02 = 01 = - 2 

Solution to the homogenous problem is simply 


x h (s) = 


x(0 ) + (s + a)x (0 ) 

(S - CTl ) 2 


or. 


x h (t) = x(0 ~)te ait + x(0 - )(l - CTiOe" 1 * 
Similarly, for the particular solution we have 

1 u v w 


xJs) = 


u 

= - + 


+ 


sG 2 + as + b) s s — cti (s — cti ) 2 


For the coefficients of s 2 and s to vanish we must have u + v = 0 and —2 cti u — o\v + w =0 with u = 1 jb. 
Hence, 

1 1 1 


w = — 
o-i 


Or, in the time domain 


x p (t) =(\- ~\e ait + -te a A H(t ) 

W CT i / 

InMATLAB.timeresponsesofasystemcanbeobtainedfromtheirLaplacetransformsdirectly. Response 
of the outputs to an input U defined over the time points T can be obtained using the command 
Y = lsim(NUM,DEN,U,T); 

where the arguments NUM and DEN arc arrays containing coefficients of the numerator and denominator 
polynomials in s arranged in descending powers of .v. 


2.5 Stability 

Dynamicstability ischaracterizedbythe responseofasystem tononzeroinitialconditions. Initialconditions 
arc equivalent to an impulsive forcing function (i.e. Dirac delta function u ( t) = 8(t)). 

For a first-order system x + ox =0 and x (0 _ ) = x a , we have the homogeneous solution x it) = x Q e~ at . 
It is simple to imagine that for a > 0, the response x(t) to initial conditions x 0 would tend toward zero (thus 
stable) as t —> 00 . On the other hand, if a < 0 the response would tend to blow up. 

For a second-order system x + ax + fix = 0, we have solutions of the form 
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• When a 2 /4 < b , 


*(7) = e at [usin cot + vcoscot ] 


• When « 2 /4 > 

x(f) = Mg' 71 ' + ve mJ 

• When « 2 /4 = b, 

x(t) — Me 0 "' + vte al 

In any case, the argument a in the exponential function e at term must be less than or equal to zero. In the 
case of a 2 /4 > b, both terms oq and oq must be less than or equal to zero. Otherwise, the solution will blow 
up as t —»■ oo. For o equal to zero, then one has a neutrally stable system. 


2.6 Example 


Consider the following ordinary differential equation 


d 3 y(t) d 2 y(t) 

dt 3 dt 2 


+ 17 


dy(t) 

dt 


+ 13 y{t) = 13n (t) 


(2.1) 


with initial conditions y(0") = 7, y(0“) = 0 and y(0~) = 0. Solve for the time response y(t) when 
u(t ) — 8(t) (impulse function or Dirac delta function). 

We can solve the problem using three methods: 


• Laplace method 

• Time-domain method involving the matrix exponential 

• Numerical integration method (via MATLAB) 


2.6.1 Laplace method 

Taking the Laplace transform on the differential equation, we obtain 

[Vy(s) - s 2 y((T) - sy(0~) - y(0")] + 5 [s 2 y(s) - ^(0") - y(0“)] + 

17 [sy (s) - y (0“)] + 13y(s) = 13 u(s) 

With y (0”) = 7, y(0~) = y(0~) = 0, equation 2.2 becomes 

[.v 3 + 5^ 2 + 17s + 13] y(s) - [s 2 + 5^ + 17 ] y(0“) = 13 u(s) 

or, solving for v (.v) we have 


5 2 + 55 + 17 

5 3 + 55 2 + 175 + 13 } ’ + 


13 


u(s ) 


( 2 . 2 ) 


(2.3) 


y(s) = 


5 3 + 5s 2 + 175 + 13 


(2.4) 
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For a unit impulse input u(t) = 8(t ) or u(s) = 1 and y(0 ) = 7, equation 2.4 becomes 

y(s) = 5-2 + + 17 - 7 q—-^13- 1 

7 ■ 5 3 + 5s 2 + 17^ + 13 5 3 + 55 2 + 175 + 13 

75 2 + 355 + 132 
(5 + l)[(5+2) 2 + 3 2 ] 

Performing the partial fraction expansion on equation 2.5, we have 

R\ R 2 R 2 

y(s) = - 1 --- 1 --- 

^ 5 + 1 5 + 2 — 73 5+2 + 73 


Is 2 + 355 + 132 

R\ = - r -=-=- = 10.4 

[(5 + 2) 2 + 3 2 ] i= _j 

R 2 = 7>i + 3S>+132 = —1.7 — jO.f 

(5 + 1) (5 + 2 + 73) S= _ 2+J3 


Is 1 + 355 + 132 


= -1.7 + 7 0.6 = Ri 


(5 + 1) (5 + 2 — 7'3) \ s __ 2 _ j3 ■ - 

where R 2 is the complex conjugate of R 2 . Taking the inverse Laplace tranform on equation 2.6, we have 




y(t ) = R x e + R 2 e ( - 2+j3)l + R 2 e ( - 2 ~ j3)t (2.7) 

ice e a+jb = e a {cosb + jsinb ), we have 

y(t ) = 10.4e —t + e~^ 1-7 — j0.6)(cos3t + 7 'sin 3 1) + (—1.7 + j0.6)(cos3t — jsin 3 1) (2.8) 

y(t) = +2e~^(—l.lcos3t + 0.6sin 3t) (2.9) 


2.6.2 Time-domain method 

Let’s define x\(t) = y(t), x 2 (t) = y(t) and x 2 (t) = y(t), thencleai'ly 

x\(t) = x 2 {t) = y(t) (2.10) 

x 2 (t) = x 3 (t) = y(t) (2.11) 

and equation 2.1 can be re-written as 

*3 (t) + 5^3 (0 + 17x7 + 13xi(t) = 13 u(t) (2.12) 

Combining equations 2.10-2.12 into a set of three first-order differential equations which can be expressed 
in a matrix equation, 

”xi(i) _ ”0 1 0 ” ”xi(0 _ 0 

x 2 (t) =001 x 2 (t) + 0 u(t ) (2.13) 

_x 3 (t)J L —13 -17 - 5 J Ljc 3 (oJ L 13 _ 
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with initial conditions xi (0 ) = y (0 ), X 2(0 ) = v (0 ) and * 3(0 ) — y (0 ). 
Equation 2.13 can be written in a concise form as 

| x(t) = Ax(t) + Bu(t) 

I x( 0 ~) = x 0 


(2.14) 


where 


x(t) = 


~Xl (f) _ 
X2 (t) 
-X3 (t)_ 


A = 


0 1 
0 0 


-13 -17 


0 ' 
1 

-5. 


B = 


~ 0 " 
0 

.13. 


x 0 


~y(.o -)" 

j(O-) 

.y(o-). 


Solution of equation 2.13 is obtained using the method of linear supeiposition and convolution (Duhamel) 
integral. Namely, 

x(t) — e At x 0 + [ e A( ' t ~ z) Bu (r)dr (2.15) 

J 0 - 

Since u(t ) = 8(t ) equation 2.15 becomes 

x(t) = e At x 0 + f e A( ' t ~ T) B8(r)dT (2.16) 

do¬ 


or 


Since v 0 


0 

. 0 . 


, we have 


x(t) = e At x 0 + e At B = e At (x Q + 6 ) 


or 



1 

O O 

1 0 ■ 

0 1 

1 / 

r7i 


■ 0 ■ 

x(t) = e 

L —13 

-17 —5_ 

( 

1 

0 0 
_ 1 

+ 

0 

.13. 


x(f) = e L 


0 1 0 
0 0 1 
-13 -17 -5 


t r 7 n 
0 
13 


(2.17) 


(2.18) 


(2.19) 
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2.6.3 Numerical integration method (via MATLAB) 

In MATLAB, the command lsim would perform numerical integration on the linear system 

f x(t) = Ax(t) + Bu(t) 

[ y(t) = Cx(t) + Du{t ) ^ ; 

for a given initial condition x(0 _ ) = x 0 and an input function u(t ) defined in the time interval 0 < t < t max . 

More precisely, we can use the following set of MATLAB codes to perform the time responses of y{t) 
for the linear system described in equation 2.20. 


t=[0:.1:tmax]; 

u= "some function of t" (e.g., u=ones(t) for a step input) 
y=lsim (A, B, C, D, u, t, xo) 
plot(t,y) 

Below is a complete listing of the m-file for this example problem. 

%Using the Laplace method 
t=[0:.1:5]; 

y=10.4*exp(—t)+2*exp (-2*t) .* (-1.7*cos(3*t)+0.6*sin (3*t)); 

%Using the state-space method 
A= [0, 1, 0; 0, 0, 1; -13, -17, -5] ; 

B=[0;0;13]; 

C=[1,0,0]; 

D=0; 

xo=[7;0;0] ; 

[n, m]=size (t); 

%Solution from the exponential matrix 
yexp= [ ]; 

for i = l:length (t) 
ti=t (i); 

ytemp= expm(A*ti)*(xo+B); 

yexp=[yexp;ytemp']; 

end 

%MATLAB command for solving time responses 

%For an impulse input the new initial condition becomes x(o+)=x (0-)+b 
% and the input u(t)=0 for t>0+ 
xoplus=xo+B; 
u=zeros(n,m); 

ysim=lsim(A,B,C,D,u,t,xoplus); 

%Plot the responses for comparison 

plot(t,y,'o',t,yexp) 

grid 

xlabel('Time (sec)') 
ylabel('y(t)') 

title('Exponential matrix method') 
pause 

%Plot the responses for comparison 

plot(t,y,'o',t,ysim) 

grid 
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xlabe 

ylabe 

title 


Exponential matrix method 
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Chapter 3 

Principles of Static and Dynamic Stability 


In most design situation, static and dynamic stability analysis plays a significant role in the determination of 
the final airplane design configuration. The decision is based according to the requirements defined in FAR 
Part 23 which states that 

“the airplanemust besafely controllableand maneuverableduring — (l)take off; (2)climb; 
(3)levelflight; (4)dive; and(5)landing(poweron, off)(with wingflapsextendedandretracted)” 

Stability of such a vehicle is also a major consideration in selecting a particular design configuration. The 
airplane must be longitudinally, directionally and laterally stable for airworthiness and minimal pilot work¬ 
load. If the airplane turns out to have undesireable flying qualities, then some of these requirements must 
subsequently be met by the use of stability augmentation systems. This requires careful design of a control 
system that feedbacks sensed aircraft motion variables to the appropriate control surfaces (e.g. elevator, 
aileron and rudder). The topic of feedback synthesis of flight control systems for stability augmentation and 
autopilot designs is the subject of AA-517 and a continuation in AA-518. In the present course, we will 
only examine the fundamental behaviour of flight vehicle and its inherent flight characteristics without the 
influence of artificial feedback control. 

The general notion of stability refers to the tendency of the vehicle to return to its original state of 
equilibrium (e.g., trim point) when disturbed. There are basically two types of stability: 

• Static stability refers to the tendency of an airplane under static conditions to return to its trimmed 
condition. Clearly we assume that there exists an equilibrium point about which static stability is 
investigated. The evaluation of static stability involves purely static (i.e. steady-state) equations from 
force and moment balance applied to a vehicle disturbed from its equilibrium. Conditions for stability 
arc governedby the directionof the forces andmoments thatwill restorethevehicletothe originaltrim 
states. Figure 3.1 shows the three possible cases of static stability. Clearly, from these illustrations, we 
determine stability from the direction of the restoring force. In Figure 3.1 (a) component of the gravity 
force tangential to the surface will bring the ball back to its original equilibrium point. However, in 
static stability analysisthereisnomentionon how andwhen theballwillreturntoitsequilibriumpoint. 
For example, without the benefit of friction, the ball will oscillate back and forth about the equilibrium 
pointandthereforewillnever reachtheequilibriumstate. Totreatthisproblemconcerningthedynamic 
behaviour of this ball rolling on a curved surface, one will need to first develop its equation of motion 
and then analyze the stability of its motion when released from a perturbed position. Stability of the 
ball in motion is then determined by the phenomena of dynamic stability. 


23 
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O_0_ 

(c) Neutrally Stable 

Figure 3.1: Three Possible Cases of Static Stability 


• Dynamic stability isgovernedbythefactthevehiclewillreturntoitsoriginalequilibriumconditionafter 
some interval of time. As discussed in the previous section, analysis of dynamic stability would entail 
a complete modeling of the vehicle dynamics and its responses when perturbed from the equilibrium 
state. Figure3.2 shows typicalresponsesof a dynamicallystable, unstable andneutrally stablesystem. 
It is important to observe from the above examples that a dynamically stable airplane must always be 
statically stable. On the other hand, a statically stable airplane is not necessary dynamically stable. 
Detailed study of dynamic stability of a flight vehicle will be performed following the development of 
the general equations of motion of a rigid-body airplane in Chapter 7. 
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(a) Dynamically Stable 



(b) Dynamically Neutrally Stable 




(c) Dynamically Unstable 


Figure 3.2: Three Possible Cases of Dynamic Stability 
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Chapter 4 

Static Longitudinal Stability 


A study of airplane stability andcontrol is primarily focused on momentsabout the airplane center of gravity. 
A balanced (i.e., trimmed) airplane will have zero moment about its center of gravity. The total moment 
coefficient about the center of gravity is defined as 


C M cg 


Meg 

qSc 


(4.1) 


where S is the wing planform area, c is the mean aerodynamic chord and is the dynamic pressure 
corresponding to the freestream velocity V^. There are numerous places where moments can be generated 
in an airplane (Figure 4.1) such as moments contributed by the wing, the fuselage, the engine propulsion, 
the controls (e.g., elevator, aileron, rudder, canard, etc...) and the vertical and horizontal tail surfaces. Note 
that the gravity force does not contribute any moment to the airplane since it is, by definition, applied at 
the center of gravity. The aerodynamic center for the wing is defined as the point about which the moment 
Mac (or its moment coefficient Cm,oc ) is independent of the angle of attack. This point is convenient for the 
derivation of the moment equation since it isolates out the part that is independent of the angle of attack. 


4.1 Notations and Sign Conventions 

Hereweintroducetheconmionlyusednotationsfordisplacements,velocities,forcesandmomentsinstability, 
control and flight mechanics. The origin of the axis system defined by the x, y, ^-coordinates is assumed 
fixed to the center of gravity of the airplane (see Figure 4.2). It will move and rotate with the aircraft. The x 
displacement has a positive forward direction, the y displacement has a positive direction to the right-wing 
directionwhile the z displacements pointedpositively downward. Therespective componentsofthe aircraft 
velocity V in the x, y and z directions are (u, v, w) respectively. The total force F applied to the airplane 
has components (X, Y, Z) while the respective moment components are (L, M, N). Note that all the forces 
and moments are assumed to apply at the center of gravity. 

We will examine a simple airplane configuration in our analysis of longitudinal static stability. The basic 
airplane consists simply of a wing and tail configuration only. This simple configuration will illustrate well 
the basic fundamentals in stability and control analysis. 


4.2 Stick-Fixed Stability 

The forces and moments of a wing-tail configuration is shown in Figure 4.3. Without loss of generality, the 
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Ligure 4.3: Lorces and Moments Applied to a Wing-Tail Configuration 
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horizontal axis is assumed to coincide with the zero-lift line of the wing. Relative to this reference line, the 
tail is shown to have a positive incidence angle. Note that in our development, we adopt the same standard 
convention for all angle definitions (i.e. according to the right-hand rule). The angle of attack of the wing 
with respect to the zero-lift line is defined as a w . At the tail, the angle of attack is reduced by an angle e due 
to the downwash at the wing. The airplane is in equilibrium when sums of all the forces and moments about 
the center of gravity are zero. 

In the longitudinal axis, we have 


SF, = W 
= 0 


(L w cosa w + D w sina w ) — [L t cos(a w — e) + D t sin (a w — e)] 


(4.2) 


and 

EM 


= M a 
Ma 

= 0 


h a 

[L t cos (a w - e) + D t sin (a w - e)]l t + [ L,sin (a v 


- D w cosa w )z w + 

e) - D,cos (a w - e)]z, 


(4.3) 


To simplify our analysis, one can usually assume that the angle of attack a w is small and use the following 
approximations for cosa w = 1 and sina w = a w where a w is in radians. Then equations (4.2) and (4.3) 
become 

W = (L w + D w a w ) + \L, + D t (a w — e)] (4.4) 


and 


M a c,w T (L w -f- D w a w )(h h a c,w)c 4" (L w oiui D w )zw~\~ 

M ac j — [L t + D t (a w — e)]l t + [L t (a w — e) — D t ]zt = 0 


(4.5) 


Let’sintroducethefollowingdefinitionsfornon-dimensionalforcesandmomentsatthewingandtailsurfaces, 


M t 

L t 


= qS 


dC I n 

da 


— qSa vjOLyj 

— qScC M ac i v 
= q,S t C Ll 

= <LS t d -^a t 
= q t S t ^-(it+a w -e) 
= q t S t a,{i, + a w - e) 
M a c,t — qtStC t C Macl 


(4.6) 


Furthermore, we define the following 


Vt = 


<7r 


(4.7) 


~ e ° + da aw 
= "f ^a&w 


(4.8) 
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where e„ is the downwash angle when the wing is at zero lift. Both e 0 and e a are obtained from wind 
tunnel data. And the variable rj t is simply the ratio of dynamic pressure at the tail to the freestream dynamic 
pressure; this can be greater or less than unity depending on whether the tail is in the wake of the propulsion 
system or not. 

Lurther simplification can be obtained using the fact that D w a w <<C L w , D t {a w — e) <?C L t , z w = 0 and 
z t = 0. Equations (4.4) and (4.5) simplify to 

r — w 

^ L qS 

= a w a w + rj t ^a t {i t + a w - e) 

= Cl„ + CL a a w 

= CL a (oiw ®-w,zerolift ) 

where 

Cl„ = m^a t (i t - e Q ) 

Ch a = T 4 1 v T 0 €a) 


(4.9) 


(4.10) 


and 

C M cg 


where 


CMac.w 4 ” tl w (X u l(H h ac ,w) T ^ qSc C Mac j C \sc t ) T CMafus 

C 4 v 4 " C M acl tjfVqjClti.it to) + {(h ^ac,w)^w ^tCjjClfii t<x) T CM a f us }&uv 

Cm,, + Cm u oi w 


Cm„ — CM acw + q-CM acJ — rj t Vjja t (i t ~ t 0 ) 

Cm u = (h h a c,w)tlw 4 / CHCli ( 1 C a ) + CM a f lls 

- dC L '-A* 


(4.11) 


(4.12) 


CMafus identifies the contribution of the fuselage to the pitchingmoment (it is generally negligeable C M afm ~ 
0), and Vh — ^ isthchorizontaltailvolumccocfficicnt. Recall thatinequilibrium.wemusthave Cm cs = 0. 
Refening to Ligure 4.4, one can see that there arc two possible cases for an equilibrium to exist; namely. 




4.2. STICK-FIXED STABILITY 


31 


1. Cm„ > Oand Cm„ < 0: This case corresponds to a statically stable equilibrium point since for any 
small change in the angle of attack, a restoring momentis generated to bring it back to the equilibrium. 


2. Cm b < Oand C > 0: This case corresponds to a statically unstable equilibrium point since the 
moment created due to any change in angle of attack will tend to increase it further. 


There exists a location of the center of gravity, i.e. when li = h n , where the coefficient Cm„ = 0. Recall that 
Vh = (Tail volume coefficient) and /, = (h ac ,t — h)c, then equation (4.12) becomes, with h substituted 


by h n . 


( h n h acw )a w ^ltihac,t 



Of( 1 — e a) + CM afus — 0 


(4.13) 


or 


St 

[Uu> + I h n = h acw a w 


T ^£h(l £a)hac,t C M a f us 


(4.14) 


Let’s examine the total lift on the wing-tail configuration, it is given by 


L — L w -(- Lf 

— qSa w a w + q,S t a t (i t + a w - e) 

— qSci w ot w T q t S t a t {i t 

or the total lift coefficient Cl is 

C L = Vtja t (i t - e a ) + [a w + q,^a t (l - e a )]a w 
= Cl o + CL a Ot w 

Thus, the combined lift curve slope is 

$ 

C La = a w + q t ^a t ( 1 - e a ) 

Lrom the above definition of Cl„, equation (4.14) is simplified to the following 


CL a h n — h a c,w®w + (Cl„ U,/i )haci C M a f us 

or the neutral point h n is given by 

ha C ,w + — l]hac,t C Mafm 


hn — 


C La 


Cl„ 


h — h — Uw (h — h i — CMa f m 
,L n — rL ac,t V L ac,t rL ac,wJ C Lci 


Lrom the above, it can be easily shown that 


Cm b = CL a (h — h n ) 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 


Note that C/.„ > 0, thus C^ a < 0 if (h — h n ) < 0 or, the center of gravity must be ahead of the neutral point. 
The other condition Cm„ >0, where C /W „ is defined in equation (4.12), will be satisfied if the tail incidence 
angle i, is negative. The quantity (h n — h) is called the static margin. It represents the distance (expressed as 
a fraction of the mean aerodynamic chord) that the center of gravity is ahead of the neutral point. Roughly, 
a desireable static margin of at least 5% is recommended. Lor airplane with relaxed static stability, the static 
margin is negative. A stability augmentation system (SAS) is needed to fly these vehicle. 
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Example 1 

Given a light airplane with the following design parameters, 

• Wing area S w = 160.0 ft 2 . Wing span b w = 30 ft, 

• Horizontal tail area S t = 24.4/t 2 , Tail span b t = 10 ft, 

• h acJ = 2.78 

• Wing with 652 — 415 type airfoil, Cm„ c = —0.07 , h ac w = 0.27, 

• rj t = 1 and e a = 0.447. 


The lift curve slopes at the wing and tail arc obtained from the following empirical equation, 

A 

«3 n — 0 ">n - 

A + [2(A + 4)/(A + 2)] 


(4.21) 


where A — b 2 /S is the aspect ratio of the surface and no sweep. Thus, with a 2 d = 2n per radians = 0.106 
per degrees, we have 

a w = 0.073 lper degrees 
a, — 0.0642per degrees 

The total lift curve slope according to equation (4.17) is C/.„ = 0.0785 per degrees, and the neutral point is 
at h„ = 0.443. Then 

C Ma = 0.0785 (h - 0.443) (4.23) 


Calculation of C Mac , Aerodynamic Center Location and Mean Aerodynamic Chord (mac) for 
a Finite Wing 

Consider a finite wing shown in Ligure 4.5. The locus of the section aerodynamic centers defines the swept 
back angle A. The pitching moment about a line through the point A and normal to the chord line is given 
by 

fbl 2 rb /2 

MA = q / c C mac dy — q / cC/vtanAdy (4.24) 

J-b/2 ‘ " J-b/2 

Then if Xa is the distance of the aerodynamic center behind the point A, then 

M ac = M a + LX a (4.25) 

or 

C Mac = C Ma + C l ^- (4.26) 

c 

Differentiating C m uc with respect to a and using the definition of an aerodynamic center yields 

dC Ms Xa 

0 - —pA + -LC La {A.21) 

da c 

Substituting equation (4.24) into the above equation and using the fact that 
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Figure 4.5: Calculation of Wing Aerodynamic Center 


we obtain from equation (4.27), 


X A = 


- -l. r . 

C La Sj-b/2 


cCi a ytanAdy 


(4.29) 


If we assume that C/ ff is constant across the wing span, then we have 

X A = 


fo /2 cydy 

s/2 

Y _-C±_ 

A y c La 

where y is the spanwise distance from the centerline out to the centroid of the half-wing area. As a special 
case, for a linearly tapered wing, equation (4.31) becomes 


Cl 

-tan A 

(4.30) 

Cl„ 

tan A 

(4.31) 


v (1 + 2A) C la b f 

Xa =-tan A 

(1 + X) C La 6 


(4.32) 


where X — — is the wing taper ratio. The mean aerodynamic chord c of a finite wing is defined as the chord 
length that, when multiplied by the wing area S, the dynamic pressure q, and an average Cm,,,. , gives the total 
moment about the wing’s aerodynamic center. Namely, 


M ac — qS cC 


(4.33) 


Combining the above equation with equation (4.24), we have 


f b/2 9 

qScCM ac —q / c~C„ lac dv (4.34) 

J-b/2 

Thus, if the wing is straight and has constant airfoil cross section (i.e. C,„ ac is constant across the wing span), 
then we have c = c. However, if c is not constant (e.g in a tapered wing) and we assume that Cm uc = C, Uac 
and Ci are constant across the wing span, then the mean aerodynamic chord c is simply, 


c = 



(4.35) 



34 


CHAPTER 4. STATIC LONGITUDINAL STABILITY 


Hinge 



(a) Stabilizer-Elevator Configuration 
Hinge 



(b) Stabilator Configuration 
Ligure 4.6: Horizontal Tail Configurations 


This integral definition of c is used for any planform. As an example, for a linear tapered wing, we have 


2 c 0 1 T X T 
~3 1 +A. 


(4.36) 


where c a is the midspan chord (Ligure 4.5). 


4.3 Stick-Free Stability 

Wehaveseenintheprevioussectionthekeyelementsinstaticstabilityanalysisforastick-fixedconfiguration. 

It was assumed that the position of the tail or elevator surface has been fixed by the pilot holding onto the 
control stick, i.e. to holdthe surfaceintrimmedposition the pilotmust exert aconstant force dueto anonzero 
moment at the elevator hinge. This may not be desireable for long duration flight. Of course, nowadays for 
highperformanceandlarge-sizeairplane, theproblemisalleviatedwith theuseofpowerassistedcontrolsand 
seldom there arc unassisted control linkages between the pilot controls and the respective control surfaces. 

Nevertheless, it would still be necessary for small-size airplanes to investigate the issue of stick-free 
stability. It turns out that the effect of freeing the control surface amounts to a reduction in static stability 
in a certain configuration (e.g. stabilizer-elevator). Let’s examine the two basic configurations of horizontal 
tail surfaces: stabilizer-elevator and stabilator as shown in Ligure 4.6 . 


Horizontal Stabilizer-Elevator Configuration 

Let’s consider the moment H e about the hinge line of the elevator and the corresponding elevator hinge 
moment coefficient C/ !f , defined as, 

C he = — He (4.37) 

l/lpV^SeCe 
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The elevator hinge moment coefficient C/ !f is found to be a function of the tail angle of attack a, and of the 
elevator deflection 8 e . As an approximation, one can write 


C he = 


dC he 

da, 


a, + 


dC he 
38 e ' 


(4.38) 


where dC^Jda, and dC],Jd8 e are assumed constant and determined empirically (i.e. they vary with the 
configurationoftheplanformofthestabilizer-elevator). Withtheconventionthatapositiveelevatordeflection 
is down , these derivative coefficients are usually negative thus producing a negative hinge moment for any 
positive change in either a, or 8 e . 

Clearly, the free elevator will reach an equilibrium position when its hinge moment is zero for any tail 
angle of attack a,. Let’s denote this angle as 8 efree which is determined by setting C/ !f , equal to zero. 


C he = 0 = 


dC h , 

da, 


a, + 


3 c h e „ 

Op r 

ds e £free 


(4.39) 


This equation allows us to solve for 8 efree in terms of the angle of attack at the tail a ,. The tail lift coefficient 
derived from equation (4.6) is then modified to include the effect of a free elevator as follows, 


dC L , 

C L , = a,a, + ——8 e 

OOp 

However, since for a stick-free case, 8 e = 8 efree , equation (4.40) becomes 

dChp 

„ _ dC L , Hit „ 

Cl, — a,a, ac, lf a? 

e ~dS7 


C L , = F e a,a, 


where 


F, = 1 


1 3C L, da, 


= 1 — r 


3 C he 
da, 


a, d8 e dCh * ac '«’ 

dS e d s e 


where r = t-A is the elevator effectiveness (see Figure 5-33 on page 250 of Perkins & Hage) in 


(4.40) 


(4.41) 

(4.42) 

(4.43) 


3Clj _ dCL.t da, 
dS e da, d8 e 


(4.44) 


and a, is the lift-curve slope of the tail. The variable F,. is called the free elevator factor, and it is usually less 
than unity. Stability analysis for the stick-free case proceeds exactly as in the stick-fixed case. The results 
arc obtained simply by substituting a, in equation (4.17) by F e a t . Namely, 


(4.45) 


The lift curve slope for a stick-free case is always less than that of a stick-fixed case. From the above result 
for CL a free , the neutral point h n is given by 


htl f ree 


h ac ,w + - 1 ]hac,t 

CL a free 

a w 



(4.46) 
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or 


free ~ h 


ac,t „ (hac,t 

C L a free 


hac,w^ 



(4.47) 


Notice that since C^ u j ree < C[ u f ixec i , we deduce that the neutral point for a stick-free case is ahead of the 
neutral point of a stick-fixed case (i.e. h nfree < h nfjxed ); hence the stick-free case is less statically stable than 
the stick-fixed case for a given center of gravity position. 

Lrom the above, it can be easily shown that 


M a free — C L a free ^n/ rec ) 


(4.48) 


Example 2 


Using results from Example 1 and assuming that C/, f = —0.31a, — 0.68A (where a t and 8 e are in radians) 
with an elevator control effectiveness of 1.616 per radians. Then 


1 3 C Lt dC h Jda, 
at dS e dCh e /dS e 

i i.6i6per rad (-0.31) _ 0 sn 

0.0642per deg 57.3 deg/rad (—0.68) 


(4.49) 


The lift curve slope C[ u is given by 


t--L a free — tl w + TJ t ^ L l? fl',(l 

= 0.073 lper deg + 1 x y^0.80 x 0.0642per deg(l - 0.447) (4.50) 

= 0.0774perdeg < Cl u fixed — 0.0785per deg 


The neutral point is at 


h„ free = 2.78 


0.073 lper deg 
0.0774per deg 


(2.78 - 0.27) = 0.4094 < h nfixed 


= 0.443 


(4.51) 


Then C Ma free = 0.0744(/i - 0.4094). 


Horizontal Stabilator Configuration 

With this configuration, the elevator deflection is mechanically linked to the horizontal stabilator deflection 
as follows, 

8 e = k e i t + 8 a (4.52) 


The deflection 8 a is used to provide zero stick force at trim. The hinge moment at the horizontal tail is given 
by 


C ht = 


dC ht 

da t 


ot t + 


dC hl 
d8 e 1 


(4.53) 


Recall that the tail angle of attack a, in a wing-tail configuration is given by a, = i t + a w — e (as in equation 
(4.6)). Thus the floating incidence angle i t at the horizontal stabilator is obtained by letting C/,, = 0, 


C h , = 0 = 


9C/ "r uu 

—— (it + a w 

da t 


-o + 


dC hl 

d8 e 


(k e it + So) 


(4.54) 
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or 


dc h , 

it = B e { ——(1 - € a )a w + 
aa. 


dC h , 

8 S e 




where the constant B e is defined as 


Be 


-1 


3 Ch, , 3 Ck, L 
da, ^ d S e Ke 


(4.55) 

(4.56) 


With the above tail incidence angle expressed as a function of a w and 8 a , one can then determine the 
corresponding tail lift coefficient as follows, 


C Lt = a,a t + —^8 e (4.57) 

00 e 


After some simple algebra that proceeds roughly along the following line, 



dC Ll 

Cl, — + a xu e) + „„ (keit + <5 0 ) 

oo e 

(4.58) 


Cl, — (7/ (1 e a )(x w + G e 8 a 

(4.59) 

where 

1 dC L , dC h , 

Fe = 1 + (1 +- ke)B e -A 

a t dd e da, 

(4.60) 

and 

^ , , 9C +, , D dC ht , dC Lt 

Ge — (,a t + k e )B e + 

udg GOif OOe 

(4.61) 


Note that in this case, the free elevatorfactor F e can be greater than unity; hence resultingin an improvement 
on static margin for the stabilator configuration. 


Example 3:[Anderson] 

For a wing-body combination, the aerodynamic center lies 0.05c ahead of the center of gravity. The moment 
coefficient about the aerodynamic center is CM ac wb = —0.016. If the lift coefficient is Cl w1> = 0.45, what is 
the moment coefficient about the center of gravity? 

Note that C Mcg wb = C Mac wb + C Lwb (h - h ac<wb ) where h - h aCilvb = 0.05, C Lwb = 0.45 and C Mac wb = 
-0.016. Thus, C Mcg , wb = -0.016 + 0.45(0.05) = 0.0065. 

Example 4:[Anderson] 

A wing-body model is tested in a subsonic wind tunnel. The lift is found to be zero at a geometric angle 
of attack a — —1.5°. At a = 5°, the lift coefficient is measured as 0.52. Also at a — 1.0° and 7.88°, the 
moment coefficients about the center of gravity are measured as —0.01 and 0.05, respectively. The center 
of gravity is located at he = 0.35c. Determine the location of the aerodynamic center and the moment 
coefficient about the aerodynamic center Cm„, wb ■ 

Knowing the lift coefficients at different angles of attack ( CL wb = 0 at a = —1.5° and C= 0.52 at 
a — 5°) one can deduce the lift curve slope (as a linear approximation) a wb as follows, 

dC L . 0.52 - 0 

a wb = —^ g —— = 0.08per deg (4.62) 

da 5 — (—1.5) 
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Measuring the moment coefficients about the center of gravity at two different angles of attack and at the 
same time we know from previous calculation the lift curve slope, one obtains from 

Cm, cg , wb = CM acwb T &wb&wb(.h h ac ,wb ) (4.63) 

the following two linear equations in two unknowns Cm,,, wb and (h — h ac ^ wb ), 

-0.01 = C Macwb + 0.0m + C5)(h-h ac , wb ) 

0.05 = CM ac wb T 0.08(7.88 + 1-5)(/z h acwb ) ^ ’ 

Lrom equations (4.64), we can solve for Cm„ c wb and (/; — h ac , w b ) as 

C Mac , wb = -0.032, (h - h ac<wb ) = 0.11 (4.65) 

Since h = 0.35, then h ac ^ wb = 0.35 — 0.11 = 0.24. 


Example 5:[Anderson] 

Consider the wing-body model in Example 4 above. The area and chord of the wing are S = 0.1 m 2 and 
c = 0.1/77, respectively. Now assume that a horizontal tail is added to the model. The distance of the airplane 
center of gravity to the tail’s aerodynamic center is l t = 0.17/77, the tail area is S t = 0.02m 2 , the tail-setting 
angle is i, = —2.1°, the tail lift slope is a t =0.1 per degrees, and from experimental measurement e 0 = 0° 
and || = c a — 0.35. If a = 7.88°, what is the moment coefficient Cm, .. for this airplane model? Does this 
airplane have longitudinal static stability and balance? Lind the neutral point. 

Lrom equation (4.11), we have 

eft S t Ct 

Cm cs = C M ac „ ; + q Sc C Mact - Tj, V H a,(i t - e 0 ) + {(h - h ac<w )a w - rj t V H a t ( 1 - e a )}oi w (4.66) 

We further assume that the tail has a symmetric airfoil shape where Cm,,,. , = 0. Lrom previous example, we 
have CM ac , wb = —0.032, a w = 0.08, a w = 7.88° + 1.5° = 9.38" and (h — h ac .w ) = 0.11. Lurthermore 


V H 


1 ( assumed ) 

S,l t _ 0.02(0.17) _ A 04 
SC — 0.1(0.1) — 


(4.67) 


Thus 


C Mcg = -0.032- l(0.34)(0.1)(-2.7-0) + {0.11(0.08) - 1 (0.34)(0.1)(1 — 0.35)}9.38 = -0.065 (4.68) 
Lor longitudinal static stability, we examine Cm u as given in equation (4.12), 

Cm a = (h - h ac , w )a w - rj t V H a t ( 1 - e a ) (4.69) 


(4.70) 


C Ma = 0.11(0.08) - 1 (0.34)(0.1)(1 - 0.35) 

= —0.0133 < 0(Statically stable) 

Is the model longitudinally balanced? To find out we need to determine Cm„ (defined in equation (4.12) and 
from which we derive the equilibrium angle of attack. 

qt^tCt 


Cm 0 — Cm„ c w + 


qSc 


-CMact - r h V HCl( (i t - € 0 ) 


(4.71) 
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or 

C Mo = -0.032 - 1(0.34) (0.1)(—2.7 - 0) = 0.0598 (4.72) 

Thus, the equilibrium angle of attack is obtained by letting Cm cs = 0 in equation (4.11), or 


Cm,.,, — C M„ T C M a &equilibrium 

=> ^-equilibrium = -C M JC Ma = -(0.0598)/(-0.0133) = 4.4962° 


(4.73) 


This angle of attack is within reasonable limits; hence the airplane can be balanced and at the same time it 
is also statically stable. 

The neutral point is given by equation (4.19) as 


hn 


hac,wb T I a " b C\H ac f 


(4.74) 


where 

C La = ^ 1/7) 4“ 1 € a ) 

= 0.08 + 1 (0.02)/(0.1)(0.1)(1 - 0.35) 
haC't = h + l t /c = 0.35 + 0.17/0.1 = 2.05 


Thus, 


hn 


0.24 +[^ - 1J2.05 

0W3 

0.08 


0.493 


One can verities the above result using equation (4.20), namely 


0.093 


(4.75) 


(4.76) 


h fi a 

0.35 - 0.493 
-0.143 


C M JC La 

(—0.0133)7(0.093) 

-0.143 


(4.77) 


In the following we discuss some other effects that enter into our analysis of the longitudinal static stability. 

4.4 Other Influences on the Longitudinal Stability 

4.4.1 Influence of Wing Flaps 

Changes in the wing flaps affect both trim and stability. Themain aerodynamic effects due to flap deflections 
are: 

• Lowering the flaps has the same effect on Cm„ wb as an increase in wing camber. That is producing a 
negative increment in A Cm„ mb ■ 

• Theangleof wing-bodyzero-liftischanged tobemorenegative. Sincethe tailincidence i t ismeasured 
relative to the wing-body zero lift line, this in effect places a positive increment in the tail incidence 
angle i t . 

• Change in the spanwise lift distribution at the wing leads to an increase in downwash at the tail, i.e. 
e„ and jjC may increase. 
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4.4.2 Influence of the Propulsive System 

The incremental pitching moment about the airplane center of gravity due to the propulsion system (Ligure 
4.7) is 

A M cg = T Zp + N p l p (4.78) 


where T is the thrust and N p is the propeller or inlet normal force due to turning of the air. Another influence 
comes from the increase in flow velocity induced by the propeller or the jet slipstream upon the tail, wing 
and aft fuselage. 

In terms of moment coefficient, 


A C M , 


T Zp N p 1 p 

qS c qS c 


(4.79) 


Since the thrust is directed along the propeller axis and rotates with the airplane, its contribution to the 
moment about the center of gravity is independent of a w . Then we have 


AC Mo = 


T_ Zp_ 
qS c 


(4.80) 


and 

A C Ma = N prop S ^^ d ^( 1 - e a ) (4.81) 

Sc da 

where the propeller normal force coefficient dC^ p /da and the downwash (or upwash) e a are usually de¬ 
termined empirically (Ligure 4.8). N prop is the number of propellers and S prop is the propeller disk area 
(= jtD 2 /4) and D is the diamctcrof the propeller. Notethat a propeller mountedaft ofthec.g. isstabilizing. 
This is one of the advantages of the pusher-propeller configuration. Note that n in Ligure 4.8 is the propeller 
angular speed in rps. 
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pY 2 D 2 


Figure 4.8: Propeller Normal Force Coefficient Cn pu = 9C g“ fl<fe f(T) 
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4.4.3 Influence of Fuselage and Nacelles 

The pitching moment contributions of the fuselage and nacelles can be approximated as follows (Perkins & 
Hage p. 229, Equation (5.31)), 

K f WlLf 

CM afuse , age =-—-(per degrees) (4.82) 

where W/ is the maximum width of the fuselage or nacelle and Lf is the length. The empirical pitching 
moment factor Kf is given in Ligure 4.9 (NACA TR 711). 



Ligure 4.9: Kf as a Lunction of the Position of the Wing c/4 Root Chord 


4.4.4 Effect of Airplane Flexibility 

Flexibility of an airframe under aerodynamic loads is evident in any flight vehicle. The phenomenon that 
couples aerodynamics with structural deformations is studied under the subject of aeroelasticity. There are 
two types of analysis: 

• Static behaviour: Herethesteady-statedeformationsof thevehicle structurcareinvestigated. Phenom¬ 
ena such as aileron reversal, wing divergence and reduction in static longitudinal stability fall under 
this category. 

• Dynamic behaviour: The major problem of interest is associated with the phenomena of dynamic 
loading, buffeting and flutter. 

Let’s study, for example, the effect of fuselage bending on the tail effectiveness. It can be seen that the angle 
of attack at the tail is reduced by the fuselage bending according to the following equation. 


a, = oiwb + i t - e - kL t 


( 4 . 83 ) 
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The tail lift coefficient (with S e = 0) is 


C L , = a,<x t = a,(a wb + i t - e - kL t ) 


(4.84) 


or 


Cl, - a,(a wb + i t - e - kq t S t C L ,) 


(4.85) 


Solving for Cl,, we get 

Cl, = Gt — (a wb + i,-e) (4.86) 

1 + kr),qS,a t 

Thus the tail effectiveness is reduced by a factor 1 /[I + kq t qS t a t \ that decreases with increasing speed V in 
the dynamic pressure q. This decrease in the tail lift curve slope will cause the neutral point to move forward 
(i.e. .reduced static stability). 

Similarly, it can be shown that the elevator effectiveness is decreased due to fuselage bending since 


C L , = a t (a wb + i t -e - kL,) + 



(4.87) 


or solving for C/ ( , we obtain 

Of {plwb + it — €) + 

C, = -^— (4.88) 

1 + kq t qS t a, 

Thus the elevator effectiveness is reduced by the same factor 1 /[I + kq t qS t a t ]. 


4.4.5 Influence of Ground Effect 

When the airplane is near the ground to within 20% of the wing span, the wing and tail lift curve slope will 
increase by about 10%. At the same time, the downwash is reduced to about half of the normal value, which 
requires a greater elevator deflection to hold the nose up. However, static stability is usually improved by 
the ground effect. 

The aircraft must have sufficient elevator effectiveness to trim in ground effect with full flaps and full 
forward c.g. location, at both power off and full power. 
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Chapter 5 

Static Longitudinal Control 


We have studied in Section 4 the concept of longitudinal stability of an airplane in trim. It was shown that 
static stability is primarily governedby the sign of the derivative of the moment coefficient aboutthe airplane 
center of gravity with respect to the angle of attack, i.e Cm„, being negative. All the above analysis relies on 
the fact that one can trim the airplane. The question is what arc the controls that allow us to trim the airplane. 


5.1 Longitudinal Trim Conditions with Elevator Control 

For a steady level flight it is easily seen that the airplane velocity in trim is given by 


F trim — 


2W 


psc Ln 


(5.1) 


Thus if the pilot wants to fly at a lower velocity V < V tr im, then from equation (5.1), we must have C 
(or the angle of attack) increased in order to offset the decrease in dynamic pressure. But increasing the 
angle of attack away from trim would generate for a statically stable airplane a negative pitching moment 
that tends to bring the angle of attack back to the original trim point (Figure 4.4). It would therefore be 
impossible to change speed if nothing else is changed about the airplane. It turns out that there are basically 
two ways to achieve a change in the trim angle of attack. The control concepts arc illustrated in Figure 5.1. 
One possibility is to change the slope of the moment coefficient curve as indicated in Figure 5. la. Thus, by 
decreasing the slope Cm u (he more negative), one can achieve a smaller trim angle of attack and hence one 
is able to fly at a faster velocity. If we examine equation (4.12), the only way to modify Cm„ is to change the 
location of the airplane center of gravity that shows up in both the variables h and Vh- This principle is used 
extensively in modern hand gliding craft but it is clearly not practical for large fixed wing airplanes. The 
alternative is tochange the value of Cm„ as indicatedin Figure5.1b. It will be shown below that by deflecting 
the elevatorin the horizontal tail one can translate the moment coefficient curve upward and downward while 
without affecting its slope. 

Let’sexaminetheeffectofdeflectingtheelevatoronthetailliftcoefficientcurve. Usingasignconvention 
of positive elevator deflection being downward (or using the right-hand rule for angle), it is clear that a 
deflected elevator causes the lift curve to shift upward and to the left as shown in Figure 5.2 . The lift curve 
slope remains unchanged. Now, if we assume that the tail is at a constant angle of attack a t , say a t i in Figure 
5.2, then an increase in elevator deflection would lead to an increase in tail lift along the vertical dashed line 
(Figure 5.2). If we plot the tail lift coefficient Cl, as a function of 8 e when the tail is at a given tail angle of 
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Ligure 5.1: How to Change Airplane Trim Angle of Attack 



Ligure 5.2: Tail Lift Coefficient vs Tail Angle of Attack 


* ▼ 
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Figure 5.3: Tail Lift Coefficient vs Elevator Deflection 


attack a, (held constant), we would have a curve much like the one given in Figure 5.3 (where we assume 
that the slope stays nearly constant and does not change with 8 e ). This slope of the tail lift coefficient with 
respect toelevatordeflectiondenoted by 9 Cl, /d8 e iscalled the elevatorcontrolejfectiveness . This quantifies 
the effectiveness of the elevator as a control surface. It can be seen that this constant 9Cz, r /35 e is always 
positive. With the above definition, one can from here on express the tail lift coefficient as a function of two 
independent variables a, and 8 e . In the form of a first-order Taylor series expansion, we have 


Cl, = 


dC L , 

3o 8 e 


a, + 


9 C Ll 
93 e a, 


(5.2) 


or since a ( = then 

Cl, = a t a t + ——(5.3) 
00 e 

Substitutingequation(5.3)intoequation(4.1 l)forthepitchingmomentcoefficientaboutthecenterofgravity, 
we have 


S,c, dC Ll 

Cm c „ — C Macw T 3t~ Cm„ c , + Li w a w (h h ac w ) lyVLiicitWt T ~rr 3 e ) 
Sc 38 e 


(5.4) 


SC Me, 


Then the rate of change of Cm c due to elevator only is defined as ... CH . From equation (5.4), we obtain 


dC Mcg 

dS e 


= -r) t V H 


9 C Ll 
dS e 


(5.5) 


Notethatsince isalwayspositive, wededuce that C ^ cg isalways negative. Thus, anincrementalchange 
in Cm cs for a given elevator deflection 8 e is simply, 

A C Mcg = -n,V H ^8 e (5.6) 

So by deflecting the elevator one can shift the moment coefficient curve downward by an amount ACm,, : 
given in equation (5.6). This confirms the behaviour depicted in Figure 5.1 b, where elevator control can 
be used to change the trim point. Moreover, from equation (5.4) we can show that the slope of the moment 
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coefficient curve with respect to angle of attack is not affected (to first-order approximation) by the elevator 
deflection. Only the value of Cm b is modified by elevator deflection. Namely, 


C M cg 


(C Mo + a C Mcg ) + -§£-a w 
(C Mo - VtVn^Se) + 3 ~^a w 


(5.7) 


5.1.1 Determination of Elevator Angle for a New Trim Angle of Attack 


The problem is to find the elevator deflection 8 etr i, n such that the moment coefficient equation is balanced at 
a new angle of attack a n . Wereturn to equation (5.7) where we substitute a w by a n and using equation (5.6) 
for A C Mcg 


Cm cs = 0 = Cm„ — V t Vn 


dC Ll 

d8 e 


+ 


dC Mcg 

da w 


a n 


(5.8) 


Solving for <5 etrim , we obtain 


Cm,, + 


ac Mcg 

da w 


Un 


11<Vh 


3 C L , 
dS e 


(5.9) 


Example 6 [Anderson] 


Considerafull-sizeairplanewiththeaerodynamic chai'acteristicsdefinedfortheaiiplanemodelin Examples 
4 and 5 of Section 4.3. The full-size airplane wing area is S = 19 m 2 with a weight of W = 2.27 x 10 4 A, 
and an elevator control effectiveness of 0.04. Determine the elevator deflection angle needed to trim the 
airplane at a velocity of V = 6\m/s at sea level. 

Lirst, we need to find the airplane angle of attack to fly at V — 6\m/s. It is given by 


2W 

Cl ~ pLAS 


2(2.27 x 10 4 ) 
1.225(61 ) 2 ( 19) 


(5.10) 


Lrom Example 5, we have C/.„ = 0.093. Then the absolute angle of attack of the airplane is 




C L 

C La 


0.52 

0.093 


5.59° 


(5.11) 


Lrom equation (5.9), we have 


or 


8 


etrim 


Cm,, + 


3 CMcg 
da w 


Oi n 


il,V h 


a Cl, 
dS e 


0.0598 + (-0.0133X5.59) 
1(0.34)0.04 


-1.0696° 


(5.12) 

(5.13) 


5.1.2 Longitudinal Control Position as a Function of Lift Coefficient 

Lrom equation (4.16) we have expressed the total lift coefficient as a function of angle of attack at the wing, 
constant component of downwash and the tail incidence angle, 

C L = C La a w + rj t -^a t {i t - e„) 


(5.14) 
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or 

C L = C La u w + C Ltii (i, — e 0 ) (5.15) 

with 

C Llli = ru -ja t (5.16) 

Note that Cl„ is given in equation (4.17). From equation (5.15), we see that Cl is a linear function of angle 
of attack a w . Thus for a given Cl and tail incidence angle i t , one can solve for a,„ as 


oc w 


Cl C L ,j t {it £o) 

c La 


(5.17) 


The other equation of importance is the one for the pitching moment about the airplane center of gravity as 
given in equations (4.11)-(4.12), 

Cm cs = Cm„ + C Mu a w (5.18) 

where Cm„ and Cm,, arc as defined previously in equations (4.12) where Cm„ can also be expressed as 
Cm,, = Cl„ (h — /;,,)■ Substituting equation (5.17) and equation (4.12) into equation (5.18) the pitching 
moment equation is now a function of the total lift coefficient Cl and the tail incidence angle i t . Namely, 


Cm cs = Cm cu 


S,c, 

+ rit-^ 1 C Mai 
Sc 


ri,V H at(it - e a ) + C Ma 


C L ~ Cl,,, (i t ~ C>) 
C La 


(5.19) 


For an aiiplane in trim, one must have Cm c = 0 in equation (5.19). Then one can solve for the tail incidence 
angle at a particular Cl trim, 

it = A r + B t C L (5.20) 


where 


At — Co + 


CMac w + ‘TT CM at 


C M^L, h 

c L a 


+ViVho, 


and 


■ CM ac , w ~C r lt ~Jc~CM ac t 

e ° r),A at [ l ±+h-h n \ 


B, 


_ Cm„ _ 

CM a CL til +VtVLa 


h-h„ 

rh^a,['-C+h-h„] 


(5.21) 


(5.22) 


Let’s study the sign of the coefficient B,. Note that for a statically stable aiiplane, Cm,, < 0. It can be easily 
shown that the denominator term CM a CL, h + Vt Vr/citCr,,, is equal to 


S 1 

Cm„Cl, h + r\ t V H atC La = Vt^a t C La (h - h„ + C) 


(5.23) 


Thus if l t /c > h n — h, then the coefficient B, will be negative. Equation (5.22) allows one to determine 
experimentally the neutral point. This is done by measuring i, as a function of C/, for differentc.g. locations. 
The slopes of the experimentally derived curves are then plotted as a function of center of gravity locations 
(i.e. h). The neutral point is determined by extrapolation to find the value of c.g. that gives a zero slope 
in dif/dCL (see Figure 5.4). According to equation (5.20) this simply corresponds to having Cm,, = 0 or 
h = h„ in equation (4.20). 
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Ligure 5.4: Determination of Stick-Lixed Neutral Point from Plight Test 

5.2 Control Stick Forces 

Pilotsusecontrolstickforcesasoneofthemeansofevaluatingtheflyingqualitiesofanairplane. Thus,ability 
to control an airplane is quantified in terms of required maximum exerted control forces and its sensitivity 
with respect to airspeed. Lor simple mechanical control systems, the pilot controls are directly linked to the 
respective control surfaces and the forces he must exert arc proportional to the hinge moment (generated 
primarily from aerodynamics) about the pivot point at the surfaces. 

Let’s review the key equations governing the analysis of control hinge moments. 

In the design of airplane control system, the stick (or control wheel) forces must he within acceptable 
limits throughout the operating envelope (V-n diagram) of the airplane. And the gradient of these forces 
with respect to airspeed at trim point must produce the proper “feel” to the pilot. In general, the pilot tends 
to push forward in the longitudinal control to fly faster and pull on it to slow down. A requirement stated in 
EAR Part 23 poses a limit on the maximum force of 60 lbs for the stick and 75 lbs for the control wheel. 

Lor analysis, we often assume that the control force P is proportional to the hinge moment H at the 
elevator. Namely, 

P = GH (5.24) 

It can be shown that for a system in equilibrium we have Ps + H8 = 0. Or 



H 


Ligure 5.5: Longitudinal Control Stick to Stabilator 


P = 





(5.25) 



5.2. CONTROL STICK FORCES 


51 


Thus, G = —S/s is the gearing ratio which, as derived here, is totally independent of the details of the 
mechanical linkage. Since 8 is negative for a positive stick displacement s as shown in Figure 5.5, the hinge 
moment would be positive for a positive stick force; thus G is positive. 


5.2.1 Stick Force for a Stabilator 

For a symmetrical airfoil, the pitching moment at the horizontal stabilator is given by 

Cm, = Cm,„ c (= 0) + CM tolt oc t + Cm, $ e 8 e (5.26) 

From equation (5.24), we deduce the stick force for the stabilator to be 

P = Gq,S t c t {C Mlu a t + C M , Se 8 e } (5.27) 

Note that C m, can be positive or negative depending on whether the aerodynamic center of the stabilator is 
ahead or behind the pivot. The coefficient CM, Se is generally negative. The first step is to express a t and 8 e 
in terms of Cl- First we recall that the elevator deflection is linked to the horizontal tail incidence as 


8 e — k e i t T 8 0 


(5.28) 


The trim angle of attack a w is determined from the total lift coefficient and the tail angle of incidence i, from 
the pitching moment equation in balance. Namely, 


S, d C, 

C L = a w a w + r),—{a t [i, - e Q + (1 - e a )a w ] + ——(k e i t + <5 C )} 

ib o Op 


(5.29) 


S t dC L , 

Cl = C La a w + C L F e i, + tj t — [—— S 0 - a t e Q ] 

5 od e 


where Cl„ is given in equation (4.17), 


and 


1 3 C L , 

Fe = 1 + ~^T-k e 


r S ‘ 

Cl,,, = rit—ci, 

One can solve for a w in terms of C/, . e 0 as follows, 

1 S t d C L , S, 

a w = t ~{Cl ~ Vt —^77^8,, + iy—a t e 0 - C Lti F e i t } 

C J_, a ib uOf ib 

The other equation we use is the pitching moment about the airplane center of gravity, 

Cm cs — Cm,, + C Ma a w 

where 

Cm 0 — C m ucw + Tit——C M acJ — VH\a t (F e i, — e 0 ) H— 
ibc 


(5.30) 

(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 


and 


Cm u — (h h ac , w )aw Cnu, (I c a ) — Cl„ (h h n ) 


(5.36) 
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Using equations (5.28), (5.33) and the fact that in trim Cm cs = 0, we can solve for i, in terms of Cl, A, 
We obtain 


where 


and 


A, = 


r)t^a t F e ( l j; + h- h n ) 


it — A s + B s Cl 


S t c t 1 1 dC L , 

{C Mac , w + Vt—jj—CM ac ,t} + TT* 0 - —^<5 0 

Sc r g Clf i*£ OOg 


R, = 


li - h n 


r] t ^a,F e ('j + h - h „) 


Substituting i t of equation (5.37) into equation (5.33), we obtain 

®u> = A a + R q Cl 


where 


and 


1 S,c t 

Aa — —~f - - -——— {Cm„ C iW + Cm„ C ',} 

(i + h- h n )C La Sc 


Bn = 


Vi + h- h n )C La 


The stick force P given in equation (5.27) becomes 
P 


Grj t qS t c, 


— \{Cm, a , + C MlS k e )it CM, a ,( 1 — ^ a )a w — C m, to + Cm, s< ,^o\ 


(5.37) 

(5.38) 

(5.39) 

(5.40) 

(5.41) 

(5.42) 

(5.43) 


— - - — — [(Cm,„, + CM tSe ke)(A s + B s Cl) + (1 — €u)(A a + BoCl) ~ CM ta ,€o + (5-44) 

UiitqofCt 

Thus the stick force P is a linear function of Cl- By collecting all the terms we have 

P 


where 

and 


= A + BC l 

GqtqStC 

A — [(Cm, a , + CM, s k e )A s + CM ra ,(l — € a )A a — CM ta ,^o + CM tSe &o\ 


(5.45) 


(5.46) 


B = [{C Mtat + C MtS k e )B s + C Mtat (1 - e a )B a ] (5.47) 

Notice that the parameter 8 a in A can be used to achieve P = 0 at a particular - trim velocity V, rwt . That is, 

0 = A + BC Ltrim (5.48) 

or A — —BCL lHm ■ We can substitute Cl by 

W W 
Cl = — = 


qS kpV 2 S 


(5-49) 
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Using equations (5.45) and (5.49), it can be easily shown that 

= <,B(C L -C L „J 

= ?«c I ( i-%-) 

- l^lRd ^ t 

(nm 


(5.50) 


Note that V tr j m is given in equation (5.1) and W/S is simply the wing loading. 

From theabove equation, wecan obtain thegradient ofthestick force with respect to speedat V = V tr im , 


dP 

dV 


-2Grj t S t c t 


W B 


S V u 


(5.51) 


For a given trim speed, a proper stick force gradient as stated in FAR Paid 23 must be negative for all 
conditions of flight. Or B must be positive. Notice that the gradient is large if V tr j m is small. Figure 5.6 
shows a typical stick force versus speed curve described in equation (5.51). At V = 0, clearly from equation 
(5.50) we have 

P = GritStC^B (5.52) 



5.2.2 Stick Force for a Stabilizer-Elevator Configuration 

We proceed as before by first giving the hinge moment associated with the stick force P, 


Ch = Ch 0 + 


dC H 

dot, 


at + 


dC H 

dS e 


S e + 


dC H 

dS r 


S t 


(5.53) 
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where in general Ch„ = 0. Next we examine the total lift coefficient 

S t dC Ll S, 


Cl = C La d w + rj t 


c Vt „ at ^° 

u doo »b 


where we assume that the lift contributed by the trim tab is negligeable. Solve a w in terms of S e 

1 S, 3 C L , S t 

u w - —{Cl - rit — ——&e + Vt— a t e 0 ] 

C ]_, kb u Og »b 


Now we consider the pitching moment. 


Cm„ - Cm,, + Cm„o/ u 


where 


and 


Cm„ = CM acw + Vt~^—CM ac , + >1tVH[a t e 0 - 

Sc ' ob e 


C M a — (h h acw )o w 1 ^ (/ ) 

Lrom trim balance with C,v/, „ = 0, we solve for 8 e 

8 e = A e + B^Cl 

where 


A, = 


1 S t c, a t 

S, 3 Cl, A, , , —{CMcc* + 

Vt S d8 e (c + 2 1,1 ) 


dC L , 

dS e 


and 


B„ = 


h — h n 


S d s' (c ^ hn) 


Similarly, by substituting equation (5.59) in equation (5.55), we can express a w in terms of Cl, 

— 4 [i + B ci Cl 

where 


A n = 


1 S,c t 

■ 77 —-— —{CM ac , w + m-^CM^,} 

(7 + h — h n )CL a Sc 


and 


k 


B rl = 


(h + h — h n )CL a 
Lrom equations (5.24)) and (5.53), the stick force P becomes 

P 3C« s 3C» dC H „ 

— - 7 — — —(1 — € a )a w + 8 e + 3 ( ] 

do e dd t 

Thus the stick force P is a linear function of Cl- By collecting all the terms we have 

P 


C TjfCjS e c e 


(5.54) 
and Cl, 

(5.55) 

(5.56) 

(5.57) 

(5.58) 

(5.59) 

(5.60) 

(5.61) 

(5.62) 

(5.63) 

(5.64) 

(5.65) 


— A e + B e CL 


(5.66) 
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where 


and 


3 C H dC H 3 C H 

Ae = hr-^(l - + ~f~A e + 

uQtt 0&e vut 

3 C H 3 C H 

Be = [—-{l-e a )B a + —-B e ] 

OOtt od e 


(5.67) 


(5.68) 


Again the trim tab 8 t is used to achieve P — 0 at V = V tr im , then one can simplify the above equation for 
the stick force P to the following. 


P IV - V 2 

= (—)B e { 1--2-} 

v trim 


GTjtSeCg S 

From the above equation, we obtain the gradient of the stick force with respect to speed at V — V tnm , 


(5.69) 


clP 

dV 


v=v, rim 


-2Gr] t S e c e 


IV Be 

5 V trim 


(5.70) 


5.3 Steady Maneuver 


We now consider the determination of the elevator angle per g in a pull-up maneuver. In the analysis, the 
concepts of stick-fixed and stick-free maneuver margins are introduced. For an airplane in a steady pull-up 
maneuver, the lift force will exceed the vehicle weight. Namely, 

L = W( 1 + —) (5.71) 

8 

where a n is the vehicle acceleration. The lift to weight ratio is known as the load factor n or 

n = 1 + — (5.72) 

8 

For a straight level flight, a n = 0 and n — 1. 

We examine the pitching moment of the airplane in this pull-up maneuver from which we derive the 
quantity known as elevator angle per g. Again we have 


Cm c = Cm„ + Cm & oi w + A Cm (due to airplane rotation in a steady pull-up maneuver) 


where 


qr),S t a,Aa,lt 

A C M = ---= -rj t V H a t Aa t 

qSc 


The incremental angle of attack at the tail due to a constant angular velocity Q is 

l,Q 


Aa, = 


V 


A 21 1 _ 
Aa, = —q 

c 

with q being the dimensionless pitch rate variable defined as 

Qc 


<7 = 


2V 


(5.73) 

(5.74) 

(5.75) 

(5.76) 

(5.77) 
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5.3.1 Horizontal Stabilizer-Elevator Configuration: Elevator per g 

When the airplane is in straight and level flight (unaccelerated), the elevator angle and stick force to trim arc 
8 e and P respectively. When in the pull-up maneuver, the elevator is deflected to 8 e + A8 e and the stick force 
required is P + A P. The quantities A8 e /{a n /g) and A P/(a n /g) arc known as the elevator angle per g and 
the stick force per g respectively. These arc measures of airplane maneuverability; the smaller they arc the 
more maneuverable it is. 

Including the effects due to pitch rotation Q in a pull-up maneuver, the pitching moment becomes 

dCjy[ dC-M 

C Mcg = C Mo + C Ma a w + —^8 e + —^Q (5.78) 

od e o {2 

Note that Q is not dimensionless and it has the units of rad /s or deg /s. To nondimensionalize Q we use 
the variable q defined in equation (5.77). Then equation (5.78) becomes, 

dC M 

Cm cs = Cm 0 + Cm„oi w + —— 8 e + CM- q q (5.79) 

00 e 

We can derive Cm ■= from A Cm and the expression of A a t due to Q, 

C Mq = -V,a,V H — (5.80) 

c 

This term Cm- is often referred to as the pitch damping term, since it produces a negative pitching moment 
due to a change in pitch rate. 

In a steady pull-up there is still no angular acceleration in the pitch axis, thus for equilibrium Cm cs = 0 
as it is in trimmed straight level flight condition. The increment A Cm resulting from a steady maneuver is 


A Cm cs = 0 = C Ma Aa w + CM Se A8 e + Cm 


(5.81) 


where we assume that Cm 0 remains constant in the maneuver. Lrom the above equation, one can solve for 
A 8 e as 


AS e 


Cm, &oi w + C M ,q 
Cm Sc 


(5.82) 


It can be shown that since a n — QV we have 


Qc a„c 
2V ~ 2U2 


gc a n 

zv*g 


(5.83) 


Again the incremental change in angle of attack Aa w is determined from 


ACL = C La Aa u , + ^A5 e 

O Op 


(5.84) 


where A Cl is the incremental lift coefficient due to the pull-up; namely 


Cl + A Cl L a„ 

C L ~ W~ + g 

(5.85) 

AC L = —C L (for 1 g flight) 

(5.86) 
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Then 


Note that 


or 


A a w 


Cl fln C Ue 

cZ~g-~cT 


_9Cl_ S t d C Lt 
Ue d8 e Vt S dS e 


c L Se = V t 


s, dC Ll da, 
S da T dS e 


S t d C L , 

rjt - 

S da, 


T 


S, 

rj t —a t T 


where the coefficient r can be determined from Figure 5.33 (Perkins & Hage p. 250). 
Substituting equations (5.83) and (5.87) into equation (5.82) we have 


(5.87) 

(5.88) 

(5.89) 


AS e _ C Mg C L +C M ,C La (AA) 

(<*n/g) CMs e Cl u ~ CM a CL Se 


(5.90) 


Note that W — qSC l = | pV 2 SCl and W = mg where m is the mass of the airplane. Hence, Cl = • 

Also recall that Cm s = —)/, V / //C/., from equation (5.6). Then 


A8 e = _ C C k ~ h " + 477' C Mg 

0 On/g ) L La C M& C La -C Ma Cu 

or letting /x = (known as the relative mass parameter), then 

A8 e _ _ c c h ~ 

(a n /g)~ L La C Mh C La — C m „Cl, 


(5.91) 


(5.92) 


Similarly to the neutral point, there is also a particular - value of h, known as the stick-fixed maneuver point 
denoted here by h m , forwhich no (i.c. very small) elevator willbe required toproduce a finite (i.e. not small) 
acceleration. This is determined from equation (5.92) by setting A S e = 0; namely 


h m — h n C M , 

2/x 


(5.93) 


Since Cm- <0, then h m > h„ or the stick-fixed maneuver point lies aft of the neutral point. The quantity 
h m — h is known as the stick-fixed maneuver margin. And we can rewrite equation (5.92) as 


or 


A<$e _ C L C La (h m - h) 

( a n /g ) C M, e Cl u — CM a CLs e 

A8 e _ _ C L (h m - h) _ 

C L% [-il,V H - (h - h„)ri t %] 

A 8 e C L (h m — h) 

- — -f;- < 0 

(a n /g) C LlSe il,f(h n -h acJ ) 


(5.94) 

(5.95) 

(5.96) 
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5.3.2 Horizontal Stabilator Configuration: Elevator per g 

A similar derivation can be performed for the stabilator configuration and we obtain the following results, 

A i, h ~ h n + 

-— - ~C L C, - — -— (5.97 

(a n /g) a C Ml C La -C Ma C Lit 

where we replace Cm Je by Cm,, and C/, ; by C /. il . 

5.3.3 Stabilizer-Elevator Configuration: Stick Force per g 


Recall from equation (5.53) that 


9 C H dC H 0 dC H „ 

C ii — — — a, + „ „ S e + „ „ S t 
oa t U0 e vOf 


assuming that the trim tab has negligeable contribution to the total lift and moment. In a steady pull-up, 


or, assuming e a — 0, 


Then the stick force P is given by 


a, = i t + a w - e + 2-q 
c 


a, = i t + (1 - e a )a w + 2-q 
c 


(5.100) 


p dc H . 2i t _ dc H „ ac ff „ 

-p; -^ — I -J Ut + (1 — ^a) a w H- q] + + “TST - &A 

G q,q8 e c e da, c do e ad, 


(5.101) 


Let 8, be adjusted to achieve P = 0 at straight level flight (unaccelerated). Then a w = a Q + Aa u 
8 e = 8 eo + A 8 e . It then follows that 


A P dC H 21, dC H 

--7T- = hr-^Kl - e a )Aa w + —q] + —^A«y 

Gq,q8 e c e da, c dd e 


(5.102) 


Recall that 


21, _ gif a„ 

—q = T72<—) 
c v- g 


(5.103) 


Thevariables A a w and A 8 e are thoseobtainedin thepreviousanalysis forelevatoiper g asgiven inequations 
(5.87) and (5.96) respectively. Thus, equation (5.102) becomes 


A P 8Ch Cl a n Cl* git a n 9Cy 

- -— = - e ff ){-^(-) - —^AA e } + + ~^A8 e ] 

Gq,qS e c e da, C La g C La V 2 g d 8 e 


(5.104) 


After some simple manipulations, we obtain 


A P W dC H (1 - e„) 1 l, d -w~C La - (1 - e a )C Ls 

—— = Gq,S e c e (—){-—^-\—— — + -d] + ss ‘ 3 °' -i( h m - h)} (5.105) 

(a n /g) 5 da, Cl„ A c Cm^Cl„ — Cm„Cl^ 
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Again, there is a position of center of gravity h for which AP = 0. This is known as the stick-free maneuver 
point denoted by h' m . It is called the stick-free maneuver point since it corresponds to the fact that the pilot 
can let go of his control stick (i.e. stick free). Equation (5.105) can be rewritten as 


AP 

C a„/g) 


dCjLc r 

. ax. 


r- c 3S, 

Gr) t S e c e (—){ --—— 

^ ^ Ms, ^ L. 


da, u 


€a)Cc s 


Cm u Cl s 


(h' m - h)} 


(5.106) 


where 


h m — h m + 


dC H (l-e a ) 
da, C La 


1 l t CM Se CL a — 

+ ^~ c a -§?C La - a -^(l-e a )C Lsi 


(5.107) 


One would create a catastrophic situation if we load the airplane such that h = h' m . In this case, the pilot 
would inadvertently generate extremely large inertia loads on the airplane by exerting little or no control 
force. Such a situation rarely occurs if the vehicle has an adequate static margin since the stick-fixed and 
stick-free maneuver points are always aft of the neutral point. 


5.3.4 Stabilator Configuration: Stick Force per g 

A similar analysis as the one performed for the stabilizer-elevator configuration can also be done for the 
stabilator configuration. Details arc left to the students. (They can be found in the reference book by B. W. 
McCormick) 

Remarks: 

• Stick force per g is a linear function of /?. For a normal range of c.g. (i.e. h < h ' m ) locations, the 
force gradient is positive. That is the pilot must exert a pull (positive) force on the stick to maneuver a 
pitch up maneuver. The airplane is highly maneuverable if this force gradient is small (i.e. when the 
c.g. location is near the stick-free maneuver point). In this case, the airplane may even be statically 
unstable. 

• Stick force per g is directly proportional to wing loading (W/S). 

• It is independent of Cl or V apart from Mach and Reynold effects. 
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Chapter 6 

Lateral Static Stability and Control 


The concept of static stability and control we studied for the longitudinal axis can also be applied to the 
lateral axis. The key motion variables in the lateral axis correspond to sideslip (with sideslip angle /J, or side 
velocity v), roll (with roll rate p) and yaw (with yaw rate r). Primary controls arc rudder 8 r and ailerons 8 a . 

6.1 Yawing and Rolling Moment Equations 

The lateral motion of an airplane is described in terms of two tightly coupled motions: yaw 
axis (i.e. directional) and roll about the a- body axis (i.e. lateral). 

Using Figure 6.1, let’s first identify the motion variables, the control effectors and their 
used in the lateral-directional stability analysis. 

Note that all along we adhere strictly to the right-hand rule for sign conventions: 

• x-body axis pointed forward, 

• y-body axis pointed to the right wing, 

• z-axis pointed down toward the earth, 

• Positive yaw motion (clockwise) (i.e. yaw rate r > 0 with yaw angle i//), 

• Positive roll motion (right wing down) (i.e. roll rate p > 0 with roll angle r/j), 

• Positive sideslip /3 — sin ~ x {v/V) corresponding to a positive u-componcnt in side velocity. 

• Aileron control positive with right aileron down ( 8 ar > 0 again consistent with the right-hand rule 
where a positive rotation about the positive y-body axis) and left aileron up (8 a i > 0). In general, we 
define aileron angle 8 a to be 8 a = 8 ar + 8 a /. Notice that this definition of positive aileron will produce 
a negative rolling moment. 

• Positive rudder deflection ( 8 r > 0) is to the left (again according to the right-hand rule a positive 
rotation about the positive z -body axis). With this definition, a positive rudder will generate a negative 
yawing moment and a positive side force. This is similar to the definition that a positive elevator 
deflection would pitch the airplane down, i.e. producing a negative pitching moment and a positive 
contribution to lift. 


about the z-body 
sign conventions 
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i x 



Ligure 6.1: Definition of the Lateral Directional Motion of an Airplane 
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6.1.1 Contributions to the Yawing Moment 

There arc numerous components that contribute to the yawing moment in the airplane when perturbed in the 
lateral motion variables. Note that in straight and level flight, we usually have ft — r — p = 0 in trimmed 
flight. 

The airplane yawing moment about the center of gravity can be written as 


N C g — Nwing T A fuselage T A verticaltail T A propulsion T A rudder T A aileron (6.1) 

or in terms of moment coefficients, 

_ N _ 

^N qSb ^^wittg ^ Nfuselage ^Nverticaltail Npropulsion ^ Nrudder ^Naileron ( 6 - 2 ) 

• Wing Contribution N w i ng : Yawing moment generated at the wing is developed mainly from perturbed 
motions in sideslip ft and roll rate p in the lateral axis. Let’s examine some of these effects: 


- Due to sideslip, there is an increase in drag on one side of the wing that is more perpendicular to 
the flow and thereby would produce a yawing moment. If the wing is swept aft, then this yawing 
moment is stabilizing (i.e. producing a positive yawing moment to a positive sideslip). An 
empirical formula for the wing yawing moment coefficient due to sideslip, C,v„, f , = dC^/dft, 
is given by 


t 1 tan A A 

C N . = Cj {-[-] [cos A- 

4i tA tcA(A + 4cosA) 2 


A 2 6(h aC:W - h)sin A 
8CO.SA A 1 


(6.3) 


where A is the wing aspect ratio, A is the wing sweep angle, Cl is the total lift coefficient, li ac ,w 
is the location of the wing aerodynamic center in percent chord and h is the location of c.g. in 
percent chord. 

- Due to roll rate, the right wing will move down and thereby see an increase in angle of attack of 
py / V applied to a wing section located at a distance y from the root. This increase in angle of 
attack will tilt the lift vector forward on the right wing, while on the left wing the inclination is 
to the real - . As a result, a differential yawing moment is created and is given by 


dN w = —q(cdy )C/(—)(2y) (6.4) 

or 

dN w = -2 ^q(cC/y 2 dy) (6.5) 

Thus integrating from 0 to b/2, we obtain the incremental yawing moment due to roll rate p as 


AN W = 



(6.6) 


or, in terms of the yawing moment coefficient, we have 


p f h A 2 

A C Nw = -2 -A— / cCiy 2 dv 
SbV Jo 


(6.7) 
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Lor a linearly tapered wing with taper ratio X = c t /c 0 and assuming that C/ is constant across 
the wing span, we can show that 


A C Nm 


Cl 

12 


1 + 3X 
1 T X 


( 6 . 8 ) 


or the yawing moment coefficient contributed by the wing due to roll rate is given by 


C Nu , p 


C L n +3X\ 
12 \ 1 + A. / 


(6.9) 


in terms of the dimensionless roll rate p = pb/2V. 


• Luselage (and Nacelle) Contribution Nf use iage '■ Yawing moment due to sideslip is a function of the 
fuselage (or nacelle) volume, length and width as follows, 


9Cat, 


C N 


fuselage 


, „ Volume Df 

-1.3 — -( —f -) (per radians) 

Sb Wf 


( 6 . 10 ) 


where Wf and Df arc respectively themaximum widthanddepth ofthe fuselage. Clearly, the fuselage 
produces a negative contribution to the lateral stability, i.e., making the vehicle less stable in the yaw 
axis. 


• Vertical Tail Contribution N ver ti ca it a ii : The vertical tail plays as significant a role in the lateral motion 
as the horizontal tail in the longitudinal motion. The effects on the vertical tail due to sideslip, roll rate 
and yaw rate are described below. 


- Due to a sideslip /3, the vertical tail produces a side force F v which in turn will result in yawing 
moment about the center of gravity as follows. Note that the side force F v is given by 


l\ — Lj v C[S yClvfiv — 


-q v qS v a v ^-/3 = -q v qS v a v { 1 - 
op 


( 6 . 11 ) 


The yawing moment coefficient produced by this side force is 

AC W „ = C Nv ^ = 

qSb 

Using equation (6.11), we deduce 

Cn v # = 


( 6 . 12 ) 


(6.13) 


where rj v is the ratiobetween the dynamic pressure at the verticaltail and the freestream dynamic 
pressure, l v is the distancefrom the aerodynamic center of the vertical tail to the center of gravity, 
ep is the side wash factor. 

- Due to roll rate p, points along the vertical tail will see an increase in angle of attack of py / V, 
thus we have an incremental yawing moment of, 


dN v = q v (c v dy)a v —l , 


(6.14) 
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where c v is the chord length at a section of the vertical tail located at a distance y from the root. 
Thus, by integrating from 0 to b v (where b v is the height of the vertical tail) and dividing by qSb, 
we obtain the following contribution to yawing moment at the vertical tail due to roll rate, 


A C Nv 


l v P f hv , 
Pv-^r—av / c v ydy 
Sb V Jo 


In terms of the dimensionless roll rate p = pb/2V, we have 


C 


Nv.i 



(6.15) 


(6.16) 


- Due to yaw rate r, the vertical tail will have a change in angle of attack of Aa , u = —rl v /V. An 
incremental side force of 

rl v 

F v = q v qS v a v — (6.17) 

is thereby produced that results in an incremental yawing moment of 

rl v 

A — F v l v = ij v qS v a v —l v (6.18) 

In terms of yawing moment coefficient obtained by dividing equation (6.18) by qSb , we have 

S yl 11 lyf 

A C Nv =-q v -^a v ^~ (6.19) 

Sb V 

or the yaw damping coefficient Cn v , primarily due to the vertical tail is given by 


S v l v l v 

c N vr = -Pv-^-civ— (6.20) 

Sb V 

In terms of the dimensionless yaw rate r — rb/2V, we have 

„ SJ V l v 2V 

c n vT =-Pv-^-a v — — ( 6 . 21 ) 

Sb V b 
or 

C Nv - r =-2q v ^a v l { = -2r, v V v a v l { ( 6 . 22 ) 

Sb b b 

S I 

where Vy — is the vertical tail volume. In Perkins & Hage, additional contribution to the 
total yawing moment due to yaw rate is included for the differential wing drag as follows, 

C D C D l v 

C Np = -^ + C Nlt - =- ^-2q v V v a v y (6.23) 

4 4 b 


• Propulsion Contribution N propu i sion : As in the longitudinal case, the propulsionsystem will contribute 
to the overall yawing moment due to unbalanced thrusts from left and right engines (e.g an engine 
failure). Also the normal forces exerted at the propeller disc due to sideslip will produce additional 
yawing moment. In summary, we can write 


A C Nn 


= A C Nn 


+ A Cv 


propulsion p • 


(6.24) 
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where 


AC 


Npropulsion < 


( T r ight Tjeft )Vp 

qSb 


where y p is the distance between the engine and the fuselage centerlines, and 


a r _ _if Sprop l p 3 Cn p 

Npropulsion p ~ n prop el> t 1 


Sb dp 


(6.25) 


(6.26) 


where N P rop is the number ofpropellers, l p is the x-distance from the propeller to the c.g. (Ligure4.7), 
the coefficient ( by symmetry consideration) can be obtained in a si mi lar fashion using Ligure 


4.8. 


• Rudder Contribution N ru dder : As seen in Section 6.3, rudder surface provides an effective way to 
control the yawing moment of the airplane. The incremental contribution to the yawing moment is 
governed by equation (6.79) from which we obtain the yawing moment coefficient due to rudder 8 r as 
given in equation (6.80). A major consideration in sizing the rudder is in the case of engine failure. 
The rudder must have enough authority to trim out the imbalance in yawing moment — Ty p due to a 
failure, for example, of the left engine. Then, since T = D in trim, 


— Dy p — qSbr) v Vva v r8 r = 0 
or 

§ - Dy P - Cr>y P 
qSbq v Vva v r q v bVva v r 

• Aileron Contribution N a u eron : As the ailerons arc deflected asymmetrically (and in equal amount) to 
produce a roll about the x-axis, drag produced at the downward deflected aileron surface is higher 
than that generated at the upward deflected aileron. In this way, a yawing moment is created by the 
aileron surfaces when the airplane is in a roll maneuver. Lor example, in a positive roll maneuver, 
the left aileron is down and the right aileron is up; thus a negative yawing moment is produced that is 
adverse to the turn coordination. In some cases, to remove the adverse yaw effect, one can increase the 
drag of the upward deflected surface by introducing a greater deflection angle at this surface. Other 
design concepts may be used to generate higher drags (e.g a Lrise aileron sticks out into the flow when 
deflected upward). 

• Spoiler Contribution N spo Her : Aspoileris anaerodynamic device, placedon the upperwing surface, to 
generate drag; thus it is a speed control effector. The spoilers are sometimes also used for roll control. 
When deflected (upward), it causes the flow to separate on the upper surface and thereby resulting in a 
loss of lift. In roll control, deflecting the right spoiler upward will result in a positive rolling moment. 
At the same time, the drag on the right spoiler would generate a desired positive yawing moment 
(hence there is no adverse yaw effect when the spoiler is used to roll the airplane). The difficulty in 
using spoilers as control effectors for roll stabilization is due to the fact its aerodynamic behaviour is 
highly nonlinear, and besides they are less effective than the aileron control surfaces since they are 
located near the wing root. In general, a loss in lift will be accompanied by a loss in altitude, and thus 
may be undesireable. 


(6.27) 

(6.28) 
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6.1.2 Contributions to the Rolling Moment 

Therollingmotionisgenerallyaffectedbythemotionvariablesinyaw r.sidclip ft androll p. Thecomponents 
that contribute mostly to the rolling moment are the wing, the vertical tail, the ailerons (located on the wing) 
and the rudder. 

The airplane rolling moment about the center of gravity can be written as 


L — L w ing + L 


fuselage 


+ Lr 


+ L 


aileron “1“ ^rudder 


(6.29) 


In terms of the moment coefficients, 


C I. qSb ^wing ~ L fuselage ~t - L ven i ca ii a ii ~ t - C'L a j/ eron ~ t~ C Ladder 


(6.30) 


Wing Contribution L,: The rolling moment produced at the wing is developed primarily from 
perturbed motions in sideslip ft, roll rate p and yaw rate r. 


- Due to sideslip, the rolling moment is primarily obtained from the wing dihedral (depicted by 
the dihedral angle T > 0 above the horizontal plane). The rate of change of rolling moment with 
sideslip, Cf f , = dCij’dft, is important to the handling qualities of an airplane. A small negative 
value of Cip is desireable. Too much dihedral will make the airplane hard to fly. When the 
airplane is in a positive sideslip, the right wing will see an increase in angle of attack of 


Act = ftTV/V = ftY 


(6.31) 


The opposite change in a occurs over the left wing. This results in a differential increment in 
rolling moment, 

dL w = —2q{cdy)a w ftYy = —2 qa w ftYcydy (6.32) 


C Lm = -2a w ftY ^ Cydy (6.33) 

Sb 

For a linearly tapered wing, we have 


Cl w 



1 +2A 
1 + X 


)Yft 


(6.34) 


or 


Cl^ 


a w 1 4 2 a 

T i + k 


(6.35) 


Another effect due to sideslip is derived from a swept-back configured wing. Figure 6.2 shows a 
swept wing in a positive sideslip. The velocity normal to the right leading edge is V cos(A — ft) 
andontheleftwing V cos(A+ft). Let C/„ bcthcscction hftcorrespondingtothenormalvelocity 
V cos(A — ft) (or V cos(A + ft)) and normal chord ccos A, 


Differential lift on the right wing — dL r = q cos 2 (A — ft)c cos ACi n ds (6.36) 

and 

Differential lift on the left wing = dL r = q cos 2 (A + ft)c cos A Ci n ds (6.37) 

The differential rolling moment is simply. 



CHAPTER 6. LATERAL STATIC STABILITY AND CONTROL 



Ligure 6.2: Effect of Sweepback on Total Lift and Rolling Moment to Sideslip 


dL V! = (dL l — dL R)y = g[cos 2 (A + fi) — cos 2 (A — / 6)]ccos A Ci n yds 


(6.38) 


where v = .s' cos A or dy = cos Ads. Integrating from 0 to b/2 we obtain, 

r-fc/2 


L w - qCi n [cos 2 (A + ft) - cos 2 (A - ft)] / cydy 

J o 


Note that the incremental lift for a swept wing is 


1 


dL = -p(Vcos A)~Ci n cdy 

Hence, integrating over the entire wing span, the total lift for a swept wing is given by 

rb/2 

L = 2q cos 2 A C/„ / cdy = qSCi n cos 2 A 

Jo 

or the wing Cl and the normal section C/„ are related by 

C L = C, n cos 2 A 

Then the rolling moment coefficient becomes 

C Lw = [cos 2 (A + j3) — cos 2 (A - /?)] /° Cydy 


s 2 A 


Sb 


(6.39) 


(6.40) 


(6.41) 


(6.42) 


(6.43) 
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If we differentiate with respect to ft and evaluate the derivative at ft = 0. we obtain 


C cydy 


C Lw b = -4 Cl tan A 


Again for a linear tapered wing, we have 


Cl w . b = 


l+2X 
3(1 + A.) 


Sb 


Cl tan A 


(6.44) 


(6.45) 


Generally, we have 

C Lw fl = -/(A, X)Cl tan A (6.46) 

where /(A, A.) is an empirically derived function of aspect ratio A and taper ratio X. It should 
be noted that wing placement on the fuselage combined with the cross-flow over the fuselage in 
sideslip (Figure 6.3) introduces additional factors in the rolling moment due to sideslip, i.e. the 
rolling coefficient Cl w 

* High wing: A Cl w ^ = —0.00016 /deg, 

* Mid-wing: A Cl w/3 =0, 

* Lowwing: A Cl wP —+0.00016 /deg, 


More lift 



Figure 6.3: Effect of Wing Placement on the Rolling Moment to Sideslip 


- Due to roll rate, the resulting effect is related to damping in roll. As the airplane rolls, a section 
on the right wing located at a distance y from the centerline will experience an increase in angle 
of attack of, 

py 

Aa = — (6.47) 

V 

Neglecting induced effects (i.e. tilting of the lift vector), the associated incremental rolling 
moment is 

py 

dL w = -2q(cdy)a w —y (6.48) 

By integrating from 0 to b/2 and using the nondimensional variable .r = y/(b/ 2), we have 


qSb 


a w A pb r l 
2 zV Jo 


C 2 

(~)x~dx 

b 


(6.49) 
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where A = b 2 /S is the wing aspect ratio. Lor a linearly tapered wing, we have 


Cl w 


ci w 1 T 3A 
~L2 1 + A 


P 


(6.50) 


or 


Cl„ 


ciyu 1 A 3A 

12 1 + A 


(6.51) 


where the dimensionless roll rate is p = pb/2V. 

The aileron control can be used to produce constant roll rate in steady-state. It is obtained from 
the following equation, 

C L J a + C L -p = 0 (6.52) 


or 



(6.53) 


- Due to yaw rate, the left wing will see a higher velocity than the right wing which is retracting 
away from the forward motion. Assuming that the wing is operating at a constant Cl and the 
section lift C/ is constant and equals to Cl- Then a differential rolling moment is produced from 
the imbalance in dynamic pressure from the two sides; namely 


clL w 


1 

2 


p[{V +ry) 2 - (V - ry) 2 ]C L cdyy 


(6.54) 


clL w = 2pVrC Lcy 2 dy 

Integrating from 0 to Zr/2, and assuming a straight wing with no taper, we obtain 

rb/2 

L w = 2pVrC LC / y 2 dy 
J o 


or 


pcC L rVb 3 
12 


^ qSbrC L 


Therefore, 


or 


Cl w = 


C L _ 

—r 
3 


dC Lw = Cl 
dr 3 


Lor a linearly tapered wing, we derive 


Cl w , f 


C l 1 + 3A 
6 1 + A 


(6.55) 

(6.56) 

(6.57) 

(6.58) 

(6.59) 

(6.60) 


• Luselage Contribution L fuselage • I n general, there is no contribution of the fuselage to the rolling 
moment, i.e. L f use ia g e = 0. 


• Vertical Tail Contribution L ver n C aitaU '■ 
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- Due to yaw rate, the angle of sideslip is decreased by rl v /V. Thus an increment in side force at 
the tail is 

F v = ii v qS v a v r -y (6.61) 

And the resulting rolling moment is 

rl v rb ly ly 

A L v = FyZv = f] v qS v a v —Zv — 2q v qS v a v — —z v = 2q v qS v a v — z v r (6.62) 
V 2 V b b 

where z v is the distance of the aerodynamic center of the vertical tail to the axis of rotation 
(x-axis) and r = i/y is the dimensionless yaw rate variable. 

Dividing equation (6.62) by qSb , we have 

Cl vF =2q v ^^a v ^ = 2q v V v ay : ^ (6.63) 

Sb b b 

- Due to sideslip, the vertical tail produces a side force F v as given in equation (6.11). The rolling 
moment coefficient produced by this side force is 

„ n FyZv qS v Cly(l — Cp)fiZv (LA\ 

A C L . = C L ^ =- = -*. -—- (6.64) 


From which we deduce 


Cl V ' P = - Gs) 


where the variables q v and ep arc as exactly those defined for equation (6.11). 

• Aileron Contribution L a u eron : Aileroncontrols areeffectiveinthe generationofrolling momentdueto 
its location from the axis of rotation (i.e. x-axis). As the right aileron is deflected, there is an increase 
in sectional lift per unit span produced on the right side. It is given by 

clL = qca w dyr8 aR (6.66) 

where r istheaileroneffectiveness (SeeFigure9-15ofPerkins&Hage). Thisresultsinanincremental 
change in the rolling moment as follows, 

dL = —qca w dyx8 aR y (6.67) 

Combining with the contribution from the left aileron, we have 

dL = -qca w dyr(S aR - S aL )y (6.68) 


Integrating over the spanwise length of the aileron, we obtain 


F — qXlw^-i^aR ^ai) 


ryi 

/ cydy 
Jyi 


In dimensionless form, where we define x = y/(b/ 2), the above equation (6.69) becomes 

L s 1 f x 2 c 

Cl = —— = —~a w z8 a A / -xdx (6.70) 

qSb 4 7xi b 

where A = b 2 /S is the wing aspect ratio and 8 a = 8 llR — 8 ai . Again for a simple linearly tapered wing, 
we obtain 

3(xt — Xj) — 2(1 — A.)(x| — Xj) 

Ci. — — a w T - ± - 1 --- 

Sa w 12(1 + A) 


(6.71) 
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• Rudder Contribution L n ,dder ■ As in equation (6.79) for the yawing moment due to rudder, the rudder 
when deflected will also produce a rolling moment, 

A L Sr = ZvVvqSvdv^Sr (6.72) 

or, the rolling moment coefficient with respect to S r is given by 

Cl Sk - Vv—^j— r a v (6.73) 


6.2 Directional Stability (Weathercock Stability) 


Ligure 6.4 shows an airplane at a positive sideslip angle ( J > > 0 or v > 0. The positive yawing moment N is 
defined according to the right-hand rule as shown. The non-dimensional yawing moment coefficient is given 
by 


C N 

N \ P V 2 Sb 


(6.74) 


The change of yawing moment with respect to sideslip ( J > is defined as 


C Nfs = 


dC N 

dp 


(6.75) 


This quantity has the same significance as the coefficient Cm„. However, for stability we see that a positive 
yawing moment would be required to bring ft back to zero. Therefore, we expect Cy, to be positive for 
directional static stability. 

Notethatwhentheairplanehasapositivesideslipthevelocityvectorisnolongerintheplaneofsynmietry. 
Thereexistsayawingmomentproducedbythefuselageandbythe sideforceontheverticaltail. Theairplane 
is directional stable if > 0. Usually the yawing moment due to the fuselage is destabilizing but its effect 
is small compared to the stabilizing moment contributed by the vertical tail. However, in most situation, the 
vertical tail is not sized by any consideration of static stability. Instead the minimum tail size is determined 
by controllability requirements in the event of an asymmetric engine failure or according to flying quality 
requirements. 

If the aerodynamic center of the vertical tail is located a distance of /,, behind the center of gravity. Then 

N = r) v qS v l v a v ( 1 - ep)fi (6.76) 


or 

Cn = = Vv~a v ( 1 - € P )P (6.77) 

qSb S b 

and thus 

Cn p = rj t V v a v (l - ep) (6.78) 

where Vy — S v l v /Sb is the vertical tail volume, ep is the sidewash factor (difficult to estimate), a v is the 
lift-curve slope for the vertical tail. 
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Figure 6.4: Airplane with a Positive Sideslip 
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6.3 Directional Control 

Effective control of yawing moment is provided by the rudder, a movable surface hinged to the vertical 
stabilizer. Incremental yawing moment created by the rudder is 


AN = —l v AL v (6.79) 

where A L v = rj t qS v a v r8 r and the rate of change of Cn with respect to 8 r is given by 


C % = -rj v Vva v x (6.80) 

where r is the effective factor which depends on the ratio of c r /c v (Lor example, see Ligure 5-33 of Perkins 
& Hage). 


6.4 Roll Stability 


Rollmomentisgeneratedbyasynmietricdeflectionof therightandleftailerons. Whenan aiiplaneisinitially 
perturbed in roll, and without the use of the aileron controls, there is no physical mechanism to provide a 
restoring moment. Thus in general Cl 0 is always zero and the conceptof static stability does not exist inroll. 
Or we can also say that the airplane simply possesses neutral static stability in roll. Recall that the rolling 
moment coefficient is defined as 



(6.81) 


Dynamic stability in roll motion is governed by the roll damping term in the stability derivative C/ ;j or 
the non-dimensional stability derivative Cl .. This term is always negative thereby providing positive roll 
damping. 


6.5 Roll Control 

Roll control is provided primarily by the asymmetric deflection of the left and right aileron surfaces. Effec¬ 
tiveness oftheaileroncontrol isdeterminedby therollingmomentcoefficient C/.„„ derivedinequation(6.71). 
A small contribution in roll control can be derived from the rudder control as defined by the term C l s given 
in equation (6.73). In some airplane, additional roll control may be derived from asymmetric deflection of 
the spoilers. 



Chapter 7 

Review of Rigid Body Dynamics 


In general a deformable body of finite dimensions may be regarded as being composed of an infinite number 
of particles, thus the system possesses an infinite number of degrees of freedom. This case would apply to a 
deformable airplane configuration if we take into consideration structural flexibility. However, in this course 
we consider the airplane to be a rigid body with a given mass and moments of inertia. It should be noted 
that for a rigid body, the system undergoes no deformation and should possess only 6 degrees of freedom, 
namely 3 translations and 3 rotations. 

To describe completely the motion of a rigid body, it is convenient to use: 

• 3 translations of a certain point in the rigid body and 

• 3 rotations about that point. 

A system of axes attached to the body arc called body axes. As shown in Figure 7.1, the motion of the body 
can be described by 

1. Translation of the origin O’ of the body axes and 

2. Rotation of the axes with respect to the inertial space. 

It should be noted that velocity of any point P in the rigid body is given by 


Vp = Vo' + a>xrp (7.1) 

Similarly, acceleration of a point P in the rigid body is given by 

ap = ao' + «xrf> + wx(®x rp) (7.2) 

Todevelopthe dynamical equationsfor a rigid body, we need firsttodefine itsl incar and angularmomentum. 
Consider Figure 7.1 where OXY Z is the inertial reference axes and O'xyz corresponds to a set of axes 
attached to the rigid body. Note that the origin O' does not necessarily coincide with the body center of mass 
C. We note the following, 

• Mass of the rigid body m is given by 

m — / dm (7.3) 

J Body 
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Figure 7.1: Motion of a Rigid Body 
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H 0 r = mrc x V o' + / r p x (co x Yp)dm 

J Body 

From here on, we conveniently locate the origin O' to be at the center of mass C, then rc 
equation (7.11) simplifies to the following 


H c = / r P x (co x Yp)dm 
J Body 


I Body 

From vector algebra, we have the following identity: r x (co x r) = a >(r • >’) — r(co • r). 
equation (7.12) becomes 


H c = f 
Je 

In Cartesian coordinates, where 


Body 


<y(rp • rp)dm 


f 

J Body 


r p(a> • rp)dm 


and 


we deduce 


and 


rp = xi + yj + zk 
co = co x i + co v j + cu-k 


rp • rp = v 2 + y 2 + z 2 


co • rp = co x x + co y y + co z z 

Substituting into equation (7.13) we obtain 


H c = / (x 2 + y 2 + z 2 )(oo x i + coyj + cu-k) - ( co x x + co y y + co z z)(x i + yj + zk) 

JBody * 


dm 


He = f Body { [(y 2 + z 2 )co x - xyco y - xzco z \ i + [(x 2 + z 2 )co y - yxco x - yzco-\ j+ 
[(.r 2 + y 2 )co z - zxco x - zyoOy\ k} dm 

Letting He = H x i + H y j + H- k, then 


Hy = 


H y = 


H- = 


/ (y 2 + z 2 )co x - xy co y - xzco z \dm = I xx co x - I^coy - I xz co z 

J Body 

/ (x 2 + z 2 )co v - xy co x - yzco z \dm = -I yx co x + l yy co y - I yz co z 

J Body L 

/ (x 2 + y 2 )co z — zxco x — zycoy \dm — -l zx co x - I zy co y + 

7 Body 


We can define an inertia matrix I to be of the following form, 


1 = 


Lv Ixy Lr 
— lyx lyy ~Iyz 

^zx ~~Izy hz 


(7.11) 
0 and 

(7.12) 
Then 

(7.13) 

(7.14) 

(7.15) 

(7.16) 

(7.17) 

(7.18) 

(7.19) 

(7.20) 

(7.21) 

(7.22) 

(7.23) 
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Notice that the matrix I is a symmetric, positive definite (i.e nonsingular) matrix whose elements have 
units of (ML 2 ). The angular momentum equation in (7.13) now can be re-written as 


H c = Ioj (7.24) 

where the vector co is the angular velocity vector of the rigid body with components (co x , co y , a> z ). The 
inertia matrix I is a constant (i.e., not time-varying) matrix since it is defined in the body-axis O'xyz ■ 
Now we can proceed to the development of the equations of motion for a rigid body from Newton’s 
laws. 


7.1 Force Equations 


We have from Newton’s laws, 


or 


dp 

I 

dt 

(7.25) 

(mVc) = F 

(7.26) 


where p = mX c is the linear momentum of the body, m is the total body mass, Vc is the velocity of the 
center of mass. Since the rigid body has an angular velocity co, then equation (7.26) becomes 


,./v, 


f dXc. 
= m V{— ) 

oxyz at 


+ co x V c ] = F 


O'xyz 


(7.27) 


Let the velocity of the center of mass Vc = ui + v j + ink and the angular velocity of the rigid body be 
co = pi + q\ + rk where oo x = p, co y = q and co- = r. Then the force equations of motion of a rigid body 
airplane are given by 


• Along the body O'x direction: 


m(ii — rv + qw) — F x = X force component 


(7.28) 


• Along the body O'y direction: 

m(v — pw + ru ) = F y = Y force component (7.29) 

• Along the body O'z direction: 

m(w — qu + pv) = F z = Z force component (7.30) 


The force components X, Y and Z on the right-hand side of the above equations are due to gravitional force, 
aerodynamic forces and propulsion forces. We will examine these in the next section. Let’s now proceed to 
the equations of motion for the rotational degrees of freedom. 
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7.2 Moment Equations 


With Newton’s laws applying to the angular momentum, we have 




= <^E) 

oxyz at 


+ co x H c = M 


O'xyz 


Using equation (7.24), equation (7.31) becomes 

Ico + co x loo = M 


(7.31) 


(7.32) 


Using the above definitions for u> and the inertia matrix I, equation (7.32) can be written in the following 
form, 

• About the body O'x direction: 


IXX P - (lyy - Izz)<l r ~ 1 yz ~ r ~) - 4 * (>'' + P c l) - hy (4 ~ >P) = M x = L ( 7 . 33 ) 

• About the body O'y direction: 

lyy 4 - (hz - r xx )rp - Izx (r 2 - P 2 ) - Ixy (P + qr) - I yz (r - pq) = M y = M ( 7 . 34 ) 


• About the body O'z direction: 

h,z r - aXX - lyy )pq - Ixy (p 2 ~ 7U “ Iyz (4 + rp) - I~ x (/) - qr) = M z = N (7.35) 


The moment components L , M and N on the right-hand side of the above equations are due to aerodynamic 
forces and propulsion forces. We will examine these in the next section. Note that there is no contribution 
from the gravitational force since these moments are taken about the center of gravity. 


7.3 Euler’s Angles 

The angular velocity components o> x (or p), a> y (or q) and co z (or r) about the body axes x, y and z cannot 
be integrated to obtain the corresponding angular displacements about these axes. In other words, the 
orientation of the rigid body in space is not known until we describe the three rotational degrees of freedom 
in terms of a set of independent coordinates. Of course, such a set is not necessarily unique. One useful 
set of angular displacements called Euler’s angles obtained through successive rotations about three (not 
necessarily orthogonal) axes as follows. 

Formostairplanedynamics,westartwithasetofinertialaxes OXY Z andperformthefollowingrotations 
in a particular order (Figure 7.2), 

1. Rotation about the Z-axis (i.e yaw) through an angle \fr => (x ], y i, zi), 

2. Rotation about the yi -axis (i.e pitch) through an angle 9 => (x 2 , }’ 2 , z 2 ), 

3. Rotation about the V 2 -axis (i.e roll) through an angle 0 ^ (.V 3 , ya, z 3 ). 
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Equations (7.40)-(7.42) arc solved (i.e integrated) to determine the orientation of the vehicle (cp, 9, xp) from 
the body angular rates ( p, q, r) derived from equations (7.33)-(7.35) when we know the externally applied 
moment components (L, M, N ) 

We can use the above transformations given in equations (7.36)-(7.38) to express the gravitational force 
into the body components as follows. 


(7.43) 


The “flat-earth” model is represented by having the gravitational force always pointed along the vector k. 
Solving for the components of the gravitational force along the vehicle body axes, we obtain 


Fx,gravity — nig S1TI 9 


(7.44) 


Fy,gra vity — mg cos 9 sin cp (7.45) 

and 

F z ,gra vity = mg cos 9 cos cp (7.46) 

These components constitute respectively parts of the force components in the right-hand side of equations 
(7.28)-(7.30). In general, the three successive rotations in equations (7.36)-(7.38) yield the relationship 
between the coordinates in the two reference frames following the Euler’s angle definition. 


~*3~ 


cos 9 cos xp 

cos 9 sin xp 

— sin 9 

X 

37 

_Z3_ 

— 

sin (p sin 9 cos xp — cos cp sin xp 
_ cos cpsin 9 cos xp + sin cp sin xp 

sin cp sin 9 sin xp + cos cp cos xp 
cos cp sin 9 sin xp — sin cp cos xp 

sin cp cos 9 
cos cp cos 9 _ 

y 

_z _ 


Similarly, one can express the coordinates (x, y, z) in terms of the coordinates (* 3 , yy, z,y) by reversing the 
above sequence of rotations. Namely, 


• — cp Rotation: 


• —9 Rotation: 


• —xp Rotation: 


* 2 ” 


“1 0 


0 " 

~*3~ 

37 

= 

0 cos cp 


— sin cp 

37 

Z2_ 


_ 0 sin cp 


cos cp 

_ £3 _ 

Xl " 


cos 9 

0 

sin# ” 

"x 2 " 

34 

— 

0 

1 

0 

37 

Z 1 _ 


_ — sin 9 

0 

cos 9 _ 

-Z2 - 

„ 1 

r 



« n n 

r „ n 


X 


cos xp 

— sin xp 

O' 

■*r 

y 

= 

sin xp 

cos xp 

0 

3 ; 2 

z _ 


0 

0 

1 _ 

-Z3- 


This yields the following complete transformation. 


(7.48) 


(7.49) 


(7.50) 


X 


"cos xp cos 9 

— sin xp cos cp + cos xp sin 9sin cp 

sin xp sin cp + cos xp sin 9coscp 

~*3~ 

y 

= 

sin xp cos 9 

cos xp cos cp + sin xp sin 9 sin cp 

— cos xp sin cp + sin xp sin 9 cos cp 

37 

z _ 


— sin 9 

cos 9 sin cp 

cos 9 cos cp 

_ £3 _ 


(7.51) 
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The above transformation can be used to determine the aircraft position in terms of its linear velocity V with 
components (u, v, w) in the body-fixed axis. First we make the observation 


* = xi + yj + zk = wi 3 + uj 3 + wk 3 


(7.52) 


or 


X 


"cos xp cos 9 

— sin xp cos cp + cos xp sin Osin cp 

sin xp sin cp + cos xp sin Ocoscp 

” u “ 

y 

= 

sin xp cos 9 

cos xp cos cp + sin xp sin 9 sin cp 

— cos xp sin cp + sin xp sin 9 cos cp 

V 

_z _ 


— sin 9 

cos 9 sirup 

cos Ocoscp 

_ w _ 


(7.53) 

Expanding the above expression, we have 

x — u cos xp cos 0 + v(— sin xp cos cp + cos xp sin Osin cp ) + w(sin xp sirup + cos xp sinOcoscp ) (7.54) 
y = u sin xp cos 9 + v(cos xp cos cp + sin xp sin0 sirup) + xv(— cos xp sin cp + sin xp sinO cos cp) (7.55) 

and 

z = —u sin 9 + v cos 9 sirup + w cos 9 cos (p (7.56) 

Let’s summarize here the equations governing the motion of a rigid body aircraft. 

• Lineal - momentum equations: 

- Along the body O'x direction: 


1 

U — T X) CJW g sin 9 + (Xaero T" A propulsion ) (7-57) 

m 

- Along the body O'y direction: 


1 

V — pw 111 -\- g sin (p COS 9 + ( Yaero T Ypropulsion ) (7.58) 

m 

- Along the body O'z direction: 


1 

w — <ju pv + g cos 9 cos <p 4 (Z aero + Zp r0 p U i S i 0n ) (7.59) 

m 

• Angular momentum equations: 

do = I ^ ! (—co x loo + M) (7.60) 

• Equations for the vehicle attitude rates: 


"0 


"1 

0 

— sin 9 

-1 

~ P 


9 

_ xp _ 

— 

0 

_0 

cos cp 
— sin cp 

cos 9 sin cp 
cos 9 cos cp _ 


<7 

_ r _ 

(7.61) 


9 = q cos (p — r sin (p 

xp = q sin (psecO + r cos (psecO 

<p — p + q sin cptan 9 + r cos cptan 9 


(7.62) 
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• Equations for Earth-relative velocities: 

- x-distance: 

x = u cos xp cos 9 + v(— sin xp cos </> + cos xp sin Osin cp) + w(sin xp sirup + cos xp sin Ocoscp) 

(7.63) 

- y -distance: 

y = u sin xp cos 6 + u(cos xf/ cos cp + sin x[r sin 6 sin cp) + w(— cos xjs sin (p + sin xp- sin 9 cos </>) 

(7.64) 

- Vertical altitude h = —z: 

h — u sin 6 — v cos 6 sin </> — w cos 0 cos cp (7.65) 


In the above equations, the external forces F and moments M (about the center of gravity) on the right-hand 
side remain tobe determined. They are derived from basic aerodynamic and propulsionforces and moments. 
In the above derivation, we made the following assumptions: 

• Rigid airframe 

• Flat Earth (i.e gravity is always pointing in the vertical k direction) 

• Axes fixed to the body with origin at the center of gravity 

• Earth-fixed reference is treated as inertial reference 
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Chapter 8 

Linearized Equations of Motion 


8.1 Linearized Linear Acceleration Equations 

Major contributions to the forces and moments in a flight vehicle arc coming from the aerodynamics of 
wings, body and tail surfaces. It would be difficult to express these in terms of the vehicle motion variables 
u, v and w. However it is much easier to express them in terms of the vehicle velocity V, angle of attack 
a and angle of sideslip p. As shown in Figure 8.1, we can express the linear velocities (; u, v, w) directly in 
terms of V, a, P through the following relations: 

V cos p coscr 

V sin/I (8.1) 

V cos p sin a 

where V is the aircraft velocity, a is the aircraft angle of attack and p the aircraft sideslip. 




Figure 8.1: Definition of Angle of Attack a and Sideslip p 
We can rewrite the linear equations of motion given in equations (7.28)-(7.30) as follows, 
u — rv — qw — g sin 6 + X/m 

v = pw — ru + g sin (p cos 0 + Y/m (8.2) 

w = qu — pv + g cos 6 cos</> + Z/m 
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The lineal - accelerations u, v and w can be derived in terms of the variables V, p and a by differentiating 
equations (8.1) with respect to time t. We obtain 

u = Vcos acos/3 — Vasin acos/3 — V ficosasinp 

v = Vsin/3 +/3Vcos (8.3) 

w — Vsin acos/3 + Vacosacosfi — Vpsinasin ft 


Or we can re-write equations (8.3) in terms of a linear system of equations, 


"cosa cos P 

— V sin a cos p 

— V cos a sin p~ 

run 


” u " 

sin p 

0 

V cos p 

a 

= 

V 

_ sin a cos p 

V cos a cos p 

— V sin a sin p _ 

ui 


_ w _ 


Solving for V, a and ft, we obtain 


or 


~V~ 


"cos a cos p 

— V sin a cos p 

— V cos a sin p ” 

a 

= 

sin p 

0 

V cos p 

J. 


_ sin a cos p 

V cos a cos p 

— V sin a sin p _ 



"cos a cos p 

sin a 

V cos f3 

sin p 

0 

sin a cos p ~ 

cos a 

V cos ft 

” it “ 

V 

cosa cos / 3 

L V 

COS ft 

~v~ 

sin a sin ft 

V J 

_ w _ 


Substitute equations (8.2) into equations (8.6) and expand into components, 


u 

v 

w 


(8.4) 


(8.5) 


(8.6) 


V = — g (sin 0 cos acos/3 — cosOsin (psinp — cos 0 cos 4>sin acos/3) 

X Y 7 

+ jpjCosacosP + jjjsinp + j^sinacosp 

a n i cos 9cos (hcos a + sin 6sin a 

a = q — pcosatanp — rsin atanp + g - —n - -n - 

v • 7 ( 8 - 7 ) 

A sin a , Z cos a 

m Vcos p ^ m Vcos p 

P — psin a — rcos a + y (sin Qcosasin p + cos Osin 4>cosP 

X Y Z 

—cos 0 cos (psin a sin P) — -^ycosasin p + -j^y-cosp — —ysinasinP 

where a, p, 6 and cp are the angles of attack, sideslip, pitch and roll respectively, X, Y and Z are the external 
forces along the x, y and z-body axes, m is the total airplane mass, g is the gravitational acceleration and V 
is the total velocity. 

From equations (8.7), we can obtain a set of linearized equations in terms of the perturbation variables 
AV, Aa, A p, A p, A cp A r, AX, AT, A Z, A 6 and A cp in a symmetric climb condition with 


• Lineal' velocities: 


• Angular velocities: 

• Force components: 


V — V 0 + AV) a — a Q + Aa, P — p o + Ap, 
P = Po + Ap, q = q„ + Aq, r = r Q + A r, 

X = X 0 + AX, Y = Y a + AY, Z = Z 0 + AZ 


( 8 . 8 ) 

(8.9) 

( 8 . 10 ) 
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• Airplane attitude angles: 


6 = 0 o + A 9, (/) = (/)„+ Ac/), xft = \j/ 0 + A xft (8.11) 

where V a is the constant aircraft trim velocity, a 0 is the trim angle of attack, 6 0 is the trim airplane pitch 
attitude, X a is the trim force component in the v-direction and Z a is the trim force component in the z- 
direction. For a symmetric climb condition, we have ft,, = p„ = q„ = r a = Y„ = (ft, = ft/,, = 0. These trim 
quantities may not be zero for other flight condition (e.g. steady level turn). Note that the variables AV, 
A a, Aft, A p, A q, A r, AX, AY, A Z, A 6 and Ac/) arc perturbation variables about the trim condition. They 
arc always treated as small quantities. 

Substituting equations ( 8 . 8 )-( 8 .1 l)into equations (8.7), we can derive the linearized equations of motion 
governing the perturbed variables AV, Aft and Act. The linearization is done by neglecting the higher-order 
terms (e.g. Aft A a ~ 0, AV Aa ~ 0, etc...) and invoking at the same time small angle approximations (i.e. 
cos Aft ~ 1, sin Aft ~ Aft, etc ...). After some lengthy manipulation, the following equations of motion for 
the perturbation variables A V, Aa and Aft arc derived, 

Ay = —gcos(0 o — a 0 )A6 H--AZ 

Ad = Aq-^sin(6 o -a o )A0- s -^AX + ^^AZ (8.12) 

Aft — sina 0 Ap — cosa„Ar + y-cos0 o A(/) + m | 7 AT 



Figure 8.2: X and Z-Force Components in terms of L, D and T 

Using Figure (8.2), we can further express the external force components X and Z in terms of the lift 
L, drag D and propulsion force T as follows. 


X 

Z 


Lsin a — Dcos a + Tcos aj 
—Lcosa — Dsin a — Tsin aj 


(8.13) 
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Again we can define L = L a + A L, D = D 0 + A D and T = T 0 + AT. The trim conditions are now 
determined from force balance with the gravity force mg. Namely, 

{ X 0 = L 0 sina 0 — D 0 cosa 0 + T 0 cosaj = mgsin 9 a 
Z a = —L 0 cosa 0 — D 0 sina 0 — T 0 sinar ——mgcosOo 


{ AX = {L 0 cosa 0 + D 0 sina 0 )Aa + sina 0 AL — cosa 0 AD + cosaj AT (8 15) 

A Z = (L 0 sina 0 — D 0 cosa 0 ) Aa — cosa„AL — sina 0 AD — sinajAT 

Substituting equations (8.15) into equations (8.12), we express the linearized equations in terms of A L, AD, 
AT and AT, 

a i> tn \\ a TnSin (On + CCt) a 1 A , cos (a 0 + cut) a t 

AV = -gcos(0 o - a 0 )A9 - - -— -Aa--^AD-\ -— —AT 

Ad = Aq — Y~sin (9 () — a„)(A9 — Aa) — To c °s{Qo + Q'o) 

. , °, , , m ° (8.16) 

sin (a 0 + a T ) A t 1 a r 

mV n 1X1 mV n 


Aft = sina 0 Ap — cosa„Ar + -A-cos9 0 A(j) + 

VO 


We also note that 


D 0 + mgsin (9 0 — a 0 ) — T 0 cos(a 0 + ut) = 0 
L„ — mgcos (9 0 — a 0 ) + T 0 sin (a„ + aj) — 0 


Equations (8.16) then reduce to 


AV = —geos (9 0 - a 0 )A9 - ^ Aa - ± AD + cos(a ° n , + ar) AT 

Ad = Aq - jrsin (B 0 - a 0 )A9 - Jfy-Aa --^AL - sm ( ^+ Ut) AT (8.18) 

A/i — sina„Ap — cosa 0 Ar + ^-cos9 0 A(f> + m y AY 

The above equations arc the linearized equations of the linear acceleration equations. To complete these 
equations, one needs to express in details the terms involved in AL, AD, AT and AT as a function of the 

vehicle motion variables and their perturbations. In the next section, we proceed to formulate the linearized 

equations of motion corresponding to the angular velocity components p, q and r. 


8.2 Linearized Angular Acceleration Equations 

According to equations (7.33)-(7.35) and examining the condition related to a steady level climb condition, 
with the angular velocities defined as 

P = Po + A p 

q = q 0 + Aq (8.19) 

r = r Q + A r 

and the moments as 

L = Tq A AT 

M = M„ + AM (8.20) 

N = N 0 + AN 
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Notice that for a steady level climb condition, p Q = q a — r a = 0 and L a = M () = N 0 = 0. Substituting 
equations (8.19) into equations (7.33)-(7.35) and retaining only the first-order terms in A p, A q and Ar, we 
obtain 

• About the O'y direction: 


IyyAq - ( I zz - I xx )ArAp - Izx(Ar 2 - Ap 2 ) - I xy (Ap - AqAr) - I yz (Ar - ApAq) = AM 

or, after neglecting all the high-order terms, 

IyyAq — IxyAp — I yz Ar = AM (8.21) 

Assuming further that the airplane has a symmetry about the O'xz plane then we have I xy — I yz = 0. 
Equation (8.21) now simplifies greatly to 

IyyAq = AM (Pitching equation) (8.22) 

• About the O'x direction: 

Ixx A p - (Iyy - 1 7Z )AqAr - I yz (Aq 2 - A r 2 ) - I xz (Ar + AqAp) - I xy (Aq - ArAp) = AL 
or, after neglecting all the high-order terms, 

I xx A p — I xz Ar = AL (Rolling equation) (8.23) 

Note that in general I xz ^ 0. 

• About the O'z direction: 

Izz Ar - ( I xx - Iyy) ApAq - I xy (A/; 2 - Aq 2 ) - I yz (Aq + ArAp) - I- x (A/3 - AqAr) = AN 
or, after neglecting all the high-order terms, 

Ar — L x Ap — AN (Yawing equation) (8.24) 


In summary, the equations describing the angular accelerations of the vehicle are given by 


Iyy 

0 

0 1 

” Aq ~ 


"AM” 

0 

IXX 

— Ixz 

Ap 

= 

A L 

L o 


Izz J 

_ Ar _ 


.AN _ 


To complete these equations, one needs to express in details the terms involved in AM, A L and AN as a 
function of the vehicle motion variables and their perturbations. The vehicle attitudes arc described using 
the Euler’s angles. In the next section, equations that describe the vehicle attitudes for small perturbation 
angles arc developed. 
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8.3 Linearized Euler’s Angle Equations 

From equations (7.40)-(7.42), and assuming the trim conditions as defined in equations (8.8)-(8.1 l),we have 
• For the pitch angle 6: 

AOcos Acf> + A \fscos(Q 0 + AO)sin Acj) — Aq 
Since cos Acj> ~ 1 and Aij/sin Acj> ~ 0, then 

A 6 = Aq (Pitch angle equation) (8.26) 


• For the yaw (heading) angle i jr. 

Ai jrcos(0 o + A6)cosA(j) — AO sin Acj) — A r 


or 

Further simplification yields 


A ijs(cos6 0 — sin0 o A0) = A r 


1 

Axjc — -A r (Yaw angle equation) 

cos 0 o 


• For the bank angle 0: 

Acj) — Axj/sin (0 O + AO) = Ap 
or 

Acj) — A\jf(sin0 o + cos6„A6) = A p 
Again, ignoring the high-order terms, 

Acj) — sin 0„ Aij; = Ap (Bank angle equation) 

Using equation (8.27), the above equation simplifies to 

Acj) = Ap + tan 6 0 Ar 


(8.27) 


(8.28) 

(8.29) 


8.4 Forces and Moments in terms of their Coefficient Derivatives 


In this section, the lift, drag and propulsion forces are expressed in terms of the motion variables and their 
perturbations; similarly, for moments about the vehicle axes. Results are represented in terms of the vehicle 
stability derivatives. 

Let’s define the following dimensionless motion related variables, 


P = pb/(2V 0 ) 
q = qc/(2V 0 ) 
r = rb /(2V 0 ) 
V=VIV 0 
a = ac/( 2V 0 ) 
} = $b/( 2Vo) 


(8.30) 


Equations describing forces due to aerodynamics and propulsion are given below. 
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8.4.1 Lift Force L 


L = -pV z SC L = L 0 + AL 


where Cl is the total lift coefficient. At trim, we have 


_ i 2 

Lo — Ltrim — 2^°^° SC Lfrim 


(8.31) 


(8.32) 


And the perturbation in lift A L is given by expanding equation (8.31) in terms of A V and A Cl- Namely 

AL = p 0 V 0 SC Llnm AV + l -poV^SAC L (8.33) 

Now we can express A Cl in terms of the vehicle lift coefficient derivatives, 

A C L = C Lm AM+ C Lh Ah + CL a Aa + C Lp Ap 
+C Lp Ap + C Lq Aq + Cl-A f 

(8.34) 

+C L ,Aa + CL- R Ap 

a p 

-\-Cl S i , A<5 e + Cl Ssp A8 sp + Cl So A8 a + Cl S i . A8 r + • ■ ■ 


Thecoefficients Cl, (andlikewise Cl comeaboutfromthefactthattheinfluenceofdownwash(orsidewash) 

a p 

e on the tail angle of attack (or sideslip) is not felt until it has been propagated from the wing (where it was 
generated) to the tail location traveling at a velocity V Q . This means that the downwash (or sidewash) can be 
expressed in terms of the wing angle of attack a (or sideslip ( J >) through a time-delay function e~ TS where 
t = y- Namely, 

6 = + -C wt, V (8.35) 

da 

and 

c = ^e- l > slv °/3 (8.36) 

Noting that the exponential function e~ l,s ^ Vo can be expanded into 


-Its 
e v « 


1-^ + 1 

To 2 



(8.37) 


Thus, the coefficient Cl, is obtained from 


A L t (s) = rjtqS t a t a t (s ) 


or 

AL,(s) = q t qS t a t (a(s) - e(s)) (8.38) 

Substituting equations (8.35) and (8.37) into equation (8.38), we otain 
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A L,(s) = rj t qS t a t {aQ) - e a - 

da 


l t s 1 / l t s 


1 V 0 + 2\V„' + "' 


a CO 


Or re-written in time domain, we have 

■ 9e 

A L t (t) = q,qS t a t I a(t) - e a - — 

da 


l, 1 / l t \ 2 

a(t) - —a(t) + - ( — ) a(t) + ... 
Vo 


2 \V„ 


From the above equation (8.39), we can deduce that 


S t 


Cl,,„ = rit—at 1 - — 


5 


96 \ 


dot J 


Similarly, 


S f fd€ l t 


CLt a V ‘ cClt \ dot V„ 


and for a better approximation one can also include terms involving a, 

„ S t del fl t ^ 2 

L, - a lh S at da 2\V„ 


(8.39) 


(8.40) 


(8.41) 


(8.42) 


or higher-order time derivatives in a, such as etc... 
In terms of the nondimensional vaiiable a, we have 


^ 2V 0 S t de 21 1 

c L- a = Cl u& -= r) t-Ot- - 

c S da c 


(8.43) 


The same procedure could be applied to the calculation of Cl- and terms involving derivatives with 

£ 

respect to [’> and higher time derivatives in sideslip ft. However, in most circumstances, the effects of sideslip 
on lift arc considered insignificant and can be neglected. 

Note that we make use of the following relationship to convert between Mach number M and velocity V 
derivatives, 

dC L _ dC L dM _ 1 dC L 
~dV ~ ~dM~dV~ ~ a~dM 

and a is the speed of sound. 

8.4.2 Drag Force D 


(8.44) 


D = -pV-SC D = D 0 + AD 


where Co is the total drag coefficient. At trim, we have 

1 


Do = Dtfim = - poVft SC D, r 


(8.45) 


(8.46) 


And the perturbation in drag AD is given by expanding equation (8.45) in terms of AV and AC p. Namely 

AD = p 0 V 0 SC Dtrim AV + 1 - Po V 2 SACd 


(8.47) 
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Now we can express A Co in terms of the vehicle drag coefficient derivatives, 

A Cd = C Dm AM + C Dh Ah + C Da Aa + C Dp Ap 
~\~C [)- A p + Cd-A^ + C Of Ar 

(8.48) 

+C D ,Aa + C D - A(5 

a p 

+Co Se A8 e + Co Ssp A8 sp + C DSa A8 a + Co Sr A8 r + ■ ■ • 

8.4.3 Side-Force T 

1 7 

Y = -pV^SCy = Yo + AY (8.49) 

where Cy is the total side-force coefficient. And the perturbation in side force AT is given by expanding 
equation (8.49) in terms of A V and ACy. Namely 

AT - p 0 V 0 SCy trim AV + l - Po V;SACy = \p 0 V^SACy (8.50) 

Sinceintrim, Y trim — 0 or Cy lrim — 0. Now wecanexpress ACy intermsofthevehicleside-force coefficient 
derivatives, 

ACy = Cy M AM + Cy h All + Cy a Aa + C Yf! Aft 
-\-Cy-Ap + Cy-Aq + Cy-Ar 

(8.51) 

+ Cy, Ad + Cy- Ap 

+ Cy Jjp A 8 sp + Cy in A 8 a + Cy Sr A8 r + ■ ■ ■ 

8.4.4 Thrust Force T 

T = l -pV 2 SC T (8.52) 

where Cy is the total thrust coefficient. And the perturbation in thrust AT is given by expanding equation 
(8.52) in terms of A V and ACy. Namely 

AT = p 0 V 0 SC Tlrim AV + X - Po V 2 SAC t (8.53) 

Now we can express ACy in terms of the thrust coefficient derivatives, 

ACy = Cy-AV + Cy a Acx + Cy s ^ Ai8 t j t + Cy, Ad + • • • (8.54) 


Similarly, the moments due to the aerodynamic and propulsion forces are given below. However, for 
simplicity, we assume that the thrust force is assumed to apply at the airplane center of gravity. Hence it has 
no effects on the vehicle pitching, rolling and yawing moments. 
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8.4.5 Pitching Moment M 


1 9 

M = -pV-ScCu (8.55) 

where Cm is the total pitchingmoment coefficient. And the perturbation inpitching moment AM isgiven by 
expanding equation (8.55) in terms of AV and A Cm- Note that for a vehicle in trim, M, n - m = CM„ im = 0. 
Namely 

AM = Po V 0 ScC Mlnm AV+ \p 0 VjScAC M = \p„V^ScAC M (8.56) 

Now we can express A Cm in terms of the vehicle pitching moment coefficient derivatives, 


A Cm = Cm m AM + Cm,, A/j + Cm„ Aa + Cm,s A^i 
+Cm[, A/3 + CMqAq + C M f Ar 

— , (8.57) 

+Cm- a Ad + Cm-- Aft 

a p 

TCm^ A8 e + Cm Ssp A8 sp + Cm So A8 a + Cm S i . A8 r + ■ ■ ■ 

As in the calculation of C/-, we can use the same approach to obtain Cm-. ■ Namely, 

AM, = — AT,/, (8.58) 


Thus 

or 


C Ma = --C Lil 

c 


, de l, 
Cm t, = -r\ t V H a t — — 
3a V„ 


In terms of the nondimensional variable d, we have 


2V 0 de 21, 

Cm- u = Cm & - — —ViVucit- - 

c da c 


(8.59) 


(8.60) 


8.4.6 Yawing Moment N 


1 9 

N = -pV 2 SbC N (8.61) 

where Cn is the total yawing moment coefficient. And the perturbation in yawing moment AN is given by 
expanding equation (8.61) in terms of AV and ACy . Note that for a vehicle in trim, A,„ m = Cy (rim = 0. 
Namely 

AN = p 0 V 0 SbC Ntrim AV + l - Po V^SbAC N = l -p n V^SbAC N (8.62) 

Now we can express A Cn in terms of the vehicle yawing moment coefficient derivatives, 


ACy — Cy M AM + Cn,, Ah + Cy a Aa + Cy^A/3 
+Cy ( ,A/3 + Cy- Aq + Cy, Ar 
+Cy, Aa + Cy. A/i 

+Cy Jc A 8 e + Cy Jjp A 8 sp + Cy Ja A 8 a + Cy Jr A 8 r + ■ ■ ■ 


(8.63) 
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8.4.7 Rolling Moment L 


1 9 

L = -pV 2 SbC L (8.64) 

where Cl is the total rolling moment coefficient. And the perturbation in rolling moment A L is given by 
expanding equation (8.64) in terms of AV and A Cl- Note that for a vehicle in trim, L trim = c Lmm = o. 
Namely 

A L = p 0 V 0 SbC Ltrim AV + l -p 0 VtSb\C L = ^ Po V~SbAC L (8.65) 

Now we can express A Cl in terms of the vehicle rolling moment coefficient derivatives, 


A C L = C Lm AM+ CL,Ah + C La Aa + CL ll AI3 
+C Lj ,Ap + C L - q Aq + Cl-A f 
+Cl, Aa + C L -,A~p 

a p 

+CLj e A8 e + Cl Ssp A8 sp + Cl Sq A8 a + Cl S i . A8 r + • ■ ■ 


( 8 . 66 ) 


In the most general cases, dynamics of a flight vehicle are fully described by 9 highly coupled nonlinear 
ordinary differential equations governing the motion variables {V, a, fi, /;, q, r, 8,4>, xfr}. Linearization 
about a trim condition reduces them to a set of 9 highly coupled (but) linear ordinary differential equations 
for the perturbed variables {A V, A a, A/i, A p, Aq, Ar, A 6, A<p, Ax//}. When further simplification 
can be achieved by taking into considerationthe airplane symmetry and the decoupled effects in aerodynamic 
forces and moments. This usually leads to the separation of the the basic equations of motion of an airplane 
into two distinct sets: one set corresponds to the longitudinal motion for the variables {A V, Aa, Aq, Ad}, 
and the other set to the lateral motion for the variables {Af8, Ap, Ar, Ar/;, Ai//}. 
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Chapter 9 


Linearized Longitudinal Equations of Motion 


In this chapter, we examine the flight dynamics characteristics associated with motion in the longitudinal 
axis. The assumptions made in the analysis arc that the effects of lateral motion on the aerodynamic and 
propulsion forces and moments associated with the lift L, drag D and thrust T forces arc negligeable. Of 
course, in the model linearization, wealso assume that the motion of the vehicle is undergoing small changes 
in the variables A V, Acr, A q and Ad along with small inputs in the controls A 8 e , A8 sp and A 8 t h ■ 

Simple approximation models for the longitudinal dynamics are developed that provide further insights 
intothefrequency separationbetween thephugoidandtheshort-periodmodes. Timeresponsesofthe vehicle 
motion inthe longitudinalaxis arealsodevelopedillustratingthe effectivenessofthecontrols. Criticaldesign 
parameters affecting these responses arc delineated. Flying qualities of the vehicle arc subsequently defined 
in terms of these fundamental response characteristics. 

Following the linearization performed in Chapter 8, the equations of motion governing the longitudinal 
dynamics arc for the motion variables {A V, Acr, Aq, AO). Thus, in general, it is described by a set of 4 
linear ordinary differential equations obtained from equations (8.18), (8.25) and (8.26). Namely, 


AF = -gcos(9 0 - a 0 )A 9 - A a - |AD + cos{a ] n + a <> ] AT 

Ac* = Aq - fain (9 0 - a 0 )A9 - J^-Aa - Jp-AL - sin ( g + AT 

Aq = jLaM 

lyy 

A 9 = Aq 


The variables 9 a , a 0 , V 0 , L„ and D 0 are determined from the trim conditions involving usually the solutions 
of a set of nonlinear algebraic equations. 

Using expressions for L a , A L, D a , AD and AT as derived in equations (8.32), (8.33), (8.46), (8.47) and 
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(8.53), equations (9.1) become 

AV = gCOS (0 o 0 l o )A 6 yfi 2 _PoV 0 SC L(trim)Ao/ jpj J PoV 0 SC Oitrim ) A V + 2 P 0 V 0 S 

A V + C Dh A h + C Da Aa + C D - A q + C D , Ad + C D&e A S e + C Dssp A S sp j 
+ cos(aT + a 0 ] \_ poV 2 s [c T -AV + C Ta Aa + C Tsih A8 th \ 

o 11? 1 r 

A o' = A q y-sin (0 o o' o )A0 niV 2 SC D(trim ) Acy mV {Po Vo SC l (trim ) A V 

+ 2Po^ 2 ^[% £ AV + C L;i A/ J + C La Aa + C L -A^ + C L ,Ad + C Ljf A5 e (9-2) 

+C Lssp A-V]) - sin %+ ao) $p 0 V?S [C T -AV + C Ta Aa + AS*] 

Ag = -i-i p 0 U 0 2 Sc [%^AU + Cm,, A/j + C M „ Aa + C M ^ Ag + C M , Ad 

+CMj e A<S e + Cm &p A^ p j 
A0 = Ag 


Foranaiiplanetrimmedatlevelflight, wehave Ah = 0. Usingthedefinition ofthenondimensionalvariables 
A V, A q and Ad in equations (8.30), we obtain 


" 1 
0 
0 

.0 


PoVoSc ^ 

> + T5r Ct; 

PoVoSc ^ 

41 yy Cm « 

0 


0 

0 

1 

0 



" A V - 



Ad 



Aq 



_ Ad . 



- £s ^T ( 2Cz) + Cd m M~ C T y cos(a r + a 0 )) 
(2 C L + C Lm M + C t . sin (a T + « 0 )) 

Z/1 yy 

0 


{C L - C Da + C Ta cos(a T + a 0 )) 
(Cd + C La + C r „ sin(a r + a 0 )) 

PoV J Sc 


21 x 


~Cm„ 


PqVqSc_^ 

4m 

l _ Po Sc 


4m 


Cl, - 


PoVoSc ^ 

4 / Cm, 


-g cos (00 - a 0 ) 
4 - sin(9„ - ol 0 ) 

y O 

0 

0 


AU 

Aff 

Aq 

L A0 J 


+ 


PoV„ 2 S 

- 2m CDs, 

_ PoVoS r 

2m Ls ‘ 
p 0 V 2 Sc 

^ Ms. 


21 


yy 


0 


PoV 2 S 


Cd, 


2m 

_ Po VpS /-> 

2m Lss p 

HpjL 


21 


yy 


Ssp 


0 


P °2m S Ct '»«* cos(q T + a ») 

~^^T C T Sth sin(a r + a 0 ) 

PoV 2 Sc n 

^ Ms,h 


21 


yy 


0 


' A<5 e 
A 8 sp 
L A 8 th 


(9.3) 
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We can solve for {A V, Ad, A q, A0} as follows, 


-Ay- 


Ad 


A q 


. A9. 



Po Vo Sc _ 

Am Cd « 

1 + r c ^ 


0 O' 
0 0 


0 - po Y‘j ,Sc ~ C M : 1 o 

.yy 

0 0 0 1J 


(lC D + C Dm M - Ct- cos(ar + av)^) 
-§| (2C L + C Lm M + Cr- sin(a r + a,)) 


,P»U 

2m 


PqVqSc 
21 v 


Cm m M 


0 


(C L - C Da + C Ta cos (a T + a 0 )) C D . 

-^f(CD + c L a + C Ta sm(a T + a 0 )) 1 - ^C L - 

PoVjSc, 


21 y 


-C M „ 


PoVoSc r' 

4 / 


o 

p c y-s 

lm~ < ~' Ds ‘ 

_ PoVoS ^ 

2m 

PoVqSc 


1 


PoVqS 
2 m ° Ss p 


-g cos(9 0 - a 0 ) 

-4- sin(0 o - a 0 ) 

0 

0 


av 

Aa 
A q 

L A9 J 


+ 


”277 


-Cm s , 


PoVqS 

2m 

PoVqSc 


”277 


Cl S! 

Cm s , 


P °2m Ct ™ cos ( a T + a o) 
~^T CT m sin(Q!r +a 0 ) 

PoV;Sc , 


”277 


■Cm,, 


0 


0 


0 


A 8 e 
A 8 sp 
L A5 (/! J 


In abbreviated notation, we write the above equation in the following form 


(9.4) 


-Ay- 


~ AV~ 

Ad 


Acr 

A q 

= F 

A q 

. A 6 _ 


. A 9 . 


+ G 


' A<5 e 
A S sp 
L a 8 th J 


(9.5) 


where the matrices F and G can be deduced from the above equation directly. 

Inmostsituations,wehave Cm m = Oandforlevelflight y 0 = 0 o —a o = 0. thcnthccharactcristiccquation 
forthelongitudinalequationisa4 -orderpolynomialin s ofthefollowingform s 4 + Bs 3 + Cs 2 + Ds + E — 0 
where 222 

E = - P °4 n J° C j2 C L + C Lm M + C Tv sin (a T + a 0 ) } gC Ma (9.6) 

Thus, a necessary condition for the quartic characteristic equation to have stable roots is for the coefficient 
E to be positive. In this case, clearly if the aircraft is statically stable then E will be positive; in another 
words, if the aircraft is statically unstable (i.e. Cm„ > 0) then E is negative and we know with certainty that 
the aircraft is dynamically unstable. This situation illustrates the fact that a statically unstable airplane is 
dynamic ally unstable; however when the airplane isstatically stable, we cannot tell whether it isdynamically 
stable until we solve for the roots of the quartic characteristic equation. In most problems, we actually do 
not calculate the characteristic equation explicitly, rather we develop the corresponding system matrix E 
(defined above) and solve for its eigenvalues (They are then exactly the roots of the characteristic equation). 

A MATLAB-mfile toformulatethe longitudinal equationsof motion isgiven below for thedesign model 
described in Section 11.1. 


rho=0.00230990; 
Vo=556.29559; 
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S=60 8; 
c=15.95; 
b=42.8; 

Ixx=28700; 

Iyy=165100; 

Izz=187900; 

Izx=-52 0; 
lxy=0;lyz=0; 

M=0.5; 

CL=0.20709; 

CD=0.01468; 
g=32.17095; 

Weight=45000;m=Weight/g; 

%Lift Coefficient Derivatives 
CLq=-17.2322; 

CLM=7.45058e-6; 

CLalpha=4.8706; 

CLalphadot=17.2322; 

CLelev=0.572957; 

%Drag Coefficient Derivatives 
CDq=0; 

CDM=0; 

CDalpha=0.37257; 

CDalphadot=0; 

CDelev=4,38308e-2; 

%Pitch Moment Coefficient Derivatives 
CMq=3.8953; 

CMM=-7.05586e-6; 

CMalpha=-0.168819; 

CMalphadot=-ll.887; 

CMelev=-0.695281; 

%Thrust Coefficient Derivatives 
CTv=0; 

CTalpha=0; 
alphaT=0; 

%Trim angle of attack 
alphao=0.18105*pi/180; 
thetao=alphao; 

%Matrix A 

A=[1,rho*Vo*S*c*CDalphadot/4/m, 0,0; 

0,l+rho*S*c*CLalphadot/4/m,0,0; 

0,-rho*Vo*S*c A 2 *CMalphadot/4/Iyy,1,0; 

0 , 0 , 0 , 1 ] ; 

Bl=[-rho*Vo*S/2/m*(2*CD+CDM*M-CTv*cos(alphaT+alphao)) ; 
-rho*S/2/m*(2*CL+CLM*M+CTv*sin(alphaT+alphao)) ; 



rho*Vo*S*c*CMM*M/2/Iyy; 0] ; 

B2=[rho*Vo A 2*S/2/m*(CL-CDalpha+CTalpha*cos(alphaT+alphao) 
-rho*Vo*S/2/m*(CD+CLalpha+CTalpha*sin(alphaT+alphao)) 
rho*Vo A 2 *S*c*CMalpha/2/Iyy; 

0 ] ; 

B3=[rho*Vo*S*c*CDq/4/m; 
l-rho*S*c*CLq/4/m; 
rho*Vo*S*c A 2*CMq/4/Iyy; 

i] ; 

B4=[-g*cos(thetao-alphao) ; 

-g/Vo*sin(thetao-alphao); 

0; 

0 ] ; 

B=[B1,B2,B3,B4]; 

C=[-rho*Vo A 2*S*CDelev/2/m; 

-rho*Vo*S*CLelev/2/m; 
rho*Vo A 2 *S*c*CMelev/2/Iyy; 

0 ] ; 

%Longitudinal equations of motion 
F=inv(A)*B; 

G=inv(A)*C; 

%Phugoid and Short-Period modes 
eigx(F); 

%Elevator pulse inputs 
xo=G*pi/180; 
tl=[0:.5:100]; 
u=zeros(tl); 

yl = lsim(F,G,eye(4),zeros(4,1),u,tl,xo) ; 
t2=[0:.01:10]; 
u=zeros(12); 

y2 = lsim(F,G,eye(4),zeros(4,1) ,u,t2,xo) ; 
clg; 

subplot(221); 
plot(tl,yl(:,1)); 
grid; 

xlabel('Time (sec)') 

ylabel('Velocity (ft/sec)') 

subplot (223); 

plot (tl,180*yl(:,4)/pi); 

grid; 

xlabel('Time (sec)') 

ylabel('Pitch Attitude (deg)') 

subplot(222); 

plot(t2, 180*y2(:,2)/pi); 
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grid; 

xlabel('Time (sec)') 

ylabel('Angle of Attack (deg)') 

subplot(224); 

plot (t2,180*y2(:,3)/pi) ; 

grid; 

xlabel('Time (sec)') 

ylabel('Pitch Rate (deg/sec)') 

pause 

%Short-period approximation: Delete the velocity and theta equations 
F sp=F(2 : 3,2 : 3) ; 

Gsp=G(2:3,1); 
damp(Fsp); 

%Elevator pulse inputs 
xo=Gsp*pi/180; 
tl=[0:.01:10]; 
u=zeros (tl); 

yl = lsim(Fsp,Gsp,eye(2),zeros(2,1) , u,11, xo) ; 
clg; 

subplot (211) ; 
plot(tl,yl(:,1)); 
grid; 

xlabel('Time (sec)') 

ylabel('Angle of Attack (deg)') 

subplot(212); 

plot(tl,180*yl(:,2)/pi) ; 

grid; 

xlabel('Time (sec)') 

ylabel('Pitch Rate (deg/sec)') 

pause 

%Phugoid approximation 
freqphg=g/Vo*sqrt(2) 


Running this MATLAB m-Iile generates the following longitudinal models: 


■AV 


~ AV ~ 

Aa 


Aa 

A q 

= F 

A q 

. A e. 


_ Ad . 


+ GA8 e 


where the matrices F and G arc given by 


F = 

-8.1994e-03 
-1.9451e-04 
6.9573e-04 
0 


-2.5708e+01 
-1,2763e+00 
1.0218e+00 
0 


0 

1.0000e+00 
-2.4052e+00 
1.0000e+00 


-3.2171e+01 
0 
0 
0 


(9.7) 
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-6.80 94e+00 
-1.4968e-01 
-1.4061e+01 
0 


Characteristic roots of the longitudinal model are given by the eigenvalues of the system matrix F. This can 
be calculated using the MATLAB function damp(F). 


Eigenvalues 


0.0012693 

0.0012693 

-0.68348 

-3.0037 


0.103 92i 
-0.10392i 
01 
01 


Damping Frequency(rad/sec) 

0.012 0.10392 (Phugoid mode) 

0.012 0.10392 (Phugoid mode) 

1.000 0.68348 (Short-Period mode) 

1.000 3.0037 (Short-Period mode) 


where the damping ratio £ of a complex root (s = o + j <w) is defined as follows 

—a 

s/o 2 + ur 

Notice that — 1 < K < 1- A negative damping ratio means that the motion is dynamically unstable. 

In most longitudinal aircraft motion, we can distinguish two basic modes: the phugoid mode and the 
short-period mode. The phugoid mode is the one that has the longest time constant and is usually lightly 
damped. Responses of the aircraft that are significantly affected by the phugoid mode are the velocity and 
the pitch attitude as seen in Fig. 9.1 . There is very little response seen in the angle of attack A a when 
the aircraft is excited in the phugoid mode. Note that the period of the velocity AV and pitch attitude 
A q responses is roughly equal to the period of the phugoid mode determined from the approximation 
Tphugoid = sflrtVolg = V2jt556. 29559/32.17095 = 76.8(sec) as discussed in Section 9.1. Damping 
of the phugoid mode is predominantly governed by the drag coefficient C o and its derivative C d m in the 
equation for AV. More precisely, the term (2 C& + Cn M M ) or (1Cd + C o v V a ) is the damping factor in the 
speed equation. 

The short-period mode is displayed in the motion of the airplane angle of attack A a and pitch rate A q 
(Fig. 9.1). It has a relatively short time constant (hence the name short-period). In most situation, this mode 
is relatively well damped (if not, then one must provide stabilization of thismodes using feedback control for 
flight safety since the pilot cannot control this mode). Theshort-period mode is usuallyidentified by a pair of 
complex roots, but in some flight conditions it can be in terms of two real roots as seen in the above example 
problem where 53 = —0.68 rad/sec and 54 = —3.0 rad/sec. In some control design problems, one would like 
to obtain a simplified second-order dynamic model for the longitudinal aircraft motion that captures the fast 
motion (i.e. the short-period mode) only. 


9.1 Phugoid-Mode Approximation 

The phugoid frequency c 0 phugoid i n a level-flight trim condition (i.e. y f , = 6 0 — a 0 = 0) can be estimated 
from the speed and pitch attitude equations where we neglect the variation in Ao/ (i.e. A a = 0) and the drag 
effects. Furthermore, we assume CY,- = 0 and C 7 ;, = 0 and for a fixed elevator A 8 e = 0. In this case, the 
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Time (sec) Time (sec) 


Figure 9.1: Longitudinal Aircraft Responses to a 1-deg Elevator Impulsive Input 


motion of the airplane is governed entirely by the exchange of kinetic and potential energies. And we have 


'AV = -gAO 
Aa = 0 = Aq — po ^ L AV 
AO — Aq = ^f^AV 
or 

'Avir o —gi r a v" 

AO - PqSCl Q AO 

L J L m V J L J 

Characteristic equation of the above equation is simply equal to 


(9.8) 

(9.9) 


= s 2 + 


PogSC L 


= ,- + 2^=0 


with characteristicroots at si t 2 = ±./ s/2g / V„ which correspondtosinusoidal motions inthe aircraftvelocity 
and pitch attitude. Thus, the phugoid mode frequency should be roughly equal to o) p h ug oid = \/ 2 g / V a 
(rad/sec). 

Applyingtothegivennumericalexample, we have o> p i nigr ,i t i = 0.0818rad/secwithaperiodof T p h ug oki = 
277 / (o p hugoid = 76.8 sec. 

Another interpretation of the phugoid oscillation is as follows. The motion exhibited in the phugoid 
mode typifies the exchange of potential and kinetic energies of the aircraft. Recall that when the airplane is 
treated as a point mass m, the total energy is given by 



+ ingh = constant 


(9.11) 
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Differentiating this equation with respect to time, we obtain 

mV V + mgh — 0 
or 

h 

V = -g— = -gy 

Recall from (way back!) equation (1.9), we have 

mVy — L — W 

where L is the total lift given by 

1 


and 


L — 2 PoV SCL(trim) 


IV = mg = -p a V~SC L( 

trim ) 


Then for small velocity perturbations A V in V — V () + A V, we obtain 


AV = -gy = 


G P " (V ” + av,2xCi “ \ p - v:;SCl 


(9.12) 

(9.13) 

(9.14) 


AV = —g — (p 0 V 0 AV SC L ) = 
mV 0 


mV} V 2 


1 


2 - Po V;SC L AV 


AV = -2^-AV 
V} 


(9.15) 


Thus, the perturbed velocity AV has a second-order harmonic motion with frequency a> 0 = Vlg / V„ 
(rad/sec). 


9.2 Short-Period Approximation 


A short-period approximation model is developed based on the following observations: 

• There is a significant frequency separation between the phugoid and short-period modes. Usually 
an order of magnitude difference in frequency between the phugoid and short-period modes, e.g. 

t*)phugoid — 9. 1 rad/sec and co s j-iortperiod — 3 rad/sec. 

• Thevelocityoftheaircrafthasnosignificantcomponentsintheshort-periodmode. Inanotherword,the 
velocitycan be assumed nearly constant when theairplane responds toan excitationinthe short-period 
mode. 


Thus, the longitudinal equations of motion can be simplified by simply removing (i.e. deleting the variables 
AV and Ad) in the original equations. That is, we obtain a set of equations involving only Aa and A q. 
Namely, 


A a 


+ O' 

. Aq _ 


PoV„Sc 2 Cm t 1 

4 1 1 

L lyy 


(Co + C La + C Ta sin(a r + «„)) 

P0VgSC 


2/ v 


-C A 


1 _ PpSC y~, 

4 m L i 

PoVqSc_ 2 


4/v 


-C* 


Aa 
. . 


+ 
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PoViSc, 


PqVqSc f 


P° 2 m C T S, h c °s (oct + oe 0 ) 

PljVgSC „ 


For the above problem, we obtain the following short-period approximation model 


"Ad _ —I 
_A q_ _ 1.1 


.2763 1 

0218 -2.4052 


Act (rad) 
A q 


' -0.1497 
.-14.0611 


where Aa, A<$ ( , are in radians and A q in radians/sec. Characteristic roots of the short-period model are given 
below. They arc almost the same as those obtained in the full ( 4 th -order) longitudinal model. 

Eigenvalues Damping Frequency(rad/sec) 

-0.68299 1.000 0.68299 


Figure 9.2 shows the responses of the short-period approximation model to a 1-deg elevator impulse 
input. Note that these responses of Aa and A q match closely those obtained for the full (i.e. 4’ 1 ' -order) 

p 0 ScC L- 

longitudinal model. In general, C/,,, >> Cp , Ct„ ~ 0 and —^— - << 1, then the above short-period 
approximation has the following characteristic equation 


2 , PoV 0 Sc- 

s + ~1I 
z - 1 yy 


Iyy C La \ , PoV L a Sc^ C Ma 

- i C m < ) s + ^rr Ct » r cr 
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Chapter 10 

Linearized Lateral Equations of Motion 


In this chapter, we examine the flight dynamics characteristics associated with motion in the lateral axis. 

The assumptions made in the analysis arc that the effects of longitudinal motion on the aerodynamic and 
propulsion forces and moments associated with the lift L, drag D and thrust T forces are negligeable. 
Hence, the motion in the lateral axis is decoupled from the longitudinal dynamics. Of course, in any model 
linearization, we also assume that the motion of the vehicle is undergoing small changes in the variables A/S, 

Ap, A r and A</> along with small inputs in the controls A 8 a and A<$,-. 

Equationsgovc rn ingthc lateral mot ionofanaircraftarcassociatcdwiththemot ion variables A/S, A/;, A r, Ar/x 

Hence, the lateral dynamic model is described by a set of 4 linear ordinary differential equations as derived 
in equations (8.18), (8.25) and (8.28). Namely, 


A/S = sina^A p — cosa 0 Ar + A- cos6 0 A(f> H- A- AY 

V o Wl V o 


I XX 

bxz 

Ap 


' A L ' 

_ ~lzx 

1 ZZ. - 

. Ar _ 


.AN _ 


A 0 = A p + tan^Ar 


( 10 . 1 ) 


Using equations (8.50), (8.65) and (8.62) for AT, A L and AN, and retaining only the contributions due to 
lateral motion variables, we obtain 


• Sideslip equation for A/3: 


A$ — sinovAp — cosa 0 Ar H- cos6 0 A(f> 

Vo 


1 1 


+ --PoVtS (c Y pAf3 + C Y[ Ap + Cy-Ar + Cy, A~p+Cy Sa A8 a + Cy Sr AS r ) (10.2) 

mV 0 2 ^ p 7 


1 b \ g 

1- Cy- - AjS — sino! 0 Ap — cosa 0 Ar H-cos 9„Acj) 

mV 0 f>2 V 0 ) V 0 

+ ^r\p° V o S {Cy p A/3 + Cy-Ap + Cy- Ar + Cy Sa A8 a + Cy Sr A8 r ) 


Dividing by (1 


rCy, ), we obtain 


(10.3) 


A (3 = 


2m V, 




p 0 V 0 S ( p 0 bS 

-^—Cy^Ap + (sm « 0 + -^-Cy, ) A p 
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PobS \ 2 PoVoS PoVnS 

+ I — cosa 0 + —— C Yf Ar + — cos0 o A0 + — - C Ysa AS fl H— - C Ygr A8 r 

4 m ) V 0 2m 2m 

Angular velocities A p and Ar: 


1 


I XX Ixz 

Ap 

_ — Izx Izz - 

. Ar _ 


2 PoV„Sb 


b 


Cl, A/3 + CL f ,2y~^P + CL ? 2y-Ar + C t ia A8 a + Cl Sk A8 r 

b 

-2V 0 

• Kinematic equation for Ar/;: 


L Cn p Afi + Cff- A p + Cjy- Ar + Cjv Jn A<5 n + C/v Jr AS r J 


A0 — A/; + tan 0 O Ar 


(10.4) 


(10.5) 


( 10 . 6 ) 


Rc-arrang i ngthcabovecquat ions, wcobtainaseto 141 incarordinarydirfcrcntialcquat ions in 


as follows. 


-Ar 


1 - 

0 

0 

0 

Ap 


0 

l\\- 

~lxz 

0 

Ar 


0 

~lxz 

^zz 

0 

_Acp_ 


0 

0 

0 

1_ 


PoVoS 

2m 


C Ye 


PoV 2 Sb 
2 

PoVjSb 


2 Cl p 


C N/t 


0 


since o + ^-C Yp 


PoV Q Sb 2 

4 

Po V n Sb 2 

4 

1 


Cl, 

C Nr , 


■ cos a 0 + 
PoV a Sb 2 

4 

PoV„Sb 2 

4 

tan 0 o 


PobS 

4m 


Cy- 4jr- COS 
y o 


Cu 

c Nf 


PqVqS 

2m 

PoV 2 Sb 

2 

PoV 2 Sb 

2 

0 


Cy, a 

c Lsa 

C N Sa 


PqVqS 

2m 

PoV 2 Sb 

2 

PoV 2 Sb 

2 

0 


Cl s , 

C Ns , 


A S a 
A Sr 


A P 
A p 
Ar 


+ 


{Ay0, Ap, Ar, Acj)} 


(10.7) 


In abbreviated notation, we write the above equation in the following form 


-a r 


-A/3- 

Ap 

= F 

Ap 

Ar 


Ar 

_ A 


_A <f>_ 


+ G 


A 8 a 
A 8 r 


( 10 . 8 ) 


where the matrices F and G can be deduced from the above equation directly. 

Note that we have added an additional equation for the heading variable Ax// as given in equation (8.27). 
Namely, 

1 

At jr - -Ar (10.9) 

cos0 o 

A MATLAB-mfile to formulatethe lateral equationsof motion for the designmodel described inSection 
11.1 are given below. 



Ill 


clg;clear; 
rho=0.00230990; 

Vo=556.29559; 

S=60 8; 
c=15.95; 
b=42.8; 

Ixx=28700; 

Iyy=165100; 

Izz=187900; 

Izx=-520; 
lxy=0;lyz=0; 

M=0.5; 

CL=0.20709; 

CD=0.01468; 
g=32.17095; 

Weight=45000;m=Weight/g; 

%Side Force Coefficient Derivatives 
CYp=0; 

CYr=0; 

CYbeta=-0.97403; 

CYbetadot=0; 

CYrud=-l.5041e-l; 

CYail=-l.1516e-3; 

%Yawing Moment Coefficient Derivatives 
CNp=-3.3 7 21e-2; 

CNr=-4.0471e-l; 

CNbeta=l.2 996e-l; 

CNrud=-6.97 63e-2; 

CNail=2.1917e-3; 

%Rolling Moment Coefficient Derivatives 
CLp=-0.2; 

CLr=0.15099; 

CLbeta=-0.13345; 

CLrud=-2.3859e-3; 

CLail=2.6356e-2; 

%Trim angle of attack 
alphao=0.18105*pi/180; 
thetao=alphao; 

%Matrix A 

A= [l-b*CYbetadot/(2*m*Vo A 2) ,0,0,0 
0,Ixx,-Izx,0 
0,-Izx,Izz,0 
0 , 0 , 0 , 1 ] ; 

Bl=[rho*Vo*S*CYbeta/2/m; 

0.5*rho*Vo A 2*S*b*CLbeta; 
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0.5*rho*Vo A 2 *S*b*CNbeta;0]; 

B2=[sin(alphao)+rho*b*S*CYp/4/m; 

0.25*rho*Vo*S*b A 2*CLp; 

0.25*rho*Vo*S*b A 2*CNp; 

1 ] ; 

B3=[-cos(alphao)+rho*b*S*CYr/4/m; 

0.25*rho*Vo*S*b A 2*CLr; 

0.25*rho*Vo*S*b A 2*CNr; 
tan(thetao)]; 

B4=[g/Vo*cos(thetao) ; 

0; 

0; 

0 ] ; 

B=[B1,B2,B3,B4]; 

C=[rho*Vo*S*CYail/2/m, rho*Vo*S*CYrud/2/m; 

0.5*rho*Vo A 2*S*b*CLail, 0.5*rho*Vo A 2*S*b*CLrud; 

0.5*rho*Vo A 2*S*b*CNail, 0.5*rho*Vo A 2*S*b*CNrud; 

0 , 0 ] ; 

%Lateral equations of motion 
F=inv(A)*B; 

G=inv(A)*C; 

%Add the heading equation 

Fx=[[F,zeros(4,1)];[0,0,1/cos(thetao),0,0]]; 

Gx—[G;[0,0]]; 

%Dutch roll, spiral and roll modes 
eigx(F); 
eigx(Fx); 

%Aileron pulse inputs 

xo=-Gx(:,1)*pi/180; %Switch sign for + right aileron down 
tl=[0:.1:30]; 
u=zeros(tl) ; 

yl=lsim(Fx,Gx(:,1),eye(5),zeros(5,1),u,tl,xo); 
clg; 

subplot (221) ; 

plot(tl, 180*yl(:,1)/pi); 

grid; 

xlabel('Time (sec)') 

ylabel('Sideslip (deg)') 

title('Aileron 1-deg Pulse Input') 

subplot(223); 

plot(tl,180*yl(:,2)/pi); 

grid; 

xlabel('Time (sec)') 

ylabel('Roll Rate (deg/sec)') 

subplot (222); 
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plot(tl, 180*yl(:,3)/pi); 
grid; 

xlabel('Time (sec)') 
ylabel('Yaw Rate (deg/sec)') 
subplot (224); 

plot (tl,180*yl(:,4)/pi,tl,180*yl(:,5)/pi); 
grid; 

xlabel('Time (sec)') 

ylabel('Roll/Heading Angles (deg)') 
pause 

%Rudder pulse inputs 
xo=Gx(:,2)*pi/180; 
tl=[0:.1:30]; 
u=zeros(tl); 

yl = lsim(Fx,Gx(:,2),eye(5),zeros(5,1),u,tl,xo) ; 
clg; 

subplot (221) ; 

plot(tl, 180*yl(:,1)/pi); 

grid; 

xlabel('Time (sec)') 

ylabel('Sideslip (deg)') 

title('Rudder 1-deg Pulse Input') 

subplot(223); 

plot (tl,180*yl(:,2)/pi); 

grid; 

xlabel('Time (sec)') 

ylabel('Roll Rate (deg/sec)') 

subplot (222); 

plot (tl,180*yl(:,3)/pi); 

grid; 

xlabel('Time (sec)') 
ylabel('Yaw Rate (deg/sec)') 
subplot (224); 

plot (tl,180*yl(:,4)/pi,tl,180*yl(:,5)/pi); 
grid; 

xlabel('Time (sec)') 

ylabel('Roll/Heading Angles (deg)') 


Running this MATLAB m-file generates the following linearized lateral dynamic model: 


-Aj8- 


~ Ap~ 

Ap 


Ap 

A r 

= Fp 

A r 

Acp 


A (j) 

_ A xjr _ 


_ Ai fr _ 


+ G P 


A8 a 

A8 r 


where the matrices F p and G p are given by 


( 10 . 10 ) 
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Fp = 

-2.7202e-01 
-4.3366e+01 
6.552 9e+00 
0 
0 

Gp = 

-3.2161e-04 
8.5397e+00 
8.4854e-02 
0 
0 


3.1599e-03 
-2.4923e+00 
-5.7313e-02 
1.0000e+00 
0 

-4,2005e-02 
-7.1067e-01 
-3.4512e+00 
0 
0 


-1.0000e+00 
1.8964e+00 
-7.7588e-01 
3.1599e-03 
1.0000e+00 


5.7830e-02 
0 
0 
0 
0 


0 

0 

0 

0 

0 


Characteristic roots of the lateral model arc given by the eigenvalues of the system matrix F p . This can be 
calculated using the MATLAB function damp(F). 


Eigenvalues 

0.00000e+00 0.OOOOOe+OOi 
-5.71628e-02 0.OOOOOe+OOi 
-3.27498e-01 2.73177e+00i 
-3.27498e-01 -2.73177e+00i 
-2.82802e+00 0.OOOOOe+OOi 


Damping 
1.000 
1.000 
0.119 
0.119 
1.000 


Frequency(rad/sec) 


0.00000e+00 
5.71628e-02 
2.75133e+00 
2.75133e+00 
2.82802e+00 


(Heading mode) 
(Spiral mode) 
(Dutch-roll mode) 
(Dutch-roll mode) 
(Roll mode) 


There arc 4 basic modes associated with the lateral motion: 


• The heading mode corresponds to a root at the origin (s = 0). This mode is simply associated with 
the integral of yaw rate for the heading angle At fr. 

• The spiral mode (s = —0.05716 rad/sec) is a slow mode that is associated with a real root depicting 
predominantly motion in the roll attitude A<p. Its value is significantly affected by the damping in roll 
from the term Ci p (i.e. rolling moment due to roll rate). At some flight condition, this mode may 
even be unstable; since it is a slow mode, the pilot can interact and correct satisfactorily for the spiral 
instability. 

• The Dutch-roll mode {s — —0.327498 ± /2.73177 rad/sec) is an oscillatory mode with significant 
components in the yaw A r and the roll Ac/) variables. This mode did not have adequate damping 
(£ = 0.12), and in general flying qualities may dictate the need of a lateral stability augmentation 
system to improve the Dutch-roll damnping via a yaw-damper feedback control design. 

• The roll mode (s = —2.82802 rad/sec) is usually associated with a real root which is located far to the 
left, i.e. very stable. The motion is predominantly in roll rate A p and settles down quite quickly. 

ShowninFiguresflO. l)and(10.2)aretimeresponsesinthelateralmotiontoaseparatelyappliedimpulse 
input at the aileron and rudder control surfaces respectively. 



Aileron 1-deg Pulse Input 



10 20 30 

Time (sec) 



10 20 30 

Time (sec) 



Time (sec) 
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Figure 10.1: Lateral Responses to a 1-deg Aileron Impulse Input 
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Time (sec) 
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Figure 10.2: Lateral Responses to a 1-deg Rudder Impulse Input 



Chapter 11 


Flight Vehicle Models 


11.1 Generic F-15 Model Data (Subsonic) 


WING AREA 

608.000 ( 

FT* *2) 

WING SPAN 

42.800 ( 

FT) 

MEAN CHORD 

15.950 ( 

FT) 

VEHICLE WEIGHT 

45000.000 ( 

LB) 

IXX 

28700.000 ( 

SLUG-FT**2) 

IYY 

165100.000 ( 

SLUG-FT**2) 

IZZ 

187900.000 ( 

SLUG-FT**2) 

IXZ 

-520.000 ( 

SLUG-FT**2) 

IXY 

0.000 ( 

SLUG-FT**2) 

IYZ 

0.000 ( 

SLUG-FT**2) 

COEFFICIENT OF LIFT 


0.20709 

COEFFICIENT OF DRAG 


0.01468 

LIFT 

(LBS) 

45001.75037 

DRAG 

(LBS) 

3190.16917 

ALTITUDE 

(FT) 

1000. 

MACH 


0.50000 

VELOCITY 

(FT/SEC) 

556.29559 

EQUIVALENT AIRSPEED 

(KTS) 

324.60706 

SPEED OF SOUND 

(FT/SEC) 

1112.59123 

GRAVITATIONAL ACCEL 

(FT/SEC**2) 

32.17095 

NORMAL ACCELERATION 

(G-S) 

1.00026 

LOAD FACTOR 


1.00013 

DYNAMIC PRESSURE 

(LBS/FT**2) 

357.41711 

DENSITY 

(SLUG/FT**3) = 0.00230990 

WEIGHT (0ALTITUDE) 

(LBS) 

44995.73786 

BETA 

(DEG) 

0.00000 

ALPHA 

(DEG) 

0.18105 

PHI 

(DEG) 

0.00000 

THETA 

(DEG) 

0.18105 

ALTITUDE RATE 

(FT/SEC) 

0.00000 

GAMMA 

(DEG) 

0.00000 

ROLL RATE 

(DEG/SEC) 

0.00000 

PITCH RATE 

(DEG/SEC) 

0.00000 
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YAW RATE 
THRUST 

CONTROL VARIABLES 

ELEVATOR 

THROTTLE 

SPEED BRAKE 

RUDDER 

AILERON 

DIFFERENTIAL TAIL 


CHAPTER 11. FLIGHT VEHICLE MODELS 

(DEG/SEC) = 0.00000 

(LBS) = 3188.85962 


0.05994 
0.06643 
0 . 00000 
0 . 00000 
0 . 00000 
0 . 00000 



11.1. GENERIC F-15 MODEL DATA(SUBSONIC) 


119 


NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES 


SIDE 




DRAG 

LIFT 

FORCE 

ZERO COEFFICIENTS 

1.08760D-02 

1.57360D-01 

2.12070D-18 

ROLL RATE 


0.00000D+00 

0.00000D+00 

0.00000D+00 

PITCH RATE 


0.00000D+00 

-1.72322D + 01 

0.00000D+00 

YAW RATE 


0.00000D+00 

0.00000D+00 

0.00000D+00 

MACH NUMBER 


0.00000D+00 

7.45058D-06 

0.00000D+00 

ALPHA 

(RAD) 

3.7257 0D-01 

4.87 0 60D + 00 

0.00000D+00 

BETA 

(RAD) 

0.00000D+00 

0.00000D + 00 

-9.74030D-01 

ALTITUDE 

(FT) 

0.00000D+00 

0.00000D+00 

0.00000D+00 

ALPHA DOT 


0.00000D+00 

1.72322D + 01 

0.00000D+00 

BETA DOT 


0.00000D+00 

0.00000D+00 

0.00000D+00 

ELEVATOR 


4.38308D-02 

5.72 957D-01 

0.00000D+00 

SPEED BRAKE 


6.4 934 6D-02 

3.74 913D-02 

0.00000D+00 

RUDDER 


0.00000D+00 

0.00000D+00 

-1.50410D-01 

AILERON 


0.00000D+00 

0.00000D+00 

-1.15160D-03 

DIFFERENTIAL 

TAIL 

0.00000D+00 

0.00000D+00 

-7.93150D-02 
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CHAPTER 11. FLIGHT VEHICLE MODELS 


NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES 


ROLLING 

PITCHING 

YAWING 

MOMENT 

MOMENT 

MOMENT 


ZERO COEFFICIENTS 

-2.75542D-20 

4.22040D-02 

-6.80788D-20 

ROLL RATE 


-2.00000D-01 

0.00000D+00 

-3.37210D-02 

PITCH RATE 


0.00000D+00 

3.89530D + 00 

0.00000D+00 

YAW RATE 


1.50990D-01 

0.00000D+00 

-4.04710D-01 

MACH NUMBER 


2.59775D-13 

-7.05586D-06 

2.32475D-12 

ALPHA 

(RAD) 

0.00000D + 00 

-1.68819D-01 

0.00000D+00 

BETA 

(RAD) 

-1.33450D-01 

0.00000D+00 

1.29960D-01 

ALTITUDE 

(FT) 

0.00000D+00 

0.00000D + 00 

0.00000D+00 

ALPHA DOT 


0.00000D+00 

-1.18870D+01 

0.00000D+00 

BETA DOT 


0.00000D+00 

0.00000D + 00 

0.00000D+00 

ELEVATOR 


0.00000D+00 

-6.95281D-01 

0.00000D+00 

SPEED BRAKE 


0.00000D+00 

-4.17500D-01 

0.00000D+00 

RUDDER 


-2.38590D-03 

0.00000D + 00 

-6.97630D-02 

AILERON 


2.63560D-02 

0.00000D + 00 

2.19170D-03 

DIFFERENTIAL 

TAIL 

4.01070D-02 

0.00000D+00 

3.05310D-02 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE 



11.2. GENERIC F-15 MODEL DATA (SUPERSONIC) 
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11.2 Generic F-15 Model Data (Supersonic) 


WING 

AREA 

608.000 

(FT* *2) 

WING 

SPAN 

42.800 

(FT) 

MEAN 

CHORD 

15.950 

(FT) 

VEHICLE WEIGHT 

45000.000 

(LB) 

IXX 


28700.000 

(SLUG-FT**2) 

IYY 


165100.000 

(SLUG-FT**2) 

IZZ 


187900.000 

(SLUG-FT**2) 

IXZ 


-520.000 

(SLUG-FT* *2) 

IXY 


0.000 

(SLUG-FT**2) 

IYZ 


0.000 

(SLUG-FT**2) 


COEFFICIENT OF LIFT 


= 

0.05540 

COEFFICIENT OF DRAG 


= 

0.00308 

LIFT 

(LBS) 

= 

44992.04685 

DRAG 

(LBS) 

= 

2498.81524 

ALTITUDE 

(FT) 

= 

20000 . 

MACH 


= 

1.40000 

VELOCITY 

(FT/SEC) 

= 

1451.69421 

EQUIVALENT AIRSPEED 

(KTS) 

= 

627.54694 

SPEED OF SOUND 

(FT/SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC**2) 

= 

32.11294 

NORMAL ACCELERATION 

(G-S) 

= 

0.99780 

LOAD FACTOR 


= 

1.00172 

DYNAMIC PRESSURE 

(LBS/FT**2) 

= 

1335.83203 

DENSITY 

(SLUG/FT**3) 

= 

0.00126774 

WEIGHT (0ALTITUDE) 

(LBS) 

= 

44914.60527 

BETA 

(DEG) 

= 

0.00000 

ALPHA 

(DEG) 

= 

-1.65589 

PHI 

(DEG) 

= 

0.00000 

THETA 

(DEG) 

= 

-1.65589 

ALTITUDE RATE 

(FT/SEC) 

= 

0.00000 

GAMMA 

(DEG) 

= 

0.00000 

ROLL RATE 

(DEG/SEC) 

= 

0.00000 

PITCH RATE 

(DEG/SEC) 

= 

0.00000 

YAW RATE 

(DEG/SEC) 

= 

0.00000 

THRUST 

(LBS) 

= 

2499.44214 

CONTROL VARIABLES 

ELEVATOR 


= 

0.06772 

THROTTLE 


= 

0.05207 

SPEED BRAKE 


= 

0.00000 

RUDDER 


= 

0.00000 

AILERON 


= 

0.00000 

DIFFERENTIAL TAIL 


= 

0.00000 
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CHAPTER 11. FLIGHT VEHICLE MODELS 


NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES 


SIDE 




DRAG 

LIFT 

FORCE 

ZERO COEFFICIENTS 

1.08760D-02 

1.57360D-01 

2.12070D-18 

ROLL RATE 

(RAD/SEC) 

0.00000D+00 

0.00000D+00 

0.00000D+00 

PITCH RATE 

(RAD/SEC) 

0.00000D+00 

-1.72320D+01 

0.00000D+00 

YAW RATE 

(RAD/SEC) 

0.00000D+00 

0.00000D+00 

0.00000D+00 

VELOCITY 

(FT/SEC) 

0.00000D+00 

0.00000D + 00 

0.00000D+00 

MACH NUMBER 


0.00000D+00 

0.00000D+00 

0.00000D+00 

ALPHA 

(RAD) 

3.7257 0D-01 

4.87 0 60D + 00 

0.00000D+00 

BETA 

(RAD) 

0.00000D+00 

0.00000D + 00 

-9.74030D-01 

ALTITUDE 

(FT) 

0.00000D+00 

0.00000D+00 

0.00000D+00 

ALPHA DOT 

(RAD/SEC) 

0.00000D+00 

1.72320D+01 

0.00000D+00 

BETA DOT 

(RAD/SEC) 

0.00000D+00 

0.00000D+00 

0.00000D+00 

ELEVATOR 


4.38310D-02 

5.72959D-01 

0.00000D+00 

THROTTLE 


0.00000D+00 

0.00000D+00 

0.00000D+00 

SPEED BRAKE 


6.4 9351D-02 

3.7 4 913D-02 

0.00000D+00 

RUDDER 


0.00000D+00 

0.00000D + 00 

-1.50410D-01 

AILERON 


0.00000D+00 

0.00000D + 00 

-1.15160D-03 

DIFFERENTIAL TAIL 

0.00000D+00 

0.00000D+00 

-7.93150D-02 



11.2. GENERIC F-15 MODEL DATA (SUPERSONIC) 
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NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES 




ROLLING 

PITCHING 

YAWING 



MOMENT 

MOMENT 

MOMENT 

ZERO COEFFICIENTS 

2.77192D-19 

4.22040D-02 

-1.74198D-19 

ROLL RATE 

(RAD/SEC) 

-2.00000D-01 

0.00000D+00 

-3.37210D-02 

PITCH RATE 

(RAD/SEC) 

0.00000D+00 

3.89530D+00 

0.00000D+00 

YAW RATE 

(RAD/SEC) 

1.50990D-01 

0.00000D+00 

-4.04710D-01 

VELOCITY 

(FT/SEC) 

3.42955D-17 

-1.15312D-10 

3.06986D-16 

MACH NUMBER 


3.55618D-14 

-1.19570D-07 

3.18322D-13 

ALPHA 

(RAD) 

0.00000D+00 

-1.68819D-01 

0.00000D+00 

BETA 

(RAD) 

-1.33450D-01 

0.00000D+00 

1.29960D-01 

ALTITUDE 

(FT) 

0.00000D+00 

0.00000D + 00 

0.00000D+00 

ALPHA DOT 

(RAD/SEC) 

0.00000D+00 

-1.18870D+01 

0.00000D+00 

BETA DOT 

(RAD/SEC) 

0.00000D+00 

0.00000D + 00 

0.00000D+00 

ELEVATOR 


0.00000D+00 

-6.95279D-01 

0.00000D+00 

THROTTLE 


0.00000D+00 

0.00000D+00 

0.00000D+00 

SPEED BRAKE 


0.00000D+00 

-4.17500D-01 

0.00000D+00 

RUDDER 


-2.38590D-03 

0.00000D + 00 

-6.97630D-02 

AILERON 


2.63560D-02 

0.00000D + 00 

2.19170D-03 

DIFFERENTIAL TAIL 

4.01070D-02 

0.00000D+00 

3.05310D-02 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE 



